#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 105;
int n,m,k,mod,lst,now;
ll f[2][N],inv[N*N];
void solve1()
{
    lst = 1;
    f[now][0] = 1;
    for(int i = 1;i <= m;i++)
    {
        swap(now,lst);
        memset(f[now],0,sizeof(f[now]));
        for(int j = 0;j <= n;j++)
        {
            if(!f[lst][j]) continue;
            f[now][j] = (f[now][j] + f[lst][j])%mod;
            if(j < n) f[now][j+1] = (f[now][j+1] + f[lst][j]*inv[i])%mod;
        }
    }
    cout<<f[now][n]<<"\n";
}
void solve2()
{
    cout<<"114514\n";
}
int main()
{
    cin>>n>>m>>k>>mod;
    inv[1] = 1;
    for(int i = 2;i <= n*m;i++) inv[i] = mod - mod/i*inv[mod%i]%mod;
    if(k == 1) solve1();
    if(k == 2) solve2();
    return 0;
}