标题: 黄金连分数

标题: 黄金连分数 
 
 
    黄金分割数0.61803... 是个无理数，这个常数十分重要，在许多工程问题中会出现。有时需要把这个数字求得很精确。 
 
    对于某些精密工程，常数的精度很重要。也许你听说过哈勃太空望远镜，它首次升空后就发现了一处人工加工错误，对那样一个庞然大物，其实只是镜面加工时有比头发丝还细许多倍的一处错误而已，却使它成了“近视眼”!!
 
 
    言归正传，我们如何求得黄金分割数的尽可能精确的值呢？有许多方法。 
 
    比较简单的一种是用连分数： 
 
                  1 
    黄金数 = --------------------- 
                        1 
             1 + ----------------- 
                          1 
                 1 + ------------- 
                            1 
                     1 + --------- 
                          1 + ... 
 
                            
 
    这个连分数计算的“层数”越多，它的值越接近黄金分割数。 
 
    请你利用这一特性，求出黄金分割数的足够精确值，要求四舍五入到小数点后100位。 
 
    小数点后3位的值为：0.618 
    小数点后4位的值为：0.6180 
    小数点后5位的值为：0.61803 
    小数点后7位的值为：0.6180340 
   （注意尾部的0，不能忽略） 
 
你的任务是：写出精确到小数点后100位精度的黄金分割值。 
 
注意：尾数的四舍五入！ 尾数是0也要保留！ 
 
显然答案是一个小数，其小数点后有100位数字，请通过浏览器直接提交该数字。 
注意：不要提交解答过程，或其它辅助说明类的内容。 
 
 

解决方法

1. 黄金分割点是（根5-1）/2，求根5的高精度到100位，再大数运算-1 和 /2。

2.递归+大数

3.斐波那契 加 大数

代码

1. 大数+开方

 api没有开根号方法，手工解决。网上搜到的方法是牛顿逼近法。

方法详解：http://blog.csdn.net/zhangpengyu321/article/details/8946703

[java]  view plain copy

1. import java.math.BigDecimal;  
2. import java.math.MathContext;  
3. import java.math.RoundingMode;  
4.   
5.   
6. public class gNumber{  
7. public static void main(String[] args){  
8. //                  x^2 -2 = 0  
9. BigDecimal x1 = new BigDecimal(1);  
10. BigDecimal x2 = new BigDecimal(0);  
11. BigDecimal dif = new BigDecimal(0);  
12.   
13. BigDecimal precision =x1.divide(new BigDecimal(10).pow(100));  
14. System.out.println(precision);  
15. while(true){  
16. x2 = x1.subtract(x1.pow(2).subtract(new BigDecimal(5)).divide(x1.multiply(new BigDecimal(2)),100, BigDecimal.ROUND_HALF_EVEN));  
17.   
18. if(x1.compareTo(x2)==1){  
19. dif = x1.subtract(x2);  
20. }else{  
21. dif = x2.subtract(x1);  
22. }  
23.   
24. if(dif.compareTo(precision)==-1){  
25. //System.out.println(x1+","+x2);  
26. break;  
27. }  
28.   
29. x1 = x2;  
30.   
31. }  
32.   
33.   
34. System.out.println(x1.subtract(new BigDecimal(1)).divide(new BigDecimal(2)));  
35. }  
36. }   

结果为： 0.6180339887 4989484820 4586834365 6381177203 0917980576
2862135448 6227052604 6281890244 9707207204 1893911375

2. 大数+递归

[java]  view plain copy

1. import java.math.BigDecimal;  
2.   
3. public class Main{  
4.  public static void main(String[] args){  
5.     
6.   int i;  
7.   BigDecimal a=new BigDecimal(100);  
8.   for(i=0;i<500;i++)  
9.    a=BigDecimal.ONE.add(BigDecimal.ONE.divide(a,102,BigDecimal.ROUND_DOWN));  
10.   a = a.subtract(BigDecimal.ONE);  
11. System.out.println(a);  
12. }  
13. }  

a的初始化其实是任意的，当循环次数达到了一定的规模，a对结果误差的影响已经很小了。数学之美大概就在于此吧。