矩阵乘法

#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#include <string>
#include <iostream>
#include <stack>
#include <map>
#include <queue>
using namespace std;
#define N 15
#define INF 0x3f3f3f3f
#define MOD 9973
int n;

struct Mat
{
    int mat[N][N];
    void init()
    {
        memset(mat, 0, sizeof(mat));
    }
    void E()
    {
        memset(mat, 0, sizeof(mat)); //构造单位矩阵
        for(int i = 0; i < n; i++)
            mat[i][i] = 1;
    }
};

Mat operator * (Mat a, Mat b) //矩阵乘法
{
    Mat c;
    c.init();
    for(int i = 0; i < n; i++) { //第一个矩阵的行数
        for(int j = 0; j < n; j++) { //第二个矩阵的列数
            for(int k = 0; k < n; k++) { //第一个矩阵的列数
                c.mat[i][j] = c.mat[i][j] + a.mat[i][k] * b.mat[k][j] % MOD;
                c.mat[i][j] %= MOD; //要加上这句话
            }
        }
    }
    return c;
}

Mat operator ^ (Mat a, int k) //矩阵快速幂
{
    Mat p, tmp = a;
    p.E();
    while(k) {
        if(k & 1) p = p * tmp;
        tmp = tmp * tmp;
        k >>= 1;
    }
    return p;
}

int main()
{
    int t;
    scanf("%d", &t);
    while(t--) {
        Mat a;
        int k;
        scanf("%d%d", &n, &k);
        for(int i = 0; i < n; i++)
            for(int j = 0; j < n; j++)
                scanf("%d", &a.mat[i][j]);
        Mat ans = a ^ k;
        int sum = 0;
        for(int i = 0; i < n; i++) {
            sum = sum + ans.mat[i][i] % MOD;
        }
        printf("%d
", sum % MOD);
    }
    return 0;
}