POJ 3233

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
using namespace std;
int n,k,m;
struct mat
{
	int a[35][35];
};
mat add(mat a,mat b)
{
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			if(b.a[i][j]) a.a[i][j]=(a.a[i][j]+b.a[i][j])%m;
		}
	}
	return a;
}
mat mul(mat a,mat b)
{
	mat ret;
	memset(ret.a,0,sizeof(ret.a));
	for(int i=0;i<n;i++)
	{
		for(int k=0;k<n;k++)
		{
			if(a.a[i][k])
			for(int j=0;j<n;j++)
			{
				if(b.a[k][j])
				ret.a[i][j]=(ret.a[i][j]+a.a[i][k]*b.a[k][j]%m)%m;
			}
		}
	}
	return ret;
}
mat mpower(mat a,int li)
{
	mat I;
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			I.a[i][j]=(i==j);
		}
	}
	while(li)
	{
		if(li&1) I=mul(I,a);
		a=mul(a,a);
		li>>=1;
	}
	return I;
}
mat dfs(int li,mat A,mat I)
{
	if(li==1) return A;
	mat s=dfs(li>>1,A,I);
	s=mul(s,add(I,mpower(A,li/2)));
	if(li&1) s=add(s,mpower(A,li));
	return s;
}
int main()
{
	scanf("%d%d%d",&n,&k,&m);
	mat A,I;
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			scanf("%d",&A.a[i][j]);
			A.a[i][j]%=m;
			I.a[i][j]=(i==j);
		}
	}
	mat tmp = dfs(k,A,I);
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n-1;j++)
		{
			printf("%d ",tmp.a[i][j]);
		}
		printf("%d\n",tmp.a[i][n-1]);
	}
	return 0;
}