#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;
    auto ksm=[&](int x,int k){
        int ret=1,tmp=x;
        while(k){
            if(k&1)ret=1ll*ret*tmp%mod;
            tmp=1ll*tmp*tmp%mod;k>>=1;
        }
        return ret;
    };
    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=0;j<=m;j++){
            int coef=ksm(k,m-j)-ksm(k-1,m-j);
            for(int i=0;i<=n;i++){
                int v=f[i][j][k-1];
                if(!v)continue;
                for(int d=0,w=1;i+d<=n;d++,w=1ll*w*coef%mod){
                    int V=v;
                    V=1ll*V*C[n-i][d]%mod;
                    V=1ll*V*base[d*j]%mod;
                    V=1ll*V*w%mod;
                    g[i+d][j][k]=(g[i+d][j][k]+V)%mod;
                }
            }
        }
        for(int i=0,w=1;i<=n*m;i++,w=1ll*w*k%mod)base[i]=w;
        for(int i=0;i<=n;i++){
            int coef=ksm(c-k+1,i)-ksm(c-k,i);
            for(int j=0;j<=m;j++){
                int v=g[i][j][k];
                if(!v)continue;
                for(int d=0,w=1;j+d<=m;d++,w=1ll*w*coef%mod){
                    int V=v;
                    V=1ll*V*C[m-j][d]%mod;
                    V=1ll*V*base[(n-i)*d]%mod;
                    V=1ll*V*w%mod;
                    f[i][j+d][k]=(f[i][j+d][k]+V)%mod;
                }
            }
        }
    }
    cout<<(f[n][m][c]+mod)%mod<<"\n";
    return 0;
}