hdu 1757 矩阵快速幂 **

一看正确率这么高，以为是水题可以爽一发，结果是没怎么用过的矩阵快速幂，233

题解链接：点我

 #include<iostream>
 #include<cstring>
 const int N=;
 using namespace std;
 int k,m;
 struct Matrix{
     int map[N][N];
 };

 Matrix matrix;

 void Initiate(){
     for(int i=;i<N;i++){
         scanf("%d",&matrix.map[][i]);
     }
     for(int i=;i<N;i++){
         for(int j=;j<N;j++){
             if(i==(j+))matrix.map[i][j]=;
             else matrix.map[i][j]=;
         }
     }
 }

 //矩阵相乘
 Matrix Mul(Matrix &a,Matrix &b){
     Matrix c;
     for(int i=;i<N;i++){
         for(int j=;j<N;j++){
             c.map[i][j]=;
             for(int k=;k<N;k++){
                 c.map[i][j]+=a.map[i][k]*b.map[k][j];
             }
             c.map[i][j]%=m;
         }
     }
     return c;
 }

 //快速幂
 Matrix Pow(int n){
     Matrix t;
     if(n==)return matrix;
     if(n&)return Mul(matrix,Pow(n-));
     else {
         Matrix temp=Pow(n>>);
         return Mul(temp,temp);
     }
 }

 int main(){
     while(scanf("%d%d",&k,&m)!=EOF){
         Initiate();
         if(k<){
             printf("%d\n",k%m);
             continue;
         }
         Matrix temp=Pow(k-);
         int ans=;
         for(int i=;i<N;i++){
             ans+=temp.map[][i]*(N-i-); //最后要乘上f[9],f[8],...,f[1],f[0];
             ans%=m;
         }
         printf("%d\n",ans);
     }
     return ;
 }