程序员面试题精选--(05)求Fibonacci数列

题目：定义Fibonacci数列如下：

 

        /  0                      n=0
f(n)=      1                      n=1
        /  f(n-1)+f(n-2)          n>1

 

输入n，用最快的方法求该数列的第n项。

 

分析：在很多C语言教科书中讲到递归函数的时候，都会用Fibonacci作为例子。因此很多程序员对这道题的递归解法非常熟悉，看到题目就能写出如下的递归求解的代码。

 

方法1代码：

 

 unsigned long long Fibonacci_Solution1 ( unsigned int n ) { if ( n == 0 ) return 0; else if ( n == 1 ) return 1; else return Fibonacci_Solution1( n - 1 ) + Fibonacci_Solution1( n - 2 ); }

 

 

 

 

但是，教科书上反复用这个题目来讲解递归函数，并不能说明递归解法最适合这道题目。我们以求解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代码：

unsigned long long Fibonacci_Solution2 ( unsigned int n ) { if ( n == 0 ) return 0; else if ( n == 1 ) return 1; else { unsigned long long fibOne = 0; unsigned long long fibTwo = 1; unsigned long long fibN = 0; for (unsigned int i = 2; i <= n; i ++ ) { fibN = fibOne + fibTwo; fibOne = fibTwo; fibTwo = fibN; } return fibN; } }

 

下面我们再看看有没有其它的方法：

如果我们知道一个数列的通项公式，使用公式来计算会更加容易，能不能把这个函数的通项公式计算出来呢？

由递推公式

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代码：

unsigned long long Fibonacci_Solution3 ( unsigned int n ) { unsigned long long fibN; double sqrtof5 = sqrt( double ( 5 )); fibN = ( sqrtof5/5 )* pow ( double ( ( 1 + sqrtof5 ) / 2 ) , double ( n ) ) - ( sqrtof5/5 )* pow ( double ( ( 1 - sqrtof5 ) / 2 ) , double ( n ) ); return fibN; }

 

上面这个方法，我们可以在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次运算，并不比前面的方法要快。但我们可以考虑乘方的如下性质：

 

 

 

 

 

 

 

 

 

 

 

 

        /  a^n/2*a^n/2                   n为偶数时
aⁿ=
        /  a(^n-1)/2*a^(n-1)/2            n为奇数时

 

要求得n次方，我们先求得n/2次方，再把n/2的结果平方一下。如果把求n次方的问题看成一个大问题，把求n/2看成一个较小的问题。这种把大问题分解成一个或多个小问题的思路我们称之为分治法。这样求n次方就只需要logn次运算了。

实现这种方式时，首先需要定义一个2×2的矩阵，并且定义好矩阵的乘法以及乘方运算。当这些运算定义好了之后，剩下的事情就变得非常简单

 

方法4代码：

/*************************************************************/ // Multiply two matrices // Input: matrix1 - the first matrix // matrix2 - the second matrix //Output: the production of two matrices /*************************************************************/ struct Matrix2by2 { Matrix2by2( long long m00 = 0, long long m01 = 0, long long m10 = 0, long long m11 =0 ) :m_00 ( m00 ), m_01 ( m01 ), m_10 ( m10 ), m_11 ( m11 ) { } long long m_00; long long m_01; long long m_10; long long m_11; }; Matrix2by2 MatrixMulitply ( const Matrix2by2 matrix1, const Matrix2by2 matrix2 ) { return Matrix2by2 ( matrix1.m_00 * matrix2.m_00 + matrix1.m_01 * matrix2.m_10, matrix1.m_00 * matrix2.m_01 + matrix1.m_01 * matrix2.m_11, matrix1.m_10 * matrix2.m_00 + matrix1.m_11 * matrix2.m_10, matrix1.m_10 * matrix2.m_01 + matrix1.m_11 * matrix2.m_11 ); } Matrix2by2 MatrixPower ( unsigned int n ) { assert ( n > 0 ); Matrix2by2 matrix; if ( n == 1 ) matrix = Matrix2by2 ( 1, 1, 1, 0 ); else if ( n % 2 == 0 ) { matrix = MatrixPower( n / 2 ); matrix = MatrixMulitply( matrix, matrix ); } else if ( n % 2 == 1 ) { matrix = MatrixPower( ( n - 1)/ 2 ); matrix = MatrixMulitply( matrix, matrix ); matrix = MatrixMulitply( matrix, Matrix2by2 ( 1, 1, 1, 0 ) ); } return matrix; } unsigned long long Fibonacci_Solution4 ( unsigned int n ) { if ( n == 0 ) return 0; else if ( n == 1 ) return 1; else { Matrix2by2 matrix2powerN = MatrixPower( n-1 ); return matrix2powerN.m_00; } }

 

 

下面的对上面 四种方法的测试，你就会发现第一种方法与其它三种方法的优劣了：

 

测试代码：

int main () { unsigned int n1, n2, n3, n4; cout << "Please enter the nubmer n1, n2, n3, n4 :" <<endl; cin >> n1 >> n2 >> n3 >> n4; cout << "The Fibonacci result of Fibonacci_Solution1 " << n1 << " is: " << Fibonacci_Solution1( n1 ) <<endl <<endl; cout << "The Fibonacci result of Fibonacci_Solution2 " << n2 << " is: " << Fibonacci_Solution2( n2 ) <<endl <<endl; cout << "The Fibonacci result of Fibonacci_Solution3 " << n3 << " is: " << Fibonacci_Solution3( n3 ) <<endl <<endl; cout << "The Fibonacci result of Fibonacci_Solution4 " << n4 << " is: " << Fibonacci_Solution4( n4 ) <<endl <<endl; return 0; }

 

测试结果如下：

 

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！

啊！终于写完了，欢迎大家拍砖！