原理:
斐波那契数列定义
用矩阵形式来描述递推关系:
以上递推式可以写成矩阵形式:
观察结果
例如F(2) = 1, F(3) = 2...
而实际上f(n)就是矩阵的左上角元素就是斐波那契数列的第n项的值
例题:
1242 斐波那契数列的第N项
基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题
斐波那契数列的定义如下:
F(0) = 0
F(1) = 1
F(n) = F(n - 1) + F(n - 2) (n >= 2)
(1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, ...)
给出n,求F(n),由于结果很大,输出F(n) % 1000000009的结果即可。
Input
输入1个数n(1 <= n <= 10^18)。
Output输出F(n) % 1000000009的结果。
Input示例11
Output示例89
代码实现:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9+9;
const int N=2;
ll i,j,k;
ll tmp[N][N],res[N][N];
void multi(ll a[][N],ll b[][N],int n) {
memset(tmp,0,sizeof(tmp));
for(i=0;i<n;i++){
for(j=0;j<n;j++){
for(k=0;k<n;k++){
tmp[i][j]+=(a[i][k]*b[k][j])%mod;
}
}
}
for(i=0;i<n;i++){
for(j=0;j<n;j++){
a[i][j]=tmp[i][j];
}
}
}
void pow(ll a[][N],ll m,int n){
memset(res,0,sizeof(res));
for(i=0;i<n;i++) res[i][i]=1;
while(m){
if(m&1) //m为奇数
multi(res,a,n);
multi(a,a,n);
m/=2;
}
}
int main()
{
ll m;
ll a[N][N];
int n;
while(scanf("%lld",&m)!=EOF){
n=2;
a[0][0]=1,a[0][1]=1,a[1][0]=1,a[1][1]=0;
pow(a,m-1,n);
printf("%lld\n",res[0][0]);
}
return 0;
}