Poj3070(Fibonacci)

这题就是最水的矩阵快速幂.

求斐波纳契数列的第n项~ ~

AC代码

#include <iostream>
#include <cstdio>
#include <string.h>
#include <stdlib.h>
#define maxn 5
#define MOD 10000
using namespace std;
class Matrix{
public:
    int row,col;
    __int64 mapp[maxn][maxn];
    Matrix(int r,int c):row(r),col(c){}
    Matrix(){}
    void ini()
    {
        memset(mapp,0,sizeof(mapp));
        mapp[0][0]=mapp[0][1]=mapp[1][0]=1;
    }
    void unit() //单位矩阵的初始化
    {
        memset(mapp,0,sizeof(mapp));
        for(int i=0;i<row;i++)
         mapp[i][i]=1;
    }
    int result() const{return mapp[0][1]%10000;}
    friend Matrix operator*(const Matrix&,const Matrix&);
    void Pow(int);
};
Matrix operator*(const Matrix& M1,const Matrix& M2)
{
    Matrix M(M1.row,M2.col);
    for(int i=0;i<M1.row;i++)
     for(int j=0;j<M2.col;j++)
     {
        M.mapp[i][j]=0;
        for(int k=0;k<M2.col;k++)
         M.mapp[i][j]=(M.mapp[i][j]+M1.mapp[i][k]*M2.mapp[k][j])%MOD;
     }
    return M;
}
Matrix M(2,2);
void Matrix::Pow(int n)//矩阵快速幂
{
    //cout<<n<<endl;
    Matrix temp(2,2);
    temp.ini();
    while(n)
    {
        if(n&1)
         M=M*temp;
        temp=temp*temp;
        n>>=1;
    }
}
int main()
{
    __int64 n;
    while(scanf("%I64d",&n),n!=-1)
    {
        M.unit();
        if(n==0) cout<<"0"<<endl;
        else if(n==1) cout<<"1"<<endl;
             else
             {
                 Matrix resM(1,2);
                 resM.mapp[0][0]=1,resM.mapp[0][1]=1;
                 M.Pow(n-1);
                 resM=resM*M;
                 cout<<resM.result()<<endl;
             }
    }
    return 0;
}