bzoj1004 [HNOI2008]Cards（burnside引理+背包dp）

题目保证给出的置换再加上一个单位元之后是一个群。因此我们就是要求集合M关于置换群的轨道数，根据burnside引理，就是对每一个置换求不动点个数，我们可以dfs出每一个置换的所有轮换，然后跑一个三维背包dp来计算可行的方案数。复杂度$O(m^2n^3)$。最后求一个逆元。
网上那些什么没有不动点的都是假的，HN数据水罢了-，-
比如你看这组样例：2 2 0 1 7 3 4 1 2

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define N 70
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n1,n2,n3,m,mod,n,a[N],tot=0,w[N],cnt=0,f[N][N][N],prod[N],ans=0;
bool vis[N];
inline void dfs(int x){
    vis[x]=1;cnt++;if(vis[a[x]]) return;dfs(a[x]);
}
inline int dp(){
    memset(f,0,sizeof(f));f[0][0][0]=1;
    for(int x=1;x<=tot;++x)
        for(int i=n1;i>=0;--i)
            for(int j=n2;j>=0;--j)
                for(int k=n3;k>=0;--k){
                    if(i>=w[x]) f[i][j][k]+=f[i-w[x]][j][k];
                    if(j>=w[x]) f[i][j][k]+=f[i][j-w[x]][k];
                    if(k>=w[x]) f[i][j][k]+=f[i][j][k-w[x]];f[i][j][k]%=mod;
                }return f[n1][n2][n3];
}
inline int ksm(int x,int k){
    int res=1;for(;k;k>>=1,x=x*x%mod) if(k&1) res=res*x%mod;return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n1=read();n2=read();n3=read();n=n1+n2+n3;prod[1]=1;
    m=read();mod=read();
    for(int i=1;i<=m;++i){
        for(int j=1;j<=n;++j) a[j]=read(),vis[j]=0;tot=0;
        for(int j=1;j<=n;++j) if(!vis[j]) cnt=0,dfs(j),w[++tot]=cnt;
        ans+=dp();ans%=mod;
    }++m;for(int i=1;i<=n;++i) w[i]=1;tot=n;ans+=dp();ans%=mod;
    printf("%d\n",ans*ksm(m,mod-2)%mod);
    return 0;
}