HDU1575 Tr A

解题思路：矩阵快速幂模板题，见代码：

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 const int maxn = ;
 #define mod 9973
 int t, n, k;

 struct MT{
     int m[maxn][maxn];
 };

 MT Mul(MT a, MT b)
 {
     MT res;
     memset(res.m, , sizeof(res.m));
     for(int i = ; i < n; i++)
     {
         for(int j = ; j < n; j++)
         {
             for(int k = ; k < n; k++)
             {
                 res.m[i][j] += a.m[i][k]*b.m[k][j] % mod;
             }
             res.m[i][j] %= mod;
         }
     }
     return res;
 }

 MT Product(MT a, int k)
 {
     MT r;
     memset(r.m, , sizeof(r.m));
     for(int i = ; i < n; i++)
     {
         r.m[i][i] = ;
     }

     while(k)
     {
         if(k & ) r = Mul(r, a);
         k >>= ;
         a = Mul(a, a);
     }
     return r;
 }

 int main()
 {
     MT a;
     scanf("%d", &t);
     while(t--)
     {
         scanf("%d %d", &n, &k);
         for(int i = ; i < n; i++)
         {
             for(int j = ; j < n; j++)
             {
                 scanf("%d", &a.m[i][j]);
             }
         }

         MT b = Product(a, k);
         int ans = ;

         for(int i = ; i < ; i++)
         {
             ans += b.m[i][i]%mod;
         }

         printf("%d\n", ans % mod);
     }
     return ;
 }