#include <cstdio>
using namespace std;
typedef long long ll;
int read(){
	char c=getchar();int x=0;
	while(c<48||c>57) c=getchar();
	do x=(x<<1)+(x<<3)+(c^48),c=getchar();
	while(c>=48&&c<=57);
	return x;
}
int n,m,k,p;
const int N=103;
int inv[N*N];
int fs[N],sm;
int f[N][N*N][2];
int main(){
	n=read();m=read();k=read();p=read();
	inv[1]=1;
	if(k==1){
		for(int i=2;i<=m;++i) inv[i]=(ll)inv[p%i]*(p-p/i)%p;
		fs[0]=1;
		for(int i=1;i<=m;++i)
			for(int j=i;j;--j) fs[j]=(fs[j]+(ll)fs[j-1]*inv[i])%p;
		printf("%d\n",fs[n]);
		return 0;
	}
	for(int i=2;i<=m*m;++i) inv[i]=(ll)inv[p%i]*(p-p/i)%p;
	int sm=0;
	f[0][0][0]=1;
	for(int i=1;i<=m;++i){
		sm+=i;
		for(int j=n;j;--j)
			for(int x=0;x<=sm;++x)
				for(int t=0;t<2;++t){
					int &dp=f[j][x][t];
					dp=(dp+(ll)inv[i]*f[j-1][x][t])%p;
					if(x>=i){
						dp=(dp-(ll)inv[i]*f[j-1][x-i][t])%p;
						if(dp<0) dp+=p;
						if(t) dp=(dp+(ll)i*f[j-1][x-i][0])%p;
					}
				}
	}
	int res=0;
	for(int i=0;i<=sm;++i)
		res=(res+(ll)f[n][i][1]*inv[i])%p;
	printf("%d\n",res);
	return 0;
}