PKU 3070 Fibonacci

 经典的用矩阵解决的题目.

f(k)表示数列的第k项,那么可以很容易的得到

f(n)            1   1                           f(1)

            =            的n-1次方 *

f(n-1)         1   0                          f(0)

 

-______________________-不会打符号郁闷

 

矩阵求次方就是二分的范畴了....

要保留4位就是对10000求模

具体见代码

 

1. #include <iostream>
2. using namespace std;
3. struct Mat
4. {
5.     int mat[2][2];
6.     Mat operator*(const Mat& mat_b)
7.         {
8.             Mat tmp;
9.             int i,j,k;
10.             for(i=0;i<2;i++)for(j=0;j<2;j++)tmp.mat[i][j]=0;
11.             for(i=0;i<2;i++)for(j=0;j<2;j++)for(k=0;k<2;k++)
12.             {
13.                 tmp.mat[i][j]+=mat[i][k]*mat_b.mat[k][j];
14.                 if(tmp.mat[i][j]>=10000) tmp.mat[i][j]%=10000;
15.             }
16.             return tmp;
17.         }
18.     Mat pow(int k)
19.         {
20.             if(k==1) return (*this);
21.             Mat tmp=(*this)*(*this);
22.             if(k&1)
23.                 return tmp.pow(k/2)*(*this);
24.             else
25.                 return tmp.pow(k/2);
26.         }
27. }A;
28. int main()
29. {
30.     int n;
31.     Mat tmp;
32.     A.mat[0][0]=A.mat[0][1]=A.mat[1][0]=1;A.mat[1][1]=0;
33.     while(cin>>n,n+1)
34.         {
35.             if(!n)
36.             {cout<<0<<endl;continue;}
37.             tmp=A;
38.             tmp=tmp.pow(n-1);
39.             cout<<tmp.mat[0][0]%10000<<endl;
40.         }
41.     return 0;
42. }