#include <bits/stdc++.h>

 using namespace std;

 const int maxn=+;

 int n,m,sr,sb,sg,p,ans;

 int a[maxn],s[maxn];

 int f[][][];

 bool vis[maxn];

 int qpow(int a,int b){

     int ret=;

     for(;b;a=(a*a)%p,b>>=) if(b&) ret=(ret*a)%p;

     return ret;

 }

 int get(){

     int cnt=; memset(f,,sizeof f); memset(s,,sizeof s); memset(vis,false, sizeof vis);

     for(int i=;i<=n;i++)if(!vis[i]){ cnt++; for(int j=i;!vis[j];j=a[j]) vis[j]=true, s[cnt]++; }

     f[][][]=;

     for(int i=;i<=cnt;i++)for(int r=sr;r>=;r--)for(int b=sb;b>=;b--)for(int g=sg;g>=;g--){

         if(r>=s[i]) f[r][b][g]=(f[r][b][g]+f[r-s[i]][b][g])%p;

         if(b>=s[i]) f[r][b][g]=(f[r][b][g]+f[r][b-s[i]][g])%p;

         if(g>=s[i]) f[r][b][g]=(f[r][b][g]+f[r][b][g-s[i]])%p;

     }

     return f[sr][sb][sg];

 }

 int main(){

     scanf("%d%d%d%d%d",&sr,&sb,&sg,&m,&p);

     n=sr+sb+sg;

     for(int i=;i<=m;i++){ for(int j=;j<=n;j++) scanf("%d",&a[j]); ans=(ans+get())%p; }

     for(int i=;i<=n;i++) a[i]=i; ans=(ans+get())%p;

     ans=(ans*qpow(m+,p-))%p;

     printf("%d\n",ans);

     return ;

 }