用BigDecimal类实现Fibonacci算法

Fibonacci(N)=Fibonacii(N-1)+Fibonacci(N-2)

其中 Fibonacci(0)=0;Fibonacci(1)=1

用循环或则递归实现Fibonacci算法很简单,这里就不说了,如果要用公式实现的话,需要进行开根号和幂运算,普通的long型号只能精确到小数点之后的16位,这就意味着当N很大的时候由于小数位overflow无法进行精确运算,本人测试当N大于92时候因为数据溢出的原因就与实际结果不一样了,如果小数点后面不是精确16位,而是自己设定的话,那么运算结果会精确到自己的期望值。JAVA为我们提供了BigDecimal这个类,BigDecimal类提供了幂运算,但是没有提供开平方运算,Fibonacci公式中的开平方运算要自己实现。这里主要介绍一下自己是如何实现开平方运算的

private static BigDecimal mysqrt(BigDecimal a, int scale) {

        BigDecimal x = new BigDecimal(Math.sqrt(a.doubleValue()),
MathContext.DECIMAL64);
if (scale < 17)
return x; BigDecimal b2 = new BigDecimal("2");
for (int tempScale = 16; tempScale < scale; tempScale *= 2) {
// x = x - (x * x - a) / (2 * x);
x = x.subtract(x.multiply(x).subtract(a)
.divide(x.multiply(b2), scale, BigDecimal.ROUND_HALF_EVEN));
}
return x; }

这里用到了牛顿线性迭代算法,如果不明白的话可以在百度搜索。返回的结果是小数点后精确到scale位的数字。

现在把所有代码给出

package com.fujitsu.soft.rad.devsemi.fibonacci3;

import java.math.BigDecimal;
import java.math.MathContext; public class Fibonacci { /*
* by l.zhang
*/
public static void main(String[] args) { System.out.println(fibonacci(100)); } /*
* by l.zhang
*/ public static String fibonacci(int i) {// 大きい桁表現すればlong 2動作の問題 if (i < 0 || i > 1000) {
System.out.println("please input a new number ");
return "-1";
} else if (i == 0) { return "0";
} else if (i == 1) {
return "1";
} else {
return siki(i); } } /*
* 公式で計算 i>=2 by l.zhang
*/
private static String siki(int i) { BigDecimal temp1 = mysqrt(new BigDecimal("5.0"), 100);// Math.pow(5.0,
// 0.5)
BigDecimal temp2 = temp1.divide(new BigDecimal("5.0"));// Math.pow(5.0,
// 0.5) / 5.0;
BigDecimal temp3 = new BigDecimal("1.0");
BigDecimal temp4 = temp3.add(temp1);// (1 + Math.pow(5.0, 0.5))
BigDecimal temp5 = temp3.subtract(temp1);// (1 - Math.pow(5.0, 0.5))
BigDecimal temp6 = temp4.divide(new BigDecimal("2.0"));// (1 +
// Math.pow(5.0,
// 0.5)) / 2.0
BigDecimal temp7 = temp5.divide(new BigDecimal("2.0"));// (1 -
// Math.pow(5.0,
// 0.5)) / 2.0
BigDecimal temp8 = temp6.pow(i, new MathContext(100));// Mypow(temp6,
// i);//
// Math.pow((1 +
// Math.pow(5.0, 0.5))
// / 2.0, i) BigDecimal temp9 = temp7.pow(i, new MathContext(100));// Mypow(temp7,
// i);//
// Math.pow((1 -
// Math.pow(5.0, 0.5))
// / 2.0, i);
// double temp2 = Math.pow((1 + Math.pow(5.0, 0.5)) / 2.0, i);
// double temp3 = Math.pow((1 - Math.pow(5.0, 0.5)) / 2.0, i);
BigDecimal temp10 = temp8.subtract(temp9);
BigDecimal result = temp2.multiply(temp10); String str_res = result.toString();
System.out.println(str_res);
int index = str_res.indexOf(".", 0); return str_res.substring(0, index);
// (temp1 * (temp2 - temp3));
} private static BigDecimal mysqrt(BigDecimal a, int scale) { BigDecimal x = new BigDecimal(Math.sqrt(a.doubleValue()),
MathContext.DECIMAL64);
if (scale < 17)
return x; BigDecimal b2 = new BigDecimal("2");
for (int tempScale = 16; tempScale < scale; tempScale *= 2) {
// x = x - (x * x - a) / (2 * x);
x = x.subtract(x.multiply(x).subtract(a)
.divide(x.multiply(b2), scale, BigDecimal.ROUND_HALF_EVEN));
}
return x; } }
上一篇:Oracle经典SQL


下一篇:tomcat端口修改以及jvm启动参数设置