题目:定义Fibonacci数列如下:
/ 0 n=0
f(n)= 1 n=1
/ f(n-1)+f(n-2) n>1
输入n,用最快的方法求该数列的第n项。
分析:在很多C语言教科书中讲到递归函数的时候,都会用Fibonacci作为例子。因此很多程序员对这道题的递归解法非常熟悉,看到题目就能写出如下的递归求解的代码。
方法1代码:
但是,教科书上反复用这个题目来讲解递归函数,并不能说明递归解法最适合这道题目。我们以求解f(10)作为例子来分析递归求解的过程。要求得f(10),需要求得f(9)和f(8)。同样,要求得f(9),要先求得f(8)和f(7)……我们用树形结构来表示这种依赖关系
f(10)
/ /
f(9) f(8)
/ / / /
f(8) f(7) f(6) f(5)
/ / / /
f(7) f(6) f(6) f(5)
我们不难发现在这棵树中有很多结点会重复的,而且重复的结点数会随着n的增大而急剧增加。这意味这计算量会随着n的增大而急剧增大。事实上,用递归方法计算的时间复杂度是以n的指数的方式递增的。大家可以求Fibonacci的第100项试试,感受一下这样递归会慢到什么程度。在我的机器上,连续运行了一个多小时也没有出来结果。当你我实现测试的时候,求第40项就能感觉出它的效率不行了。
有没有更加优化的方法呢?
其实改进的方法并不复杂。上述方法之所以慢是因为重复的计算太多,只要避免重复计算就行了。比如我们可以把已经得到的数列中间项保存起来,如果下次需要计算的时候我们先查找一下,如果前面已经计算过了就不用再次计算了。
更简单的办法是从下往上计算,首先根据f(0)和f(1)算出f(2),在根据f(1)和f(2)算出f(3)……依此类推就可以算出第n项了。很容易理解,这种思路的时间复杂度是O(n)。
方法2代码:
下面我们再看看有没有其它的方法:
如果我们知道一个数列的通项公式,使用公式来计算会更加容易,能不能把这个函数的通项公式计算出来呢?
由递推公式
F( n )=F( n-1 )+F( n-2 ),知道F(n)的特征方程为
X的平方 = X +1
它的根为: x1 = (1 + sqrt(5))/2, x2 = x1 = (1 - sqrt(5))/2
所以存在A,B,使得:
F(n)= A*(x1的n次方)+B*(x2的n次方)
代入F(0)=0,F(1)=1,解得 A=sqrt(5)/5,B=-sqrt(5)/5,则:
F(n)=sqrt(5)/5 * (x1的n次方) - sqrt(5)/5 * (x2的n次方),而x1,x2已知,于是可以得到答案
哎,终于写完了这些,由于编辑器的原因,只能这样写,在纸上画一画就很清楚了!
方法3代码:
上面这个方法,我们可以在O(1)的时间内得到F(n),但公式中引入了无理数,所以不能保证结果的精度!
下面再看一种在下损失精度的情况下,时间复杂度是O(logn)的方法,其速度更快,先看一下下面的数学公式
注意到Fibonacci数列是二阶递推数列,所以存在一个2*2的矩阵A,做得:
( F(n), F(n-1) ) = ( F(n-1), F(n-2) ) *A
求解,可得 :A= (1,1,1,0)
(注:{f(n+1), f(n), f(n), f(n-1)}表示一个矩阵。在矩阵中第一行第一列是f(n+1),第一行第二列是f(n),第二行第一列是f(n),第二行第二列是f(n-1)。)
由上面的我们可以知道如下关系:
( F(n), F(n-1) ) = ( F(n-1), F(n-2) ) *A = ( F(n-2), F(n-3) ) *A的平方 = .... .... .. .. .=(F(1), F(2))*A的n-1次方
所以问题就是转化成求矩阵A= (1,1,1,0)的n-1次方!如果简单第从0开始循环,n次方将需要n次运算,并不比前面的方法要快。但我们可以考虑乘方的如下性质:
/ an/2*an/2 n为偶数时
an=
/ a(n-1)/2*a(n-1)/2 n为奇数时
要求得n次方,我们先求得n/2次方,再把n/2的结果平方一下。如果把求n次方的问题看成一个大问题,把求n/2看成一个较小的问题。这种把大问题分解成一个或多个小问题的思路我们称之为分治法。这样求n次方就只需要logn次运算了。
实现这种方式时,首先需要定义一个2×2的矩阵,并且定义好矩阵的乘法以及乘方运算。当这些运算定义好了之后,剩下的事情就变得非常简单
方法4代码:
下面的对上面 四种方法的测试,你就会发现第一种方法与其它三种方法的优劣了:
测试代码:
测试结果如下:
Please enter the nubmer n1, n2, n3, n4 :
30 30 30 30
The Fibonacci result of Fibonacci_Solution1 30 is: 832040
The Fibonacci result of Fibonacci_Solution2 30 is: 832040
The Fibonacci result of Fibonacci_Solution3 30 is: 832040
The Fibonacci result of Fibonacci_Solution4 30 is: 832040
The End!
啊!终于写完了,欢迎大家拍砖!