Tr A HDU 1575 (矩阵快速幂)

 #include<iostream>
 #include<vector>
 #include<string>
 #include<cmath>
 #include<algorithm>
 #include<cstdio>
 #include<cstring>
 #include<list>

 using namespace std;

 #define maxn 15
 int n, k;
 struct matrix//定义一个结构体，方便传递值
 {
     int m[maxn][maxn];
 };

 /*
 maxn和mod由全局定义，其中mod根据需要可以省去
 */

 matrix mul(matrix a, matrix b)    //矩阵求积, 矩阵乘法
 {
     matrix ans;
     for(int i = ; i <= n; i++)
     {
         for(int j = ; j <= n; j++)
         {
             ans.m[i][j] = ;
             for(int k = ; k <= n; k++)
             {
                 ans.m[i][j] += (a.m[i][k] * b.m[k][j]) % ;
                 ans.m[i][j] %= ;
             }
         }
     }
     return ans;
 }

 matrix quick_pow(matrix a, int b)    //矩阵快速幂
 {
     matrix ans;
     for(int i = ; i <= n; i++)
     {
         for(int j = ; j <= n; j++)
         {
             if(i == j)
                 ans.m[i][j] = ;
             else
                 ans.m[i][j] = ;//这里要初始化为单位矩阵，类比普通快速幂这里初始化为1
         }
     }
     while(b != )//方法与普通快速幂相同，只有乘法的实现不同
     {
         if(b %  == )
             ans = mul(a, ans);
         a = mul(a, a);
         b /= ;
     }

     return ans;
 }

 int main()
 {
     int T;
     cin >> T;
     while(T--)
     {
         matrix a;

         cin >> n >> k;
         for(int i = ; i<= n; ++i)
             for(int j = ; j <= n; ++j)
                 cin >> a.m[i][j];

         matrix tmp = quick_pow(a, k);
         int ans = ;
         for(int i = ; i<= n; ++i)
             ans += tmp.m[i][i] % ;

         ans %= ;    // 最后这里一定要再次取余！
         cout << ans << endl;
     }

     return ;
 }