HDU 1575 矩阵快速幂裸题

题意：中文题 我就不说了吧，。。。
思路：矩阵快速幂

// by SiriusRen
#include <cstdio>
#include <cstring>
using namespace std;
int cases,n,k,mod=9973,ans;
struct matrix{int a[100][100];matrix(){memset(a,0,sizeof(a));}}first,cpy;
matrix mul(matrix &a,matrix &b){
    matrix temp;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            for(int k=0;k<n;k++)
            temp.a[i][j]=(a.a[i][k]*b.a[k][j]+temp.a[i][j])%mod;
    return temp;
}
int main(){
    scanf("%d",&cases);
    while(cases--){
        scanf("%d%d",&n,&k),k--,ans=0;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                scanf("%d",&first.a[i][j]),cpy.a[i][j]=first.a[i][j];
        while(k){
            if(k&1)first=mul(cpy,first);
            cpy=mul(cpy,cpy),k>>=1;
        }
        for(int i=0;i<n;i++)ans=(ans+first.a[i][i])%mod;
        printf("%d\n",ans);
    }
}