#include<bits/stdc++.h>
using namespace std;
inline int qread()
{
	int a=0,f=1;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch)){(a*=10)+=(ch^48);ch=getchar();}
	return a*f;
}
int n,m,k,mod,inv[10010];
inline void fmod(int &x)
{
	if(x>=mod)x-=mod;
	return ;
}
namespace SOLVE1
{
	int f[110];
	void main()
	{
		f[0]=1;
		for(int i=1;i<=m;++i)for(int j=n;j;--j)f[j]=(1ll*f[j]+1ll*f[j-1]*inv[i])%mod;
		printf("%d\n",f[n]);
		return ;
	}
}
namespace SOLVE2
{
	int ans,sum,f[2][110][10010][2];
	void main()
	{
		f[1][0][0][0]=1;
		for(int i=1,np,tp;i<=m;sum+=i,++i)
		{
			np=i&1,tp=np^1;
			for(int j=i-1;~j;--j)
			{
				for(int k=sum;~k;--k)
				{
					for(int op=0;op<=1;++op)
					{
						fmod(f[tp][j][k][op]+=f[np][j][k][op]);
						f[tp][j+1][k][op]=(1ll*f[tp][j+1][k][op]+1ll*f[np][j][k][op]*inv[i])%mod;
						f[tp][j+1][k+i][op]=(1ll*f[tp][j+1][k+i][op]+1ll*f[np][j][k][op]*(mod-inv[i]))%mod;
						if(!op)f[tp][j+1][k+i][1]=(1ll*f[tp][j+1][k+i][1]+1ll*f[np][j][k][op]*i)%mod;
						f[np][j][k][op]=0;
					}
				}
			}
		}
		for(int i=1;i<=sum;++i)ans=(1ll*ans+1ll*f[(m&1)^1][n][i][1]*inv[i])%mod;
		printf("%d\n",ans);
		return ;
	}
}
int main()
{
	n=qread(),m=qread(),k=qread(),mod=qread();
	inv[0]=inv[1]=1;
	for(int i=2;i<=n*m;++i)inv[i]=1ll*inv[mod%i]*(mod-mod/i)%mod;
	if(k==1)
	{
		SOLVE1::main();
		return 0;
	}
	SOLVE2::main();
	return 0;
}