POJ3233

#include <iostream>#include <cstdio>#include <cstring>using namespace std;struct matrix{int a[32][32];};int m,n;matrix Sum(matrix x,matrix y){matrix temp;memset(temp.a,0,sizeof(temp.a));int i,j,k;for(i=0;i<n;i++){for(j=0;j<n;j++){for(k=0;k<n;k++){temp.a[i][j]=(temp.a[i][j]+x.a[i][k]*y.a[k][j]%m)%m;}}}return temp;}matrix Add(matrix x,matrix final){int i,j;matrix temp;for(i=0;i<n;i++){for(j=0;j<n;j++){temp.a[i][j]=(final.a[i][j]+x.a[i][j])%m;} }  return temp;}matrixMul(matrix c,int num){matrix origin ;int i,j;memset(origin.a,0,sizeof(origin.a));for(i=0;i<n;i++)origin.a[i][i]=1;while(num!=0){if(num&1)origin=Sum(origin,c);c=Sum(c,c);num>>=1;}return origin;}matrix Num(matrix x,int k){matrix temp;int i;if(k==1)return x;memset(temp.a,0,sizeof(temp.a));for(i=0;i<n;i++)temp.a[i][i]=1;temp=Add(temp,Mul(x,k>>1));temp = Sum(temp,Num(x,k>>1)) ;if ( k&1 )        temp = Add(temp,Mul(x,k)) ; return temp ;}int main(){int i,j,k;matrix temp;while(~scanf("%d",&n)){scanf("%d%d",&k,&m);for(i=0;i<n;i++){for(j=0;j<n;j++){scanf("%d",&temp.a[i][j]);}} temp=Num(temp,k);for(i=0;i<n;i++){for(j=0;j<n;j++){if(j==0){printf("%d",temp.a[i][j]);}else printf(" %d",temp.a[i][j]);}printf("\n");}}return 0;}