#include<bits/stdc++.h>
using namespace std;
int n,mod,fac[250505],inv[250505],dp[505][505],f[505][505],g[505][505],c[505][505];
int qpow(int x,int y)
{
	int res=1;
	while(y)
	{
		if(y&1)res=1ll*res*x%mod;
		x=1ll*x*x%mod;y/=2;
	}
	return res;
}
void Add(int &x,int y){x=(x+y>=mod?x+y-mod:x+y);}
int C(int x,int y){return (x<y?0:1ll*fac[x]*inv[y]%mod*inv[x-y]%mod);}
int main()
{
    cin.tie(0)->sync_with_stdio(0);
    cin>>n>>mod;fac[0]=inv[0]=c[0][0]=1;
    for(int i=1;i<=n;i++)
    {
    	c[i][0]=1;
    	for(int j=1;j<=i;j++)c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
	}
    for(int i=1;i<=n*n+n;i++)fac[i]=1ll*fac[i-1]*i%mod;
    inv[n*n+n]=qpow(fac[n*n+n],mod-2);
    for(int i=n*n+n-1;i>=1;i--)inv[i]=1ll*inv[i+1]*(i+1)%mod;
    for(int i=1;i<=n;i++)
    {
    	for(int j=1;j<=n;j++)g[i][j]=C(i*j+n-1,n);
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			f[i][j]=g[i][j];
			for(int k=1;k<j;k++)
			{
				if((k-j)&1)Add(f[i][j],mod-1ll*c[j][k]*g[i][k]%mod);
				else Add(f[i][j],1ll*c[j][k]*g[i][k]%mod); 
			} 
		}
	}
	for(int j=1;j<=n;j++)
	{
		for(int i=1;i<=n;i++)
		{
			dp[i][j]=f[i][j];
			for(int k=1;k<i;k++)
			{
				if((k-i)&1)Add(dp[i][j],mod-1ll*c[i][k]*f[k][j]%mod);
				else Add(dp[i][j],1ll*c[i][k]*f[k][j]%mod);
			}
		}
	}
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)g[i][j]=dp[n-i][n-j];
	}
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			f[i][j]=g[i][j];
			for(int k=j+1;k<n;k++)
			{
				if((k-j)&1)Add(f[i][j],mod-1ll*c[k][j]*g[i][k]%mod);
				else Add(f[i][j],1ll*c[k][j]*g[i][k]%mod); 
			} 
		}
	}
	for(int j=0;j<n;j++)
	{
		for(int i=0;i<n;i++)
		{
			dp[i][j]=f[i][j];
			for(int k=i+1;k<n;k++)
			{
				if((k-i)&1)Add(dp[i][j],mod-1ll*c[k][i]*f[k][j]%mod);
				else Add(dp[i][j],1ll*c[k][i]*f[k][j]%mod);
			}
		}
	}
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)cout<<dp[i][j]<<' ';
		cout<<'\n';
	}
	return 0;
}