#include<cstdio>
using namespace std;
#define N 105900
#define M 1059
#define K 42
#define mod 164511353
int n,k,m,tp,ri=41,as,rp;
int v1[M][K],sv[M][K],c[M][M],t2[K],rv[M],s[M][M],f[M];
int pw(int a,int p){int as=1;while(p){if(p&1)as=1ll*as*a%mod;a=1ll*a*a%mod;p>>=1;}return as;}
int vf[N][K],vg[N][K],pr[K];
void solve(int p,int q)
{
	for(int i=0;i<ri;i++)for(int j=0;j<=m;j++)v1[j][i]=0;
	for(int i=0;i<=ri;i++)vf[0][i]=vg[n+m+1][i]=i==1;
	for(int i=1;i<=n+m;i++)
	{
		int vl=mod-pw(p,n-i+mod-1);
		for(int j=0;j<ri;j++)vf[i][j]=0;
		for(int j=0;j<ri;j++)vf[i][j*p%ri]=(vf[i][j*p%ri]+vf[i-1][j])%mod;
		for(int j=0;j<ri;j++)vf[i][j]=(vf[i][j]+1ll*vl*vf[i-1][j])%mod;
		if(i%n==0)
		{
			for(int j=0;j<ri;j++)vf[i][j]=0;
			vf[i][p%ri]=1;vf[i][1]+=vl;
		}
	}
	for(int i=n+m;i>=1;i--)
	{
		int vl=mod-pw(p,n-i+mod-1);
		for(int j=0;j<ri;j++)vg[i][j]=0;
		for(int j=0;j<ri;j++)vg[i][j*p%ri]=(vg[i][j*p%ri]+vg[i+1][j])%mod;
		for(int j=0;j<ri;j++)vg[i][j]=(vg[i][j]+1ll*vl*vg[i+1][j])%mod;
		if(i%n==0)for(int j=0;j<=ri;j++)vg[i][j]=j==1;
	}
	int s1=1;for(int i=1;i<=n*(q-1);i++)s1=s1*p%ri;
	int s2=pw(p,n*q);
	for(int i=0;i<=m;i++)
	{
		int l=i+1,r=i+n,s3=pw(s2,i);
		for(int p=0;p<ri;p++)for(int q=0;q<=ri;q++)
		{
			int nt=p*q*s1%ri;
			v1[i][nt]=(v1[i][nt]+1ll*vf[r][p]*vg[l][q]%mod*s3)%mod;
		}
	}
}
int r1[K][K],r2[K][K];
int is[N*10],si[K];
void dfs(int nw,int si,int pr)
{
	if(nw==k){is[si]=pr==(1<<k)-1;return;}
	dfs(nw+1,si,pr);
	si|=1<<nw;
	int p1=pr;
	for(int i=0;i<k;i++)p1|=1<<r2[i][nw];
	while(p1>pr)
	{
		int u=0;
		for(int i=0;i<k;i++)if(((p1^pr)>>i)&1)u=i;
		pr+=1<<u;
		for(int j=0;j<=k;j++)if((p1>>j)&1)p1|=1<<r1[u][j];
	}
	dfs(nw+1,si,pr);
}
void solve_1()
{
	for(int i=0;i<=m;i++)sv[i][1]=1;
	for(int i=2;i<=k;i++)if(k%i==0)
	{
		int su=0;while(k%i==0)su++,k/=i;
		solve(i,su);
		for(int p=0;p<=m;p++)
		{
			for(int j=0;j<ri;j++)t2[j]=0;
			for(int a=0;a<ri;a++)for(int b=0;b<ri;b++)
			{
				int nt=1ll*a*b%ri;
				t2[nt]=(t2[nt]+1ll*v1[p][a]*sv[p][b])%mod;
			}
			for(int j=0;j<ri;j++)sv[p][j]=t2[j];
		}
	}
}
void solve_2()
{
	for(int i=0;i<k;i++)for(int j=0;j<k;j++)scanf("%d",&r1[i][j]);
	for(int i=0;i<k;i++)for(int j=0;j<k;j++)scanf("%d",&r2[i][j]);
	dfs(0,0,0);
	for(int i=1;i<1<<k;i<<=1)
	for(int j=0;j<1<<k;j+=i*2)
	for(int k=0;k<i;k++)
	is[j+k]-=is[j+k+i];
	for(int i=0;i<1<<k;i++)
	{
		int ci=0;
		for(int j=0;j<k;j++)if((i>>j)&1)ci++;
		si[ci]+=is[i];
	}
	for(int i=0;i<ri;i++)vf[0][i]=vg[n+m+1][i]=i==1;
	for(int i=1;i<=n+m;i++)
	{
		for(int j=0;j<ri;j++)vf[i][j]=0;
		for(int l=1;l<=k;l++)if(si[l])
		{
			int vl=1ll*pw(l,i-n+mod-1)*(si[l]+mod)%mod;
			for(int j=0;j<ri;j++)vf[i][j*l%ri]=(vf[i][j*l%ri]+1ll*vl*vf[i-1][j])%mod;
		}
		if(i%n==0)
		{
			for(int j=0;j<ri;j++)vf[i][j]=0;
			for(int l=1;l<=k;l++)if(si[l])vf[i][l]=1ll*pw(l,i-n+mod-1)*(si[l]+mod)%mod;
		}
	}
	for(int i=n+m;i>=1;i--)
	{
		for(int j=0;j<ri;j++)vg[i][j]=0;
		for(int l=1;l<=k;l++)if(si[l])
		{
			int vl=1ll*pw(l,i-n+mod-1)*(si[l]+mod)%mod;
			for(int j=0;j<ri;j++)vg[i][j*l%ri]=(vg[i][j*l%ri]+1ll*vl*vg[i+1][j])%mod;
		}
		if(i%n==0)for(int j=0;j<ri;j++)vg[i][j]=j==1;
	}
	int s1=pw(k,1ll*n*(n-1)/2%(mod-1));
	for(int i=0;i<=m;i++)
	{
		int l=i+1,r=i+n;
		for(int p=0;p<=ri;p++)for(int q=0;q<=ri;q++)
		{
			int nt=p*q%ri;
			sv[i][nt]=(sv[i][nt]+1ll*vf[r][p]*vg[l][q]%mod*s1)%mod;
		}
	}
}
int main()
{
	scanf("%d%d%d%d",&n,&k,&m,&tp);
	for(int i=0;i<=m;i++)c[i][0]=c[i][i]=1;
	for(int i=2;i<=m;i++)for(int j=1;j<i;j++)c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
	s[0][0]=1;
	for(int i=1;i<=m;i++)for(int j=1;j<=i;j++)s[i][j]=1ll*j*(s[i-1][j-1]+s[i-1][j])%mod;
	if(tp==1)solve_1();
	else solve_2();
	for(int j=0;j<=m;j++)for(int i=0;i<ri;i++)rv[j]=(rv[j]+1ll*sv[j][i]*pw(2,i))%mod;
	f[0]=1;
	for(int i=0;i<=m;i++)
	{
		for(int j=0;j<=m;j++)as=(as+1ll*s[m][i]*f[j]%mod*rv[j])%mod;
		for(int j=m;j>=0;j--)f[j+1]=(f[j+1]+f[j])%mod,f[j]=1ll*f[j]*(mod-i)%mod;
		int tp=pw(2*(i+1),mod-2);for(int j=m;j>=0;j--)f[j]=1ll*f[j]*tp%mod;
	}
	printf("%d\n",as);
}