关于求解斐波那契数列,这是一道比较经典的题目,本文主要是对斐波那契数列求解方法的小结。
首先,定义Fibonacci数列如下:
方法1:
利用递归求解,这是最容易写出的算法,代码如下:
#include<iostream>
using namespace std;long Fibonacci(int n)
{if (n == 0)return 0;else if (n == 1)return 1;elsereturn Fibonacci(n - 1) + Fibonacci(n-2);
}int main()
{int N;cin >> N;cout << Fibonacci(N) << endl;return 0;
}
该算法的时间复杂度为O(2^N)。为什么呢?因为每一次计算Fibonacci(n)时,都需要计算Fibonacci(n-1)和Fibonacci(n-2)共2次,相当于计算了n个2相乘的次数。
方法2:
不用递归,利用迭代的方法计算(理论上讲,任何递归方法都能利用迭代方法实现),代码如下:
#include<iostream>
using namespace std;long Fibonacci(int n)
{if (n == 0)return 0;else if (n == 1)return 1;int num1 = 0;//表示F(n-2)int num2 = 1;//表示F(n-1)for(int i = 1; i < n; i++){num2 = num1 + num2;num1 = num2 - num1;}return num2;
}int main()
{int N;cin >> N;cout << Fibonacci(N) << endl;return 0;
}
迭代算法的时间复杂度为O(N),此处应注意一点,这段代码中只用了两个变量num1和num2来分别代替f(n-2)和f(n-1)。实际上也可以用一个数组来代替整个斐波那契数列,只不过空间复杂度会增加。
方法3:
利用矩阵乘法的原理,斐波那契的递推公式可以表示成如下矩阵形式,所以
求Fibonacci(n)就可以转化成求矩阵A的n-1次幂问题了。而对于矩阵幂的问题,我们有:
利用分治的算法思想可以考虑如下求解一个数A的幂。
从它的求解过程来看,这种算法的时间复杂度为O(logN),要比前两种效率都要好,尤其是在N比较大的时候更能体现出这种优势。该算法的实现需要事先定义矩阵以及矩阵相关的计算,代码如下:
#include<iostream>
#include<string>
using namespace std;//定义2×2矩阵;
struct Matrix2by2
{//构造函数Matrix2by2(long m_00,long m_01,long m_10,long m_11):m00(m_00),m01(m_01),m10(m_10),m11(m_11){}//数据成员long m00;long m01;long m10;long m11;
};//定义2×2矩阵的乘法运算
Matrix2by2 MatrixMultiply(const Matrix2by2& matrix1,const Matrix2by2& matrix2)
{Matrix2by2 matrix12(1,1,1,0);matrix12.m00 = matrix1.m00 * matrix2.m00 + matrix1.m01 * matrix2.m10;matrix12.m01 = matrix1.m00 * matrix2.m01 + matrix1.m01 * matrix2.m11;matrix12.m10 = matrix1.m10 * matrix2.m00 + matrix1.m11 * matrix2.m10;matrix12.m11 = matrix1.m10 * matrix2.m01 + matrix1.m11 * matrix2.m11;return matrix12;}//定义2×2矩阵的幂运算
Matrix2by2 MatrixPower(unsigned int n)
{Matrix2by2 matrix(1,1,1,0);if(n == 1){matrix = Matrix2by2(1,1,1,0);}else if(n % 2 == 0){matrix = MatrixPower(n / 2);matrix = MatrixMultiply(matrix, matrix);}else if(n % 2 == 1){matrix = MatrixPower((n-1) / 2);matrix = MatrixMultiply(matrix, matrix);matrix = MatrixMultiply(matrix, Matrix2by2(1,1,1,0));}return matrix;
}
//计算Fibnacci的第n项
long Fibonacci(unsigned int n)
{if(n == 0)return 0;if(n == 1)return 1;Matrix2by2 fibMatrix = MatrixPower(n-1);return fibMatrix.m00;}int main()
{unsigned int number;cin>>number;cout<<Fibonacci(number)<<endl;return 0;
}
参考:
1.http://blog.csdn.net/liyuanbhu/article/details/51703018
2.http://www.cnblogs.com/python27/archive/2011/11/25/2261980.html