BZOJ 1004 HNOI2008 Cards Burnside引理

标题效果：特定n张卡m换人，编号寻求等价类

数据保证这m换人加上置换群置换后本身构成

BZOJ坑爹0.0 条件不那么重要出来尼玛怎么做

Burnside引理……昨晚为了做这题硬啃了一晚上白书0.0 都快啃吐了0.0

Burnside引理：一个置换群下的等价类个数等于全部置换的不动点个数的平均值

没有接触过群论的建议去啃白书…… 网上的东西看不懂的

最后那个除法要用乘法逆元我懒得写EXGCD写了费马小定理0.0

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 70
using namespace std;
int r,g,b,m,n,p,ans;
int a[M],stack[M],top;
int f[21][21][21];
void DFS(int x)
{
	stack[top]++;
	int temp=a[x];
	a[x]=0;
	if(a[temp])
		DFS(temp);
}
int DP()
{
	int i,j,k;
	memset(f,0,sizeof f);f[0][0][0]=1;
	while(top)
	{
		for(i=r;~i;i--)
			for(j=g;~j;j--)
				for(k=b;~k;k--)
				{
					if(i>=stack[top]) f[i][j][k]+=f[i-stack[top]][j][k];
					if(j>=stack[top]) f[i][j][k]+=f[i][j-stack[top]][k];
					if(k>=stack[top]) f[i][j][k]+=f[i][j][k-stack[top]];
					f[i][j][k]%=p;
				}
		stack[top--]=0;
	}
	return f[r][g][b];
}
int KSM(int x,int y)
{
	int re=1;
	while(y)
	{
		if(y&1)re*=x,re%=p;
		x*=x,x%=p;
		y>>=1;
	}
	return re;
}
int main()
{
	int i,j;
	cin>>r>>g>>b>>m>>p;
	n=r+g+b;
	for(i=1;i<=m;i++)
	{
		for(j=1;j<=n;j++)
			scanf("%d",&a[j]);
		for(j=1;j<=n;j++)
			if(a[j])
				++top,DFS(j);
		ans+=DP(),ans%=p;
	}
	for(j=1;j<=n;j++)
		a[j]=j;
	for(j=1;j<=n;j++)
		if(a[j])
			++top,DFS(j);
	ans+=DP(),ans%=p;
	ans*=KSM(m+1,p-2),ans%=p;
	cout<<ans<<endl;
}