bzoj 1004 组合

代码：

//根据Burnside定理:有m个置换k钟颜色，所有本质不同的染色方案数就是每种置换的不变元素的个数的平均数。所谓不变元素就是一种染色方案
//经过置换变换后和之前一样。所以现在就是要求不变元素的个数，要想变换后和之前一样那么改置换的循环节中一定是同一种颜色，所以现在就
//是要求每个置换的循环节然后求出每个循环节染同一种颜色的方案数。因为只有3种颜色可以用三维的01背包求方案数。
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
int sr,sb,sg,n,m,p,a[70][70];
ll f[70][70][70];
ll solve(int x)
{
    bool vis[70];
    int b[70],sum=0;
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n;i++){
        if(!vis[a[x][i]]){
            b[++sum]=1;
            vis[a[x][i]]=1;
            int y=a[x][i];
            while(!vis[a[x][y]]){
                b[sum]++;
                vis[a[x][y]]=1;
                y=a[x][y];
            }
        }
    }
    memset(f,0,sizeof(f));
    f[0][0][0]=1;
    for(int h=1;h<=sum;h++)
    for(int i=sr;i>=0;i--){
        for(int j=sb;j>=0;j--){
            for(int k=sg;k>=0;k--){
                if(i>=b[h]) f[i][j][k]=(f[i][j][k]+f[i-b[h]][j][k])%p;
                if(j>=b[h]) f[i][j][k]=(f[i][j][k]+f[i][j-b[h]][k])%p;
                if(k>=b[h]) f[i][j][k]=(f[i][j][k]+f[i][j][k-b[h]])%p;
            }
        }
    }
    return f[sr][sb][sg];
}
ll pow_mod(int a,int b)
{
    if(b==0) return 1;
    ll x=pow_mod(a,b/2);
    x=x*x%p;
    if(b&1) x=x*a%p;
    return x;
}
int main()
{
    scanf("%d%d%d%d%d",&sr,&sb,&sg,&m,&p);
    n=sr+sb+sg;
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&a[i][j]);
    m++;
    for(int j=1;j<=n;j++) a[m][j]=j;
    ll ans=0;
    for(int i=1;i<=m;i++)
        ans=(ans+solve(i))%p;
    ans=(ans*pow_mod(m,p-2))%p;
    printf("%lld\n",ans);
    return 0;
}