矩阵快速幂hdu 1575

矩阵快速幂就是求方阵A的n次方，一定得是方阵.复杂度为logn，对角线相加就是矩阵的值
主要是套模板,下面给出我收藏的模板

struct mat{
    int m[maxn][maxn];
}unit;
//矩阵乘法
mat operator * (mat a,mat b)
{
    mat ret;
    ll x;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    {
        x=0;
        for(int k=0;k<n;k++)
            x+=mod((ll)a.m[i][k]*b.m[k][j]);
        ret.m[i][j]=mod(x);
    }
    return ret;
}
void init_unit()
{
    for(int i=0;i<maxn;i++)
        unit.m[i][i]=1;
    return ;
}
mat pow_mat(mat a,ll n)
{
    mat ret=unit;
    while(n)
    {
        if(n&1)
        {
            ret= ret*a;
        }
        n>>=1;
        a=a*a;
    }
    return ret;
}

#include<iostream>
#include<algorithm>
#include<iomanip>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<string>
typedef long long ll;
using namespace std;
#define mod(x) ((x)%9973)
int n;
int t;
int x;
int k;
const int maxn = 11;
struct mat{
    int m[maxn][maxn];
}unit;
//矩阵乘法
mat operator * (mat a,mat b)
{
    mat ret;
    ll x;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    {
        x=0;
        for(int k=0;k<n;k++)
            x+=mod((ll)a.m[i][k]*b.m[k][j]);
        ret.m[i][j]=mod(x);
    }
    return ret;
}
void init_unit()
{
    for(int i=0;i<maxn;i++)
        unit.m[i][i]=1;
    return ;
}
mat pow_mat(mat a,ll n)
{
    mat ret=unit;
    while(n)
    {
        if(n&1)
        {
            ret= ret*a;
        }
        n>>=1;
        a=a*a;
    }
    return ret;
}
int main()
{

    cin>>t;
    init_unit();
    while(t--)
    {mat a;

        cin>>n>>k;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
            {
            cin>>a.m[i][j];
            }
            a = pow_mat(a, k);
            ll ans=0;
            for(int i=0;i<n;i++)
            ans+=a.m[i][i];
        ans%=9973;
        cout<<ans<<endl;
    }
    return 0;
}