01 洛谷P1962 斐波那契数列

题目背景

大家都知道，斐波那契数列是满足如下性质的一个数列：

• f(1) = 1

• f(2) = 1

• f(n) = f(n-1) + f(n-2) (n ≥ 2 且 n 为整数)

题目描述

请你求出 f(n) mod 1000000007 的值。

输入输出格式

输入格式：

·第 1 行：一个整数 n

输出格式：

第 1 行： f(n) mod 1000000007 的值

输入输出样例

输入样例#1： 

5

输出样例#1： 

5

输入样例#2： 

10

输出样例#2： 

55

说明

对于 60% 的数据： n ≤ 92

对于 100% 的数据： n在long long(INT64)范围内。

n 的规模太大，需用矩阵快速幂。

附代码：

[cpp]  view plain  copy

1. #include<iostream>  
2. #include<algorithm>  
3. #include<cstdio>  
4. #define MAXN 5  
5. #define MOD 1000000007  
6. using namespace std;  
7. struct node{  
8.        long long a[MAXN][MAXN];  
9. }ans,base,s;  
10. inline long long read(){  
11.        long long date=0,w=1;char c=0;  
12.        while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}  
13.        while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}  
14.        return date*w;  
15. }  
16. node operator *(const node &x,const node &y){  
17.      node ret;  
18.      for(int i=1;i<=2;i++)  
19.      for(int j=1;j<=2;j++){  
20.              ret.a[i][j]=0;  
21.              for(int k=1;k<=2;k++){  
22.                      ret.a[i][j]+=x.a[i][k]%MOD*y.a[k][j]%MOD;  
23.                      ret.a[i][j]%=MOD;  
24.                      }  
25.              }  
26.      return ret;  
27. }  
28. void mexp(long long k){  
29.      for(int i=1;i<=2;i++)s.a[i][i]=1;  
30.      base.a[1][2]=base.a[2][1]=base.a[2][2]=1;  
31.      while(k){  
32.               if(k&1)s=s*base;  
33.               base=base*base;  
34.               k>>=1;  
35.               }  
36. }  
37. int main(){  
38.     long long k;  
39.     k=read();  
40.     if(k==1||k==2){  
41.                    printf("1\n");  
42.                    return 0;  
43.                    }  
44.     mexp(k-2);  
45.     ans.a[1][1]=ans.a[1][2]=1;  
46.     ans=ans*s;  
47.     printf("%d\n",ans.a[1][2]);  
48.     return 0;  
49. }