Tr A （矩阵快速幂（模板））

【题目来源】：https://vjudge.net/problem/HDU-1575
【题意】
中问题意不在解释。
【思路】
很水的一个快速幂题，而且又是给出了一个矩阵，所以想到矩阵快速幂。。
然而，并没有在敲完代码后直接过，因为在代码里我用到了两个矩阵：ans和base，他们俩的含义就和快速幂的一样，但是呢，本题是T组输入，然后我把他们俩用结构体定义在了外面，base是一直在更新（因为是把矩阵直接输入的），但是ans没有清0，所以又浪费了一点时间，，
【代码】

#include<set>#include<map>#include<stack>#include<cmath>#include<queue>#include<cstdio>#include<string>#include<cstring>#include<iostream>#include<limits.h>#include<algorithm>#define mem(a,b) memset(a,b,sizeof(a))using namespace std;const int mod=9973;typedef unsigned long long ll;typedef long long LL;int n,k;struct mat{    int a[14][14];    mat()    {        mem(a,0);        for(int i=1; i<=10; i++)            a[i][i]=1;    }} base;mat operator*(mat s,mat t){    mat r;    mem(r.a,0);    for(int i=1; i<=n; i++)        for(int j=1; j<=n; j++)            for(int p=1; p<=n; p++)            {                r.a[i][j]=r.a[i][j]+s.a[i][p]*t.a[p][j];                if(r.a[i][j]>=mod)                    r.a[i][j]%=mod;            }    return r;}void pow_mat(mat &ans){    while(k)    {        if(k&1)            ans=ans*base;        base=base*base;        k>>=1;    }}int main(){    int T;    scanf("%d",&T);    while(T--)    {        mat ans;        scanf("%d%d",&n,&k);        for(int i=1; i<=n; i++)            for(int j=1; j<=n; j++)                scanf("%d",&base.a[i][j]);        pow_mat(ans);        int sum=0;        for(int i=1; i<=n; i++)        {            sum+=ans.a[i][i];            if(sum>=mod)                sum%=mod;        }        printf("%d\n",sum);    }}