[BZOJ1004][HNOI2008][Burnside引理][DP]Cards

[Problem Description]

小春现在很清闲,面对书桌上的N张牌,他决定给每张染色,目前小春只有3种颜色:红色,蓝色,绿色.他询问Sun有多少种染色方案,Sun很快就给出了答案.进一步,小春要求染出Sr张红色,Sb张蓝色,Sg张绝色.他又询问有多少种方案,Sun想了一下,又给出了正确答案. 最后小春发明了M种不同的洗牌法,这里他又问Sun有多少种不同的染色方案.两种染色方法相同当且仅当其中一种可以通过任意的洗牌法(即可以使用多种洗牌法,而每种方法可以使用多次)洗成另一种.Sun发现这个问题有点难度,决定交给你,答案可能很大,只要求出答案除以P的余数(P为质数).

[Algorithm]

Polya定理，DP

[Analysis]

这是一个经典的Burnside引理的应用。唯一的亮点在于这里它的颜色个数是限定的，所以不能够用Polya定理去做。用一个三维的DP去求每一个置换不动点的方案数即可。最后求平均的时候由于要用到mod的除法，所以用逆元

[Pay Attention]

置换要考虑不变的那一种

[Code]

/**************************************************************
    Problem: 1004
    User: gaotianyu1350
    Language: C++
    Result: Accepted
    Time:152 ms
    Memory:1332 kb
****************************************************************/
 
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <iostream>
using namespace std;
 
const int MAXN     = 66;
const int MAXCOLOR = 22;
const int MAXM     = 70;
 
int f[MAXCOLOR][MAXCOLOR][MAXCOLOR];
int ca, cb, cc, m, n, MOD;
bool check[MAXN];
int xunhuan[MAXN], xi[MAXM][MAXN];
 
void exgcd(int a, int b, int &d, int &x, int &y)
{
    if (!b)
        d = a, x = 1, y = 0;
    else
        exgcd(b, a % b, d, y, x), y -= x * (a / b);
}
 
inline int inv(int a)
{
    int x, y, d;
    exgcd(a, MOD, d, x, y);
    if (d != 1) return -1;
    else return (x + MOD) % MOD;
}
 
inline void up(int &a, int b)
{
    a = (a + b) % MOD;
}
 
int main()
{
    scanf("%d%d%d%d%d", &ca, &cb, &cc, &m, &MOD);
    n = ca + cb + cc;
    for (int i = 1; i <= m; i++)
    for (int j = 1; j <= n; j++)
        scanf("%d", &xi[i][j]);   
    for (int i = 1; i <= n; i++)
        xi[m + 1][i] = i;
    m++;
    int ans = 0;
    for (int i = 1; i <= m; i++)
    {
        memset(f, 0, sizeof(f));
        memset(check, 0, sizeof(check));
        memset(xunhuan, 0, sizeof(xunhuan));
        int tot = 0;
        for (int j = 1; j <= n; j++)
            if (!check[j])
            {
                xunhuan[++tot] = 1; 
                check[j] = true;
                int now = xi[i][j];
                while (now != j)
                {
                    xunhuan[tot]++;
                    check[now] = true;
                    now = xi[i][now];
                }
            }
        f[0][0][0] = 1;
        for (int i = 1; i <= tot; i++)
            for (int a = ca; a >= 0; a--)
            for (int b = cb; b >= 0; b--)
            for (int c = cc; c >= 0; c--)
            {
                if (a >= xunhuan[i])
                    up(f[a][b][c], f[a - xunhuan[i]][b][c]);
                if (b >= xunhuan[i])
                    up(f[a][b][c], f[a][b - xunhuan[i]][c]);
                if (c >= xunhuan[i])
                    up(f[a][b][c], f[a][b][c - xunhuan[i]]); 
            }
        up(ans, f[ca][cb][cc]);
    }
    ans = (ans * inv(m)) % MOD;
    printf("%d\n", ans);
}