hdoj1575 Tr A（矩阵快速幂）

简单的矩阵快速幂。最后求矩阵的秩。

 #include<iostream>
 #include<cstring>
 using namespace std;
 const int Mod=;
 int n,k;
 struct mat{
     int a[][];
 };
 mat res,c;
 mat mat_mul(mat x,mat y){
     mat ans;
     memset(ans.a,,sizeof(ans.a));
     for (int i=;i<n;i++)
         for (int j=;j<n;j++)
         for (int k=;k<n;k++){
             ans.a[i][j]+=x.a[i][k]*y.a[k][j];
             ans.a[i][j]%=Mod;
         }
      return ans;
 }
 int mat_pow(int n,int k){
     for (int i=;i<n;i++){
         for (int j=;j<n;j++){
             cin >> res.a[i][j];
         }
     }
     c=res;
     k--;
     while (k){
         if (k&) res=mat_mul(res,c);
         c=mat_mul(c,c);
         k>>=;
     }
     int ans=;
     for (int i=;i<n;i++){
         ans=(ans+res.a[i][i])%Mod;
     }
     return ans;
 }
 int main(){
     int t;
     cin >> t;
     while (t--){
         cin >> n >> k;
         cout << mat_pow(n,k) << endl;
     }
     return ;
 }