#include<bits/stdc++.h>
using namespace std;
const int N=505;
int n,P,f[N][N],g[N][N],c[N][N];
int fac[N*N],inv[N*N];
inline int md(int x){
	return x>=P?x-P:x;
}
int C(int n,int m){
	return 1ll*fac[n]*inv[m]%P*inv[n-m]%P;
}
void init(){
	for(int i=0;i<N;++i)c[i][0]=1;
	for(int i=1;i<N;++i)for(int j=1;j<=i;++j)c[i][j]=md(c[i-1][j]+c[i-1][j-1]);
	fac[0]=inv[0]=fac[1]=inv[1]=1;
	for(int i=2;i<N*N;++i)fac[i]=1ll*fac[i-1]*i%P,inv[i]=1ll*(P-P/i)*inv[P%i]%P;
	for(int i=2;i<N*N;++i)inv[i]=1ll*inv[i-1]*inv[i]%P;
}
int main(){
	scanf("%d%d",&n,&P);
	init();
	for(int i=1;i<=n;++i)for(int j=1;j<=n;++j)f[i][j]=C(i*j+n-1,n);
	for(int i=1;i<=n;++i)
		for(int j=1;j<=n;++j)
			for(int k=1;k<j;++k)
				f[i][j]=(1ll*(P-f[i][k])*c[j][k]+f[i][j])%P;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=n;++j)
			for(int k=1;k<i;++k)
				f[i][j]=(1ll*(P-f[k][j])*c[i][k]+f[i][j])%P;
	for(int i=1;i<=n;++i)for(int j=1;j<=n;++j)g[n-i][n-j]=f[i][j];
	for(int i=n-1;~i;--i)
		for(int j=n-1;~j;--j)
			for(int k=j+1;k<n;++k)
				g[i][j]=(1ll*(P-g[i][k])*c[k][j]+g[i][j])%P;
	for(int i=n-1;~i;--i)
		for(int j=n-1;~j;--j)
			for(int k=i+1;k<n;++k)
				g[i][j]=(1ll*(P-g[k][j])*c[k][i]+g[i][j])%P;
	for(int i=0;i<n;++i)for(int j=0;j<n;++j)printf("%d%c",g[i][j]," \n"[j==n-1]);
	return 0;
}