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