flyinghearts《编程之美》读书笔记连载（3）

作者：flyinghearts

 

2.9 Fibonacci 序列

 

计算 Fibonacci 序列最直接的方法就是利用递推公式 F(n+2)=F(n+1)+F(n) 。而用通项公式来求解是错误的，用浮点数表示无理数本来就有误差，经过n 次方后，当n 相当大时，误差能足够大到影响浮点数转为整数时的精度，得到的结果根本不准。

 

用矩阵来计算，虽然时间复杂度降到O(lgN) ，但要用到矩阵类，相当麻烦。观察：

 

       F(n+2)=F(n)+F(n-1) ＝2*F(n-1)+F(n-2)=3*F(n-2)+2*F(n-4)

 

用归纳法很容易证明 F(n) = F(k)*F(n+1-k) + F(k-1)*F(n-k) ，利用该递推公式和原递推公式，要计算F(n) ，只要计算F([n/2]) 和F([n/2]+1) ，时间复杂度为 O(lgN) 。如：要计算F(58) ， 由 58 -> 29,30 -> 14,15 -> 7,8 -> 3,4 -> 1,2 可知只要算5 次。可以用一个栈保存要计算的数，实际上，将n 的最高位1 （假设在第k 位）左边的0 去除掉后，第m 次要计算的数就是：第k 位到第k-m+1 位这m 个位组成的值为t(m) ，则第m-1 次组成的值为t(m-1) ，则t(m)=2*t(m-1)+( 第k-m+1 位是否为1) 。若第m-1次计算得到了 f(k)和f(k+1)，则：

 

第k-m+1 位为1： f(k),f(k+1) -> f(2*k+1),f(2*k+2)

 

第k-m+1 位为0： f(k),f(k+1) -> f(2*k),f(2*k+1)

 

具体公式见下面代码。

 

下面是计算F(n) 最后四位数（某道ACM题）的代码。

 

扩展阅读：

《编程之美》豆瓣

《编程之美》互动网

《编程之美》，IT人求职面试必读

 

原贴地址：http://blog.csdn.net/flyinghearts/archive/2010/05/18/5605933.aspx