首页 > 代码库 > 用BigDecimal类实现Fibonacci算法
用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; } }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。