ll@(h,w,k,m),hw=h*w;
unionFind u('m',hw);
mint().setmod(m);
mint a,z,c[hw+1];
rep(i,hw){
	a+=mint(1)/(i+1);
}
rep(i,hw+1){
	c[i]=mint(i)**k;
}
rep(b,1<<hw){
	u.init();
	rep(j,hw){
		if(b&1<<j){
			if(j%w && b&1<<j-1){
				u(j-1,j);
			}
			if(j>=w && b&1<<j-w){
				u(j-w,j);
			}
		}
	}
	rep(j,hw){
		if(b&1<<j){
			if(u(j)==j){
				z+=c[u.size(j)];
			}
		}
	}
}
wt(a*z/mint(2)**hw);