【BZOJ】【P1004】【HNOI2008】【Cards】【题解】【Burnside引理+dp】

传送门：http://www.lydsy.com/JudgeOnline/problem.php?id=1004

等价类计数用Burnside定理,需要计算不动点个数

dp计算

Code:

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
int Sa,Sb,Sc,m,n;
LL p;
int N;
int w[3601],fac[101];
int dp[21][21][21],a[61];
int power(int x,int k,int p){
int ans=1;
for(;k;k>>=1){
if(k&1)ans=(ans*x)%p;
x=(x*x)%p;
}return ans;
}
int inv(int x){
return power(x,p-2,p);
}
int ans;
int main(){
cin>>Sa>>Sb>>Sc>>m>>p;n=Sa+Sb+Sc;
fac[0]=1;for(int i=1;i<=n;i++)fac[i]=fac[i-1]*i%p;
ans=fac
*inv(fac[Sa])%p*inv(fac[Sb])*inv(fac[Sc])%p;
for(int i=1;i<=m;i++){
for(int j=1;j<=n;j++)scanf("%d",&a[j]);
int vis[61];memset(vis,0,sizeof vis);
for(int j=1;j<=n;j++){
if(!vis[j]){
N++;
for(int x=j;!vis[x];x=a[x]){
vis[x]=1;w
++;
}
}
}
}
dp[0][0][0]=1;
for(int i=1;i<=N;i++)
for(int j=Sa;j>=w[i];j--)
for(int k=Sb;k>=w[i];k--)
for(int l=Sc;l>=w[i];l--)
dp[j][k][l]+=dp[j-w[i]][k][l]+dp[j][k-w[i]][l]+dp[j][k][l-w[i]],
dp[j][k][l]%=p;
cout<<(LL)(ans+dp[Sa][Sb][Sc])*inv(m+1)%p<<endl;
return 0;
}