HDU1575 Tr A

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

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

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

MT Mul(MT a, MT b)
{
    MT res;
    memset(res.m, 0, sizeof(res.m));
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            for(int k = 0; 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, 0, sizeof(r.m));
    for(int i = 0; i < n; i++)
    {
        r.m[i][i] = 1;
    }

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

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

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

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

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