時間限制：1 秒

記憶體限制：128 兆

特殊判題：否

提交：283

解決：168

題目描述：

A為一個方陣，則Tr A表示A的跡（就是主對角線上各項的和），現要求Tr(A^k)%9973。

輸入：

資料的第一行是一個T，表示有T組資料。
每組資料的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)兩個資料。接下來有n行，每行有n個數據，每個資料的範圍是[0,9]，表示方陣A的內容。

輸出：

對應每組資料，輸出Tr(A^k)%9973。

樣例輸入：

2
2 2
1 0
0 1
3 99999999
1 2 3
4 5 6
7 8 9

樣例輸出：

2
2686

間。

#include<iostream>
using namespace std;
#define MOD 9973

class Matrix
{
public:
	int m[11][11];
	int n;//方陣的行數
	int tr;
	Matrix()
	{
		tr = 0;
		for (int i=0; i<11; i++)
		{
			for (int j=0; j<11; j++)
			{
				m[i][j] = 0;
			}
		}
	}
	Matrix(int n)
	{
		tr = 0;
		this->n = n;
		for (int i=0; i<n; i++)
		{
			for (int j=0; j<n; j++)
			{
				m[i][j] = 0;
			}
		}

	}
	void getT(int n)
	{

		tr=0;
		for (int i=0; i<n; i++)
		{
			tr = (tr + m[i][i])%MOD;
		}
	}

	Matrix operator * (const Matrix a) const
	{
		Matrix res;
		res.n = n;
		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] = (res.m[i][j] + (m[i][k] * a.m[k][j]) % MOD) % MOD;
				}
			}
		}
		return res;
	}
};


int main()
{
	freopen("input.in", "r", stdin);
	freopen("output.out","w", stdout);

	int num;
	while (cin >> num && num!=0)
	{
		
		int n, k, i, j, x, y, z;
		for (i=0; i<num; i++)//每組一次
		{

			//輸入部分
			cin >> n >> k;
			Matrix mtx(n);
			for (x=0; x<n; x++)
			{
				for (y=0; y<n; y++)
				{
					cin >> mtx.m[x][y];
				}
			}

			//計算部分
			int tmp = k;
			Matrix res(n);
			for (z=0; z<n; z++)
			{
				res.m[z][z] = 1;
			}
			while (k>0)
			{
				if (k%2 == 1)
				{
					res = res * mtx;
				}
				mtx = mtx*mtx;
				k/=2;
			}
			res.getT(n);
			cout << res.tr << endl;
		}//end of for
	}//end of while

	return 0;
}