Tr A
每組數據的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)兩個數據。接下來有n行,每行有n個數據,每個數據的範圍是[0,9],表示方陣A的內容。
#include<cstdio>
#include<cstring>
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;
for(int i = 0; i < n; i++)
{
for(int j = 0; j < n; j++)
{
res.m[i][j]=0;
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 ans;
memset(ans.m, 0, sizeof(ans.m));
for(int i = 0; i < n; i++) //初始化
{
ans.m[i][i] = 1;
}
while(k>0)
{
if(k % 2 == 1)
ans = Mul(ans, a);
k /= 2;
a = Mul(a, a);
}
return ans;
}
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 sum = 0;
for(int i = 0; i < n; i++)
sum += b.m[i][i]%mod;
printf("%d\n", sum % mod);
}
return 0;
}