#include <bits/stdc++.h>
using namespace std;
const int Maxn=105;
int f[Maxn][Maxn][Maxn],g[Maxn][Maxn][Maxn],C[Maxn][Maxn];
int main(){
    int n,m,c,mod;cin>>n>>m>>c>>mod;
    for(int i=0;i<105;C[i++][0]=1)
        for(int j=1;j<=i;j++)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
    f[0][0][0]=1;
    for(int k=1;k<=c;k++){
        static int base[Maxn*Maxn];
        for(int i=0,w=1;i<=n*m;i++,w=1ll*w*(c-k+1)%mod)base[i]=w;
        for(int j=m,b1=1,b2=1;j>=0;j--,b1=1ll*b1*k%mod,b2=1ll*b2*(k-1)%mod)
            for(int i=0;i<=n;i++)
                for(int d=0,w=1;i+d<=n;d++,w=1ll*w*(b1-b2)%mod)
                    g[i+d][j][k]=(g[i+d][j][k]+1ll*f[i][j][k-1]*C[n-i][d]%mod*base[d*j]%mod*w)%mod;
        for(int i=0,w=1;i<=n*m;i++,w=1ll*w*k%mod)base[i]=w;
        for(int i=0,b1=1,b2=1;i<=n;i++,b1=1ll*b1*(c-k+1)%mod,b2=1ll*b2*(c-k)%mod)
            for(int j=0;j<=m;j++)
                for(int d=0,w=1;j+d<=m;d++,w=1ll*w*(b1-b2)%mod)
                    f[i][j+d][k]=(f[i][j+d][k]+1ll*g[i][j][k]*C[m-j][d]%mod*base[(n-i)*d]%mod*w)%mod;
    }
    cout<<(f[n][m][c]+mod)%mod<<"\n";
    return 0;
}