#include<bits/stdc++.h>
using namespace std;
int n,m,k,mod,inv[10005];
void Add(int &x,int y){x=(x+y>=mod?x+y-mod:x+y);}
namespace Sub1
{
	int dp[105][105];
	void main()
	{
		dp[0][0]=1;
    	for(int i=1;i<=m;i++)
    	{
    		for(int j=0;j<=n;j++)
    		{
    			Add(dp[i][j],dp[i-1][j]);
    			Add(dp[i][j+1],1ll*dp[i-1][j]*inv[i]%mod);
			}
		}
		cout<<dp[m][n];
	}
}
namespace Sub2
{
	int dp[2][105][10005][2],ans;
	void main()
	{
		dp[0][0][0][0]=1;
		for(int i=1;i<=m;i++)
		{
			int p=i&1;
			for(int j=0;j<=min(i,n);j++)for(int k=0;k<=i*(i+1)/2;k++)dp[p][j][k][0]=dp[p][j][k][1]=0;
			for(int j=0;j<=min(i,n);j++)
			{
				for(int k=0;k<=i*(i+1)/2;k++)
				{
					if(dp[p^1][j][k][0])
					{
						Add(dp[p][j][k][0],dp[p^1][j][k][0]);
						Add(dp[p][j+1][k][0],1ll*dp[p^1][j][k][0]*inv[i]%mod);
						Add(dp[p][j+1][k+i][0],1ll*(mod-dp[p^1][j][k][0])*inv[i]%mod);
						Add(dp[p][j][k+i][1],1ll*i*dp[p^1][j][k][0]%mod);
					}
					if(dp[p^1][j][k][1]) 
					{
						Add(dp[p][j][k][1],dp[p^1][j][k][1]);
						Add(dp[p][j+1][k][1],1ll*dp[p^1][j][k][1]*inv[i]%mod);
						Add(dp[p][j+1][k+i][1],1ll*(mod-dp[p^1][j][k][1])*inv[i]%mod);
					}
				}
			}
		}
		for(int i=0;i<=m*(m+1)/2;i++)Add(ans,1ll*dp[m&1][n-1][i][1]*inv[i]%mod);
		cout<<ans;
	}
}
int main()
{
    cin.tie(0)->sync_with_stdio(0);
    cin>>n>>m>>k>>mod;
    inv[1]=1;
    for(int i=2;i<=10000;i++)inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
    if(k==1)Sub1::main();
    else Sub2::main();
	return 0;
}