[51NOD]-1242 斐波那契數列的第N項 [矩陣快速冪]
阿新 • • 發佈:2019-02-09
基準時間限制: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<stdio.h>
#include<vector>
using namespace std;
typedef long long LL;
typedef vector<LL> vec;
typedef vector<vec> mat;
const LL MOD = 1000000009;
mat mul(mat &A,mat &B){
mat C(A.size(),vec(B[0].size()));
for(int i=0;i<A.size();++i){
for(int k=0;k<B.size();++k){
for (int j=0;j<B[0].size();++j)
C[i][j]=(C[i][j]+A[i][k]*B[k][j])%MOD;
}
}
return C;
}
mat pow(mat A,LL n){
mat B(A.size(),vec(A.size()));
for(int i=0;i<A.size();++i) B[i][i]=1;
while(n){
if(n&1) B=mul(B,A);
A=mul(A,A);
n>>=1 ;
}
return B;
}
int main()
{
mat A(2,vec(2));
A[0][0]=A[0][1]=A[1][0]=1;
A[1][1]=0;
LL n;
scanf("%lld",&n);
A=pow(A,n);
printf("%lld\n",A[1][0]);
return 0;
}