#include<bits/stdc++.h>
using namespace std;
int n,m,C[110][110],K,mod,p[110][110],f[2][110][110];
void clean(int o){for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)f[o][i][j]=0;}
int main()
{
	scanf("%d%d%d%d",&n,&m,&K,&mod);
	for(int i=0;i<=max(n,m);i++)
	{
		C[i][0]=1;
		for(int j=1;j<=i;j++)C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
	}
	for(int i=1;i<=K;i++)
	{
		p[i][0]=1;
		for(int j=1;j<=max(n,m);j++)p[i][j]=1ll*p[i][j-1]*i%mod;
	}
	f[0][0][0]=1;int o=0;
	for(int k=1;k<=K;k++)
	{
		clean(o^1);
		for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)
		{
			int x=f[o][i][j],y=1ll*p[K-k+1][j]*p[k][m-j]%mod;
			for(int a=0;i+a<=n;a++,x=1ll*x*y%mod)f[o^1][i+a][j]=(f[o^1][i+a][j]+1ll*C[n-i][a]*x)%mod;
		}
		o^=1; 
		clean(o^1);
		for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)
		{
			int x=f[o][i][j],y=1ll*p[K-k+1][i]*p[k][n-i]%mod;
			for(int b=0;j+b<=m;b++,x=1ll*x*y%mod)f[o^1][i][j+b]=(f[o^1][i][j+b]+1ll*C[m-j][b]*x)%mod;
		}
		o^=1; 
		clean(o^1);
		for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)
		{
			int x=f[o][i][j],y=1ll*-p[K-k][j]*p[k][m-j]%mod;
			for(int a=0;i+a<=n;a++,x=1ll*x*y%mod)f[o^1][i+a][j]=(f[o^1][i+a][j]+1ll*C[n-i][a]*x)%mod;
		}
		o^=1; 
		clean(o^1);
		for(int i=0;i<=n;i++)for(int j=0;j<=m;j++)
		{
			int x=f[o][i][j],y=1ll*-p[K-k][i]*p[k][n-i]%mod;
			for(int b=0;j+b<=m;b++,x=1ll*x*y%mod)f[o^1][i][j+b]=(f[o^1][i][j+b]+1ll*C[m-j][b]*x)%mod;
		}
		o^=1; 
	} 
	printf("%d",(f[o][n][m]+mod)%mod);
	return 0;
}