BZOJ1004 Cards （burnside dp 扩展欧几里德）

题目大意

有三种颜色的卡片分别sa，sb，sc张，并给出m种相等的置换方式，求所有本质不同的排列方式。

---

题解：

根据burnside引理，本质不同的排列数（ans）=sigma（Zi）/m。其中Zi表示对于第i种置换，无论怎么置换，置换几遍都相同的排列的个数，既不动点的个数。根据Polya定理可知，对于某个置换的某个循环节，循环节内所有的元素都相同的排列可称为一个不动点。

由于有sa，sb，sc的约束，所以要对每一个置换进行01背包dp。每一个循环节是一个物品，循环节的大小就是物品体积。

最后求平均数的时候要用exgcd求出平均数。

---

代码：

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;

int p;
int add(int a,int b) {
    int c=a+b;
    return c>=p?c-p:c;
}
int mul(int a,int b) {
    int c=a*b;
    return c%p;
}
void exgcd(int a,int b,int &x,int &y) {
    if(b==0) {x=1,y=0; return;}
    exgcd(b,a%b,x,y);
    int tmp=x;
    x=y, y=tmp-a/b*y;
}

const int maxn=100;
int sa,sb,sc;
int n,m;
int a[maxn][maxn];
bool vis[maxn];
int sta[maxn],top=0;
int f[25][25][25];

int dp(int *to) {
    top=0;
    for(int i=1;i<=n;i++) if(!vis[i]) {
        int siz=1, cur=to[i];
        while(cur!=i) {vis[cur]=true, siz++, cur=to[cur];}
        sta[++top]=siz;
    }
    memset(f,0,sizeof f);
    f[0][0][0]=1;
    for(int i=1;i<=top;i++)
    for(int a=sa;a>=0;--a)
    for(int b=sb;b>=0;--b)
    for(int c=sc;c>=0;--c) {
        int s=sta[i];
        if(a>=s) f[a][b][c]=add(f[a][b][c],f[a-s][b][c]);
        if(b>=s) f[a][b][c]=add(f[a][b][c],f[a][b-s][c]);
        if(c>=s) f[a][b][c]=add(f[a][b][c],f[a][b][c-s]);
    }
    return f[sa][sb][sc];
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
#endif // ONLINE_JUDGE
    scanf("%d%d%d%d%d",&sa,&sb,&sc,&m,&p); m++;
    n=sa+sb+sc;

    for(int i=1;i<=n;i++) a[1][i]=i;
    for(int i=2;i<=m;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&a[i][j]);

    int ans=0;
    for(int i=1;i<=m;i++) ans=add(ans,dp(a[i]));

    int x,y;
    exgcd(m,p,x,y);
    while(x<0) x+=p;
    printf("%d\n",mul(ans,x));

    return 0;
}