#include<bits/stdc++.h>
using namespace std;
const int N=510,M=N*N;
int n,mod,fac[M],inv[M],f[N][N],g[N][N];
int C(int n,int m){
	return n<m?0:1ll*fac[n]*inv[m]%mod*inv[n-m]%mod;
}
signed main(){
	scanf("%d%d",&n,&mod);
	fac[0]=inv[0]=inv[1]=1;
	for(int i=2;i<M;i++) inv[i]=1ll*inv[mod%i]*(mod-mod/i)%mod;
	for(int i=1;i<M;i++)
		fac[i]=1ll*fac[i-1]*i%mod,inv[i]=1ll*inv[i-1]*inv[i]%mod;
	for(int x=0;x<=n;x++)
		for(int y=0;y<=n;y++) f[x][y]=C(n+x*y-1,x*y-1);
	for(int y=0;y<=n;y++)
		for(int x=n;x>=0;x--)
			for(int i=0;i<x;i++)
				f[x][y]=(f[x][y]+1ll*((x-i)&1?mod-1:1)*C(x,i)%mod*f[i][y])%mod;
	for(int x=0;x<=n;x++)
		for(int y=n;y>=0;y--)
			for(int j=0;j<y;j++)
				f[x][y]=(f[x][y]+1ll*((y-j)&1?mod-1:1)*C(y,j)%mod*f[x][j])%mod;
	for(int x=0;x<=n;x++)
		for(int y=0;y<=n;y++) g[x][y]=f[n-x][n-y];
	for(int y=0;y<=n;y++)
		for(int x=0;x<=n;x++)
			for(int i=x+1;i<=n;i++)
				g[x][y]=(g[x][y]+1ll*((i-x)&1?mod-1:1)*C(i,x)%mod*g[i][y])%mod;
	for(int x=0;x<=n;x++)
		for(int y=0;y<=n;y++)
			for(int j=y+1;j<=n;j++)
				g[x][y]=(g[x][y]+1ll*((j-y)&1?mod-1:1)*C(j,y)%mod*g[x][j])%mod;
	for(int x=0;x<n;x++,puts(""))
		for(int y=0;y<n;y++)
			printf("%d ",g[x][y]);
	return 0;
}