hdu 1575 （矩阵快速幂）

点击打开链接

#include"stdio.h"
#include"string.h"
#define N 9973
struct node
{
	int a[11][11];
}A,B;
int n,m;
void init()
{
	int i,j;
	memset(A.a,0,sizeof(A.a));
	for(i=1;i<=n;i++)
		A.a[i][i]=1;//单位矩阵
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
			scanf("%d",&B.a[i][j]);
	}

}
node mult(node A,node B)
{
	int i,j,k;
	node C;
	memset(C.a,0,sizeof(C.a));
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
		{
			for(k=1;k<=n;k++)
			{
				C.a[i][j]+=(A.a[i][k]*B.a[k][j])%N;
				C.a[i][j]%=N;
			}
		}
	}
	return C;
}
int main()
{
	int i;
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d",&n,&m);
		init();
		while(m)
		{
			if(m&1)
				A=mult(B,A);
			B=mult(B,B);
			m/=2;
		}
		int ans=0;
		for(i=1;i<=n;i++)
			ans+=A.a[i][i];
		printf("%d\n",ans%N);
	}
	return 0;
}