斐波那契数列算法小结

　　关于求解斐波那契数列，这是一道比较经典的题目，本文主要是对斐波那契数列求解方法的小结。

首先，定义Fibonacci数列如下：



方法1：

　　利用递归求解，这是最容易写出的算法，代码如下：

#include<iostream>
using namespace std;

long Fibonacci(int n)
{
if (n == 0)
return 0;
else if (n == 1)
return 1;
else
return 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