bzoj 2601: [Jsoi2011]同分异构体计数

Description

Antonio 最近对有机化学比较感兴趣，他想请你帮助他快速计算出某种烃类的同分异

构体的数目。 

为了表述方便，我们作出如下定义： 

  环烷烃： 具有n 个碳原子的环烷烃可以表示成一张具有n 个顶点n 条边的无向连通

简单图(基环+外向树)。每个顶点的度数不超过 4。 

  M-环烷烃：至多有m 个顶点在环上的环烷烃。（注意环上至少有 3 个顶点，因为

任意两个顶点之间至多只能有1 条边）。 

 同构：假设结构A和结构B 均具有n 个碳原子，A和B 同构当且仅当能够对A和

B 中的每个碳原子都按照 1~n 编号，使得对于编号为 v1 和 v2 的两个碳原子，他

们在 A中存在边相连当且仅当他们在 B中存在边相连。（换言之，A和 B对应的图

同构）。 

现在，给出n, m，Antonio 希望你帮助他统计有多少种互不同构的含有n 个碳原子的

m-环烷烃。由于这个数量可能很大，你只需要输出它对p 的余数。（p是一个素数）。 

在本题中，我们不考虑某结构在化学上是否能够稳定存在，也不考虑其他的异构方式。

Input

输入文件只有一行，用空格隔开的三个整数n, m, p 。保证有m <=n,p为素数。

Output

输出文件有且仅有一行，表示具有n 个碳原子的互不同构的m-环烷烃的数量，对 p的

余数。

先处理出根的度为2，其余点度<=4的无标号有根树的方案数

环有旋转和翻转两种变换，由于m>=3，构成的置换群阶为2m，用burnside引理处理

旋转k(0<=k<m)步可以形成gcd(m,k)个等价类，每个等价类包含m/gcd(m,k)个位置

旋转+翻转需要分奇偶处理：

 若m为奇数，则有m个这种置换，形成(m+1)/2个等价类，其中一个等价类包含1个位置，其余包含2个位置

 若m为偶数

　　则有m/2个置换形成m/2个等价类，每个等价类包含2个位置

　　另有m/2个置换形成m/2+1个等价类，其中两个等价类包含1个位置，其余包含2个位置

#include<cstdio>
typedef unsigned long long u64;
typedef unsigned int u32;
int n,m;
u32 P;
int gcd(int a,int b){
    for(int c;b;c=a,a=b,b=c%b);
    return a;
}
int phi(int n){
    int v=n;
    for(int i=;i*i<=n;++i)if(n%i==){
        do n/=i;while(n%i==);
        v=v/i*(i-);
    }
    if(n>)v=v/n*(n-);
    return v;
}
inline u32 fix(int a){
    return a+(a>>&P);
}
struct num{
    u32 x;
    num(u32 a=):x(a){}
    num operator+(num w){return fix(x+w.x-P);}
    num operator*(num w){return u64(x)*w.x%P;}
    void operator+=(num w){x=fix(x+w.x-P);}
};
num s[][],gs[],iv[],f0[][],f1[][],ans;
void cal(int m,int n){
    int g=gcd(n,m);
    num v=;
    for(int d=;d<=g;++d)if(g%d==)v+=f0[m/d][n/d]*phi(d);
    v+=f1[m][n]*m;
    ans+=v*iv[m*];
}
int main(){
    scanf("%d%d%u",&n,&m,&P);
    if(m>n)m=n;
    s[][]=iv[]=;
    for(int i=;i<=;++i)iv[i]=iv[P%i]*(P-P/i);
    for(int i=;i<=n;++i){
        f0[][i]=f1[][i]=s[][i-]+s[][i-]+s[][i-];
        gs[]=f0[][i]+s[][i-];
        for(int j=;j<=;++j)gs[j]=gs[j-]*(gs[]+(j-))*iv[j];
        for(int j=;j;--j){
            for(int k=n;k>=i;--k){
                for(int t=;t<=j;++t){
                    int w=k-t*i;
                    if(w>=)s[j][k]+=gs[t]*s[j-t][w];
                }
            }
        }
    }
    for(int i=;i<=m;++i){
        for(int j=i;j<=n;++j){
            for(int k=;k<j;++k)f0[i][j]+=f0[i-][j-k]*f0[][k];
            if(i&){
                for(int k=;k<j;k+=)f1[i][j]+=f0[i>>][k>>]*f0[][j-k];
            }else{
                for(int k=;k<j;++k)f1[i][j]+=f1[i-][j-k]*f0[][k];
                if(~j&)f1[i][j]+=f0[i>>][j>>];
                f1[i][j]=f1[i][j]*iv[];
            }
        }
    }
    for(int i=;i<=m;++i)cal(i,n);
    printf("%d\n",ans.x);
    return ;
}