代码:

1 #include 
      
        2 #include 
       
         3 #include 
        
          4 #include 
         
           5 using namespace std; 6 #define maxn 44 7 #define mod 10 8 int n; 9 struct matrix{10     int f[maxn][maxn];11 };12 matrix sum(matrix a,matrix b){13     int i,j;14     matrix s;15     for(i=0;i
          
           >1;43     }44     return s;45 }46 matrix work(matrix e,int k){47     matrix s,a,b;48     if(k==1)49         return e;50     a=work(e,k/2);51     s=sum(a,mul(a,pows(e,k/2)));52     if(k&1)53         s=sum(s,pows(e,k));54     return s;55 }56 int main(){57     int k;58     while( cin>>n>>k){59         if(n==0)60             break;61         int i,j,a;62         matrix e;63         for(i=0;i
           
            >a;68 e.f[i][j]=a%10;69 }70 }71 e=work(e,k);72 for(i=0;i