#include<bits/stdc++.h>
using namespace std;
int n,m,K,mod;
int qpow(int a,int b){
	int c=1;for(;b;b>>=1){
		if(b&1)c=1ll*a*c%mod;
		a=1ll*a*a%mod;
	}
	return c;
}
int dp[110][110],inv[10010];
int f[110][5110][2],f2[110][5110][2];
int main(){
	scanf("%d%d%d%d",&n,&m,&K,&mod);
	for(int i=1;i<=10000;i++)inv[i]=qpow(i,mod-2);
	if(K==1){
		dp[0][0]=1;
		for(int i=1;i<=m;i++){
			for(int j=0;j<=n;j++){
				(dp[i][j]+=dp[i-1][j])%=mod;
				if(j)(dp[i][j]+=1ll*dp[i-1][j-1]*inv[i]%mod)%=mod;
			}
		}
		printf("%d",dp[m][n]);return 0;
	}
	f[0][0][0]=1;
	for(int i=1;i<=m;i++){
		memset(f2,0,sizeof(f2));
		for(int j=0;j<n;j++)for(int k=0;k<=i*(i-1)/2;k++)for(int o=0;o<2;o++)if(f[j][k][o]){
			(f2[j][k][o]+=f[j][k][o])%=mod;
			(f2[j+1][k+i][o]-=1ll*inv[i]*f[j][k][o]%mod)%=mod;
			(f2[j+1][k][o]+=1ll*inv[i]*f[j][k][o]%mod)%=mod;
			if(!o)(f2[j][k+i][1]+=1ll*i*f[j][k][o]%mod)%=mod;
		}
		for(int j=0;j<n;j++)for(int k=0;k<=i*(i+1)/2;k++)for(int o=0;o<2;o++)
		f[j][k][o]=f2[j][k][o];
	}
	int ans=0;
	for(int i=0;i<=m*(m+1)/2;i++)if(f[n-1][i][1])
	(ans+=1ll*inv[i]*f[n-1][i][1]%mod)%=mod;
	return printf("%d",(ans+mod)%mod),0;
}