hdu 1575

#include <stdio.h>
#include <string.h>
const int mod=9973;
#define  N  11

struct matrix
{
	int m[N][N];    
}A,I;

int n;  //全局变量和局部变量出现问题 

matrix matrixmul(matrix a,matrix b)
{
	matrix c;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		{
			c.m[i][j]=0;
			for(int k=0;k<n;k++)
				c.m[i][j]+=(a.m[i][k]*b.m[k][j])%mod;    
			c.m[i][j]=c.m[i][j]%mod;       
		}
		return c;
}

matrix qiuckpow(int k)
{
	matrix ans=I,b=A;
	
	while(k>=1)
	{
		if(k&1)
            ans=matrixmul(ans,b);
		k=(k>>1);
		
		b=matrixmul(b,b);
	}
	
	return ans;  
}

int main()
{
    int t,k;
    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]);
				if(i!=j)
					I.m[i][j]=0;
				else
					I.m[i][j]=1;
			}
        }   
	  matrix ans= qiuckpow(k);
      int  sum=0;
      for( i=0;i<n;i++)
		
			sum=(sum+ans.m[i][i])%mod;
		
		printf("%d\n",sum);     
		
		
    }
    
    
    
	return 0;    
}