HDU 2256(矩阵快速幂）

相关专题链接：

http://m.blog.csdn.net/blog/cgl1079743846/10309423

题目分析：之前用库函数 fmod（x,mod），然后就wa，毕竟用double来取膜的确是不可能的

然后参考大神的代码过了。

#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
const int N=2;
const int MOD=1024;
struct Matrix
{
    int mat[N][N];
    Matrix operator *(const Matrix &m)const{
    Matrix tmp;
    for(int i=0;i<N;i++)
      for(int j=0;j<N;j++)
      {
         tmp.mat[i][j]=0;
         for(int k=0;k<N;k++)
         {
             tmp.mat[i][j]=(tmp.mat[i][j]+mat[i][k]*m.mat[k][j])%MOD;
         }
          tmp.mat[i][j] %= MOD;
      }
      return tmp;
    }
};
int pow(Matrix & m,int k)
{
    if(k==1)
    return 9;
    k--;
    Matrix ans;
    memset(ans.mat,0,sizeof(ans.mat));
    for(int i=0;i<N;i++)
    ans.mat[i][i]=1;

    while(k)
    {
        if(k&1)
        ans=ans*m;
        m=m*m;
        k>>=1;
    }
    int x = (ans.mat[0][0]*5+ans.mat[0][1]*2)%MOD;
    return (2*x-1)%MOD;
}
int main()
{
    int n,t;
    Matrix m;

    scanf("%d",&t);
    while(t--)
    {
       scanf("%d",&n);
         m.mat[0][0]=5;m.mat[0][1]=12,m.mat[1][0]=2,m.mat[1][1]=5;
        printf("%d\n",pow(m,n));
    }
    return 0;
}