Fibonacci介紹
Fibonacci定義：
avatar
Fibonacci數(shù)列默認(rèn)第一個(gè)數(shù)為0,第二個(gè)為1,其余后面的數(shù)為前兩個(gè)之和，例如：
0,1,1,2,3,5,8,13,21......
普通計(jì)算第n個(gè)Fibonacci的方式是從第三位開始求前兩位之和，直至第n位。這種計(jì)算方式的時(shí)間復(fù)雜度是O(n)。
下面用二分來計(jì)算Fibonacci
二分矩陣求Fibonacci
項(xiàng)目地址：github/DivideAndConquer/Fibonacci/
首先有下面的數(shù)學(xué)公式
avatar
也就是我們需要計(jì)算矩陣的n次方，這與利用分治計(jì)算冪是相同的道理，即
avatar
代碼如下：
#include <iostream>
#include <cstring>
using namespace std;
#define MAX_SIZE 1000
//二維矩陣類
class Marix {
public:
Marix(){};
Marix(int x1,int x2,int y1,int y2){ value[0][0]=x1;value[0][1]=x2;value[1][0]=y1;value[1][1]=y2;};
int value[2][2];
bool isSetValue;
void Show(){cout<<value[0][0]<<"\t"<<value[0][1]<<endl<<value[1][0]<<"\t"<<value[1][1]<<endl;};
};
//矩陣乘法
Marix Mult(Marix m1,Marix m2){
Marix ans;
ans.value[0][0]=m1.value[0][0]*m2.value[0][0]+m1.value[0][1]*m2.value[1][0];
ans.value[0][1]=m1.value[0][0]*m2.value[0][1]+m1.value[0][1]*m2.value[1][1];
ans.value[1][0]=m1.value[1][0]*m2.value[0][0]+m1.value[1][1]*m2.value[1][0];
ans.value[1][1]=m1.value[1][0]*m2.value[0][1]+m1.value[1][1]*m2.value[1][1];
return ans;
}
//用于保存計(jì)算過的冪的值
Marix powerMap[MAX_SIZE];
//初始化表
void InitMap(Marix a)
{
for(int i=2;i<MAX_SIZE;i++)
{
//memset(powerMap[i].value, 0x0, 4 * sizeof(int));
powerMap[i].isSetValue=false;
}
powerMap[1]=a;
powerMap[1].isSetValue=true;
}
//設(shè)置指數(shù)為n的值
void SetPowerN(Marix a,int n)
{
if(n%2==0)
{
if(!powerMap[n/2].isSetValue)
{
SetPowerN(a,n/2);
}
powerMap[n]=Mult(powerMap[n/2],powerMap[n/2]);
}
else
{
if(!powerMap[n/2].isSetValue)
{
SetPowerN(a,n/2);
}
powerMap[n]=Mult(Mult(powerMap[n/2],powerMap[n/2]),powerMap[1]);
}
}
//獲取指數(shù)n的值，不存在就調(diào)用SetPowerN
long long GetPowerN(Marix a,int n)
{
if(powerMap[n].isSetValue>0)
{
return powerMap[n].value[1][1];
}
SetPowerN(a,n);
return powerMap[n].value[1][1];
}
int main() {
Marix a(1,1,1,0);
int n;
cin>>n;
InitMap(a);
cout<<GetPowerN(a,n)<<endl;
return 0;
}
