#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 ;
	}
}
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;
	}
	return 0;
}