由题意可知给出的一定是一个置换群(但是没有单位置换),且答案就是在该置换群下等价类的个数
直接用引理就行了
由于颜色只有三种,所以可以用一个DP来求各个置换的不动点
#include<cstdio>
#include<cstring>
using namespace std;
int n,R,G,B,m,TT,ans,tot,a[65],cl[65],cnt[65],f[65][65][65];
inline int read(){
int ret=0,fh=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')fh=-fh;ch=getchar();}
while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar();
return ret*fh;
}
int inv(int x){
int ret=1,t=x,y=TT-2;
while(y){
if(y&1) ret=ret*t%TT;
t=t*t%TT,y>>=1;
}
return ret;
}
int calc(){
memset(cl,0,sizeof(cl));
memset(cnt,0,sizeof(cnt));
memset(f,0,sizeof(f));
tot=0;
for(int i=1;i<=n;i++){
if(cl[i]==0){
tot++;
for(int p=i;cl[p]==0;p=a[p]) cnt[tot]++,cl[p]=tot;
}
}
f[0][0][0]=1;
for(int i=1;i<=tot;i++)
for(int jr=R;jr>=0;jr--)
for(int jg=G;jg>=0;jg--)
for(int jb=B;jb>=0;jb--){
if(jr>=cnt[i]) f[jr][jg][jb]=(f[jr][jg][jb]+f[jr-cnt[i]][jg][jb])%TT;
if(jg>=cnt[i]) f[jr][jg][jb]=(f[jr][jg][jb]+f[jr][jg-cnt[i]][jb])%TT;
if(jb>=cnt[i]) f[jr][jg][jb]=(f[jr][jg][jb]+f[jr][jg][jb-cnt[i]])%TT;
}
return f[R][G][B];
}
int main(){
freopen("P1004.in","r",stdin);
freopen("P1004.out","w",stdout);
R=read(),G=read(),B=read();
n=R+G+B,m=read(),TT=read();
for(int i=1;i<=m;i++){
for(int j=1;j<=n;j++) a[j]=read();
ans=(ans+calc())%TT;
}
for(int i=1;i<=n;i++) a[i]=i;
ans=(ans+calc())%TT;
ans=ans*inv(m+1)%TT;
printf("%d\n",ans);
return 0;
}