//Author: Kevin
#include<bits/stdc++.h>
//#pragma GCC optimize("O2")
using namespace std;
#define ll long long
#define ull unsigned ll
#define pb emplace_back
#define mp make_pair
#define ALL(x) (x).begin(),(x).end()
#define rALL(x) (x).rbegin(),(x).rend()
#define srt(x) sort(ALL(x))
#define rev(x) reverse(ALL(x))
#define rsrt(x) sort(rALL(x))
#define sz(x) (int)(x.size())
#define inf 0x3f3f3f3f
#define pii pair<int,int>
#define lb(v,x) (int)(lower_bound(ALL(v),x)-v.begin())
#define ub(v,x) (int)(upper_bound(ALL(v),x)-v.begin())
#define uni(v) v.resize(unique(ALL(v))-v.begin())
#define longer __int128_t
void die(string S){puts(S.c_str());exit(0);}
using LL=ll;
int p;
struct Mod
{
	LL m, p;
	void init(int pp) { m = ((__int128)1 << 64) / pp; p = pp; }
	LL operator ()(LL x)
	{
		x=x - ((__int128(x) * m) >> 64) * p;
		if(x<0) x+=p;
		return x;
	}
} mod;
int n,k;
int S1;
int S2;
int A[1000100],B[1000100],C[1000100],D[1000100];
int f[444],g[1000100];
ll tr1[402][82000];
ll tr2[402][82000];
int ipw2[1001000];
ll Binom[1002][1002];
ll fact[1002],rfact[1002];
ll ksm(ll a,ll b)
{
	ll ans=1;
	while(b)
	{
		if(b&1) ans=mod(ans*a);
		b>>=1;
		a=mod(a*a);
	}
	return ans;
}
ll Coef[1002][1002];
ll Chk[1002][1002];
int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>k>>p;
	mod.init(p);
	for(int i=0;i<1002;i++)
		Binom[i][i]=Binom[i][0]=1;
	for(int i=2;i<1002;i++)
		for(int j=1;j<i;j++)
			Binom[i][j]=mod(Binom[i-1][j-1]+Binom[i-1][j]);
	ipw2[0]=1;
	for(int i=1;i<1001000;i++)
		ipw2[i]=mod(1ll*ipw2[i-1]*(p+1)/2);
	fact[0]=1;
	for(int i=1;i<1002;i++)
		fact[i]=mod(fact[i-1]*i);
	for(int i=0;i<1002;i++)
		for(int j=0;j+j<=i;j++)
			Coef[i][j]=mod(mod(mod(Binom[i][j]*Binom[i-j][j])*fact[j])*ipw2[j]);
	for(int i=0;i<1002;i++)
		for(int j=0;j<=i;j++)
			Chk[i][j]=mod(Binom[i][j]*fact[j]);
	S1=n+1;
	for(int i=0;i<=n+n;i++)
		for(int j=0;j<=n+n-i;j++)
			for(int k=0;k<=n+n-i-j;k++)
			{
				int l=n+n-i-j-k;
				if(-3*i-j+k+3*l==0)
				{
					A[S2]=i;
					B[S2]=j;
					C[S2]=k;
					D[S2]=l;
					S2++;
				}
			}
	for(int i=0;i<=n;i++)
		for(int j=0;j<S2;j++)
		{
			int x=i,y=n+n-i-i,z=i;
			int a=::A[j],b=::B[j],c=::C[j],d=::D[j];
			int A=a;
			int B=x-a;
			int C=b-B;
			int F=d;
			int E=z-d;
			int D=c-E;
			if(A<0||B<0||C<0||D<0||E<0||F<0) continue;
			ll SS1=0,SS2=0;
			if(B+D+F!=A+C+E) continue;
			if(A+C-B<0) continue;
			for(int v1=0;v1<=min(B,C);v1++)
			{
				ll coef=Coef[b][v1];
				int bp=b-v1-v1;
				coef=mod(coef*Binom[bp][B-v1]);
				coef=mod(coef*Chk[A][B-v1]);
				SS1+=coef;
			}
			for(int v2=0;v2<=min(D,E);v2++)
			{
				ll coef=Coef[c][v2];
				int cp=c-v2-v2;
				coef=mod(coef*Binom[cp][D-v2]);
				coef=mod(coef*Chk[A+C-B][D-v2]);
				SS2+=coef;
			}
			SS1=mod(SS1);
			SS1=mod(SS1*fact[F]);
			SS2=mod(SS2);
			tr1[i][j]=mod(SS1*SS2);
		}
	for(int i=0;i<=n;i++)
		for(int j=0;j<S2;j++)
		{
			tr2[i][j]=1;
			int X=i,Y=n+n-i-i,Z=i;
			int A=::A[j],B=::B[j],C=::C[j],D=::D[j];
			if(A+A>X||D+D>Z||B<A||C<D)
			{
				tr2[i][j]=0;
				continue;
			}
			tr2[i][j]=mod(1ll*Binom[X][A]*Binom[Z][D]);
			tr2[i][j]=mod(mod(1ll*Binom[X-A][A]*tr2[i][j])*Binom[Z-D][D]);
			tr2[i][j]=mod(1ll*tr2[i][j]*ipw2[A+D]);
			tr2[i][j]=mod(mod(1ll*tr2[i][j]*fact[A])*fact[D]);
			X-=A*2;
			Z-=D*2;
			B-=A;
			C-=D;
			int V1=X;
			int V2=Z;
			if(V1>B||V2>C)
			{
				tr2[i][j]=0;
				continue;
			}
			int V3=B-V1;
			int V4=C-V2;
			if(V3+V4!=Y)
			{
				tr2[i][j]=0;
				continue;
			}
			ll ways=0;
			for(int i=0;i<=min(V3,V4);i++) if(V3-i<=V1&&V4-i<=V2&&(V1-V3+i)==(V2-V4+i))
			{
				ll coef=mod(Coef[Y][i]);
				coef=mod(coef*fact[V1-V3+i]);
				coef=mod(coef*Binom[Y-i-i][V3-i]);
				coef=mod(coef*Chk[V1][V3-i]);
				coef=mod(coef*Chk[V2][V4-i]);
				ways+=coef;
			}
			ways=mod(ways);
			tr2[i][j]=mod(tr2[i][j]*ways);
		}
	f[0]=1;
	for(int i=1;i<=k;i++)
		if(i&1)
		{
			memset(g,0,sizeof(g));
			for(int a=0;a<=n;a++) if(f[a])
				for(int b=0;b<S2;b++)
					g[b]=mod(g[b]+1ll*f[a]*tr2[a][b]);
		}
		else
		{
			memset(f,0,sizeof(f));
			for(int b=0;b<S2;b++) if(g[b])
				for(int a=0;a<=n;a++)
					f[a]=mod(f[a]+1ll*g[b]*tr1[a][b]);
		}
	if(k&1)
	{
		LL ans=accumulate(g,g+S2,0ll)%p;
		cout<<ans<<'\n';
	}
	else
	{
		LL ans=accumulate(f,f+n+1,0ll)%p;
		cout<<ans<<'\n';
	}
	return 0;
}