#include<bits/stdc++.h>
using namespace std;
int n,m,k,P,inv[10005],f[105],g[2][105][10005][2],s,ans;
inline void Add(int &x,int y){
	if((x+=y)>=P)x-=P;
}
void solve1(){
	f[0]=1;
	for(int i=1;i<=m;++i)
		for(int j=i-1;~j;--j)Add(f[j+1],1ll*inv[i]*f[j]%P);
	printf("%d\n",f[n]);
}
void solve2(){
	g[0][0][0][0]=1;
	for(int i=1,p=1,q=0;i<=m;s+=i,++i,p^=1,q^=1)
		for(int j=0;j<i;++j)for(int k=0;k<=s;++k)for(int o=0;o<2;++o){
			int w=g[q][j][k][o];
			g[q][j][k][o]=0;
			Add(g[p][j][k][o],w);
			Add(g[p][j+1][k][o],1ll*w*inv[i]%P);
			Add(g[p][j+1][k+i][o],1ll*(P-w)*inv[i]%P);
			if(!o)Add(g[p][j+1][k+i][1],1ll*w*i%P);
		}
	for(int i=0;i<=s;++i)Add(ans,1ll*g[m&1][n][i][1]*inv[i]%P);
	printf("%d\n",ans);
}
int main(){
	scanf("%d%d%d%d",&n,&m,&k,&P);
	inv[1]=1;
	for(int i=2;i<=10000;++i)inv[i]=1ll*(P-P/i)*inv[P%i]%P;
	if(k==1)solve1();
	else solve2();
	return 0;
}