#include<bits/stdc++.h>
// #define int long long 
#define ll long long
using namespace std;
int n;
ll MOd;
ll g[505][505];
ll f[505][505];
ll C[505][505],l[505];
ll fc[256035],inv[256035];
ll ksm(ll x,int y){
	ll sum=1;
	while(y){
		if(y&1)sum=sum*x%MOd;
		x=x*x%MOd;
		y>>=1;
	}
	return sum;
}
int AL(int x,int y){
	int nn=x*y+n-1,mm=x*y-1;
	return fc[nn]*inv[nn-mm]%MOd*inv[mm]%MOd;
}
signed main(){
	cin.tie(0)->sync_with_stdio(0);
	//freopen("1.in","r",stdin);
	//freopen("1.out","w",stdout);
	cin>>n>>MOd;
	fc[0]=1;
	for(int i=1;i<=n*n+2*n;++i)fc[i]=fc[i-1]*i%MOd;
	inv[n*n+2*n]=ksm(fc[n*n+2*n],MOd-2);
	for(int i=n*n+2*n-1;i>=0;--i)inv[i]=inv[i+1]*(i+1)%MOd;
	for(int i=0;i<=n;++i)C[i][0]=1;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=n;++j)
			C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOd;

	for(int b=1;b<=n;++b){
		for(int x=1;x<=n;++x){
			l[x]=0;
			for(int y=1;y<=b;++y)
				l[x]=(l[x]+((b-y&1)?-1:+1)*C[b][y]*AL(x,y))%MOd;
		}
		for(int a=1;a<=n;++a){
			g[a][b]=0;
			for(int x=1;x<=n;++x)
				g[a][b]=(g[a][b]+C[a][x]*((a-x&1)?-1:1)*l[x])%MOd;
		}
	}

	for(int y=0;y<=n;++y){
		for(int a=0;a<n;++a){
			l[a]=0;
			for(int b=y;b<=n;++b)
				l[a]=(l[a]+((b-y&1)?-1:1)*C[b][y]*g[n-a][n-b])%MOd;
		}
		for(int x=0;x<=n;++x){
			f[x][y]=0;
			for(int a=x;a<n;++a)
				f[x][y]=(f[x][y]+((a-x&1)?-1:1)*C[a][x]*l[a])%MOd;
		}
	}
	for(int i=0;i<n;++i){
		for(int j=0;j<n;++j)
			cout<<(f[i][j]+MOd)%MOd<<' ';
		cout<<'\n';
	}
	cerr<<clock()*1.0/CLOCKS_PER_SEC<<' ';
	return 0;
}