1. 程式人生 > >[51NOD]-1242 斐波那契數列的第N項 [矩陣快速冪]

[51NOD]-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<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; }