poj 3070 Fibonacci (矩阵快速幂乘/模板)

题意：给你一个n,输出Fibonacci (n)%10000的结果

思路：裸矩阵快速幂乘，直接套模板

代码：

#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

typedef long long ll;
const int N=,M=,P=;
const int MOD=;
struct  Matrix
{
    ll m[N][N];
};

Matrix  A={,,
          ,};

Matrix  I={,,
          ,};

Matrix multi(Matrix a,Matrix b)
{
    Matrix ans;
    for(int i=;i<N;i++)
    {
       for(int j=;j<M;j++)
       {
          ans.m[i][j]=;
          for(int k=;k<P;k++)
          {
             ans.m[i][j]+=a.m[i][k]*b.m[k][j]%MOD;
          }
          ans.m[i][j]%=MOD;
       }
    }
    return ans;
}

Matrix power(Matrix a,int k)
{
    Matrix ans=I,p=a;
    while(k)
    {
      if(k&)
      {
        ans=multi(ans,p);
      }
      k>>=;
      p=multi(p,p);
    }
    return ans;
}

int main(int argc, char const *argv[])
{
    int n;
    while(scanf("%d",&n)!=-)
    {
       if(n==-) break;
       if(n==)
       {
           cout<<""<<endl;
           continue;
       }
       Matrix ans=power(A,n-);
       printf("%lld\n",ans.m[][] );
    }
    return ;
}