BZOJ1004: [HNOI2008]Cards（Burnside引理）

由题意可知给出的一定是一个置换群（但是没有单位置换），且答案就是在该置换群下等价类的个数

直接用引理就行了

由于颜色只有三种，所以可以用一个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;
}