R=range
H,W,K,M=map(int,input().split())
HW=H*W
p=1 << HW
found=[0]*p
count=[0]*(HW+1)
dfs=[]
for h in R(H):
	for w in R(W):
		i=h*W+w
		s=1<<i
		found[s]=1
		bs=0
		for dx,dy in zip([-1,0,1,0],[0,-1,0,1]):
			if 0<=h+dx<H and 0<=w+dy<W:
				j=i+dx*W+dy
				bs|=(1<<j)
		dfs+=[[s,bs]]
while dfs:
	s,bs=dfs.pop()
	c=sum((s>>d)&1for d in R(HW))
	count[c]+=1<<(HW-c-sum((bs>>d)&1for d in R(HW)))
	for i in R(HW):
		if(bs>>i)&1<1:continue
		h,w=i//W,i%W
		t=s|(1<<i)
		if found[t]:continue;
		found[t]=1
		bt=bs^(1<<i)
		for dx,dy in zip([-1,0,1,0],[0,-1,0,1]):
			if 0<=h+dx<H and 0<=w+dy<W:
				j=i+dx*W+dy
				if(s>>j)&1<1:bt|=(1<<j)
		dfs+=[[t,bt]]
print(sum(pow(c,K,M)*count[c]for c in R(1,HW+1))*sum(pow(HW-n+1,-1,M)for n in R(1,HW+1))*pow(2,-HW,M)%M)