1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 #include<cmath>
 6 using namespace std;
 7 int sr,sb,sg,n,m,p;
 8 int huan[110][110];
 9 bool visited[110];
10 int ge[110];
11 int f[110][110][110];
12 int ans=0;
13 void exgcd(int a,int b,int &x,int &y){
14     if(!b){  x=1,y=0;  return ;}
15     exgcd(b,a%b,x,y);
16     int t=x;  x=y;  y=t-a/b*y;
17 }
18 int dp(int x){
19     memset(visited,0,sizeof(visited));
20     int cnt=0,temp;
21     for(int i=1;i<=n;i++)if(!visited[i]){
22         visited[i]=true;
23         ge[++cnt]=1;  temp=i;
24         while(!visited[huan[x][temp]]){
25             ge[cnt]++;
26             visited[huan[x][temp]]=true;//这里容易蒙……
27             temp=huan[x][temp];//temp指向了这个连的下一个，再次循环时判断的是下一个的下一个……
28         }
29     }
30     memset(f,0,sizeof(f));  f[0][0][0]=1;
31     for(int t=1;t<=cnt;t++)
32         for(int i=sr;i>=0;i--)
33             for(int j=sb;j>=0;j--)
34                 for(int k=sg;k>=0;k--){
35                     if(i>=ge[t])  f[i][j][k]=(f[i][j][k]+f[i-ge[t]][j][k])%p;
36                     if(j>=ge[t])  f[i][j][k]=(f[i][j][k]+f[i][j-ge[t]][k])%p;
37                     if(k>=ge[t])  f[i][j][k]=(f[i][j][k]+f[i][j][k-ge[t]])%p;
38                 }
39     return f[sr][sb][sg];
40 }
41 int main(){//freopen("ddd.in","r",stdin);
42     cin>>sr>>sb>>sg>>m>>p;
43     n=sr+sb+sg;
44     for(int i=1;i<=m;i++)
45         for(int j=1;j<=n;j++)
46             scanf("%d",&huan[i][j]);
47     m++;
48     for(int i=1;i<=n;i++)  huan[m][i]=i;//注意还有不置换的情况
49     for(int i=1;i<=m;i++)
50         ans=(ans+dp(i))%p;
51     int x,y;
52     exgcd(m,p,x,y);
53     while(x<=0)x+=p;
54     cout<<ans*x%p<<endl;
55     return 0;
56 }