C :

#include<stdio.h>
int main()
{
    int t,n,k,a[11][11],b[11][11],c[11][11],i,j,m,p;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&k);
        for(i=0;i<n;++i)
            for(j=0;j<n;++j)
            {
                scanf("%d",&a[i][j]);
                b[i][j]=0;
                if(i==j)
                    b[i][j]=1;
            }
        for(m=1;m<=k;++m)
            {
                for(i=0;i<n;++i)
                    for(j=0;j<n;++j)
                    {
                        c[i][j]=0;
                        for(p=0;p<n;++p)
                        c[i][j]=c[i][j]+b[i][p]*a[p][j];
                    }
                for(i=0;i<n;++i)
                    for(j=0;j<n;++j)
                        b[i][j]=c[i][j]%9973;
           }
            for(i=j=0;i<n;++i)
                j=j+b[i][i];
            printf("%d\n",j%9973);
    }
    return 0;
}

C++ :

#include <stdio.h>
 #include <stdlib.h>
 #include <string.h>
 typedef struct Node
 {
     int m[11][11];
 }Matrix;
 Matrix init, unit;  //初始化输入矩阵，单位矩阵如果用递归写Pow函数可以不用单位矩阵
 int n, K;
 void Init( )  //初始化
 {
     scanf( "%d%d", &n, &K );
     for( int i=0; i<n; ++ i )
         for( int j=0; j<n; ++ j )
         {
             scanf( "%d", &init.m[i][j] );
             unit.m[i][j]=( i == j );
         }
 }
 
 Matrix Mul( 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];
             }
             c.m[i][j] %= 9973;
         } 
     }
     return c;
 }
 
 Matrix Pow( Matrix a, Matrix b, int k )
 {
     while( k>1 )
     {
         if( k&1 )  // k为奇数时
         {
             k --;
             b=Mul( a, b );
         }
         else   // k为偶数
         {
             k >>= 1;
             a=Mul( a, a );
         }
     }
     a=Mul( a, b );
     return a;
 }
 
 int main( )
 {
     int T;
     scanf( "%d", &T );
     while( T -- )
     {
         Matrix x;
         Init( );
         x=Pow( init, unit, K );
         int sum=0, i=0;
         n--;
         while( n >= 0 )
         {
             sum += x.m[n][n];
             sum%=9973;
             n --;    
         }
         printf( "%d\n", sum%9973 );
     }
 }