【XSY3032】画画（Burnside引理，计数）

为了方便，我们肯定是先考虑有标号图的个数，再用 Burnside 引理去重，但是用 Burnside 引理时得先考虑清楚映射集合 $X$ 是哪个集合 $A$ 到哪个集合 $B$ 的哪些映射，以及作用在 $A$ 上的置换群 $G$ 是什么。

首先考虑两张有标号图等价的定义：存在一个置换使得第一张图在置换作用后得到的图与第二张图完全相同，也就是说任意一条边 $(u,v)$ 要么在两者中同时存在，要么在两者中同时不存在。

那么自然地就能理清映射关系：$X$ 是由大小为 $\frac{n\times (n-1)}{2}$ 的边集 $A$ 到 B = { 0 , 1 } B=\{0,1\} B={0,1} 的所有映射，$0$ 表示这条边不存在，$1$ 表示这条边存在。

但是作用在 $A$ 上的置换群 $G$ 又是什么呢？题目中的置换是针对点的，貌似不是针对边的。

但我们知道，一种点置换对应着一种边置换，不同的点置换对应着不同的边置换，而且所有的点置换对应的所有的边置换也构成一个置换群。

于是 $G$ 就是所有的点置换对应的所有的边置换。

但事实上，由于所有的点置换和 $G$ 构成双射关系，我们可以通过枚举所有的点置换代替枚举 $G$ 中的边置换，并且用作用于图上的点置换代替作用于图上的边置换。

考虑 Burnside 引理：
$$\frac{1}{|G|}\sum _{g\in G}|X^g|$$
由刚刚我们说的双射关系，$|G|$ 显然为 $n!$。

然后我们枚举所有的点置换，可以将一个置换划分为若干个环，发现若由这些环的大小构成的可重集相同，那么两种点置换对答案的贡献是一样的，所以我们可以枚举这些可重集。

假设现在已经枚举了一个可重集，我们需要找到在当前情况下的不动点个数。（注意此时图上的点是有标号的）

我们先考虑环内的连边：若环的大小为 $2s$，那么有 $s-1$ 种连边方案不会影响点的奇偶性，有 $1$ 种方案会把这个环上所有的点的奇偶性都改变一次；若环的大小为 $2s+1$，则有 $s$ 种连边方案不会影响点的奇偶性。

再考虑环与环之间的连边：假设两个环 $A,B$ 的大小分别为 $a,b$，那么共有 $\gcd (a,b)$ 种连边方案，每一种连边方案都会给 $A$ 环中的每个点多连出 $\Delta A=\frac{\mathrm{lcm}(a,b)}{a}=\frac{b}{\gcd (a,b)}$ 条边，会给 $B$ 环中的每个点多连出 $\Delta B=\frac{\mathrm{lcm}(a,b)}{b}=\frac{a}{\gcd (a,b)}$ 条边。根据 $\Delta A$ 和 $\Delta B$ 的奇偶性有三种可能：两个环上的点都不改变奇偶性、某个环上的点改变奇偶性、两个环上的点都改变奇偶性。

发现无论怎么连边，环上的点改变奇偶性是一起改变的，所以我们可以把一个环缩为一个点。而对于环与环之间连边的最后一种可能，我们在两个环之间连一条边。

我们先把那些不会改变任何点的奇偶性的连边方案给统计上，现在问题变成了：给定一张图，点权最初全是 $0$，每个点 $i$ 有 $c{p}_i$ 次机会使点权异或 $1$，每条边 $i$ 有 $c{e}_i$ 次机会使两个端点的点权同时异或 $1$，问最后点权仍全是 $0$ 的方案数。

对每个连通块分开考虑，我们先给这个连通块（假设大小为 $s$）随便找一棵 dfs 树，然后给所有的树边留一次异或的机会（共留下 $s-1$ 次），然后假设连通块内除这 $s-1$ 次机会以外的其他所有的机会是否用都已经确定了，那么有合法方案当且仅当对于单点的异或次数之和为偶数，而且若有合法方案那么这 $s-1$ 次机会的使用情况也已经确定了，所以合法方案是唯一的。

所以这里贡献的方案数为 $2^{\max (\sum cp-1,0)+(\sum ce)-(s-1)}$，其中 $\sum cp-1$ 是因为前 $\sum cp-1$ 次机会任意选，然后根据奇偶确定最后一次机会是否选。

#include<bits/stdc++.h>

#define N 55

using namespace std;

namespace modular
{
    const int mod=998244353;
    inline int add(int x,int y){return x+y>=mod?x+y-mod:x+y;}
    inline int dec(int x,int y){return x-y<0?x-y+mod:x-y;}
    inline int mul(int x,int y){return 1ll*x*y%mod;}
}using namespace modular;

inline int poww(int a,int b)
{
    int ans=1;
    while(b)
    {
        if(b&1) ans=mul(ans,a);
        a=mul(a,a);
        b>>=1;
    }
    return ans;
}

inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-') f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+(ch^'0');
        ch=getchar();
    }
    return x*f;
}

int n,p2[N*N*2],fac[N],ifac[N],inv[N],gcd[N][N];
int cnt,a[N];
int fa[N],sp[N],se[N],size[N];
int ans;

int ggcd(int a,int b)
{
    return b?ggcd(b,a%b):a;
}

int find(int x)
{
    return x==fa[x]?x:(fa[x]=find(fa[x]));
}

int calc()
{
    int ans=1,tim=0;
    for(int i=1;i<=cnt;i++)
    {
        ans=mul(ans,inv[a[i]]);
        tim++;
        if(a[i]!=a[i+1]||i==cnt)
        {
            ans=mul(ans,ifac[tim]);
            tim=0;
        }
    }
    for(int i=1;i<=cnt;i++) fa[i]=i,sp[i]=se[i]=0,size[i]=1;
    for(int i=1;i<=cnt;i++)
    {
        ans=mul(ans,p2[(a[i]-1)>>1]);
        if(!(a[i]&1)) sp[i]++;
    }
    for(int i=1;i<=cnt;i++)
    {
        for(int j=i+1;j<=cnt;j++)
        {
            int g=gcd[a[i]][a[j]],A=a[j]/g,B=a[i]/g;
            if(!(A&1)&&!(B&1)) ans=mul(ans,p2[g]);
            else if((A&1)&&(B&1))
            {
                int x=find(i),y=find(j);
                se[x]+=g;
                if(x!=y)
                {
                    fa[y]=x;
                    sp[x]+=sp[y];
                    se[x]+=se[y];
                    size[x]+=size[y];
                }
            }
            else if(A&1) sp[find(i)]+=g;
            else sp[find(j)]+=g;
        }
    }
    for(int i=1;i<=cnt;i++)
        if(find(i)==i)
            ans=mul(ans,p2[max(sp[i]-1,0)+se[i]-(size[i]-1)]);
    return ans;
}

void dfs(int sum,int lst)
{
    if(sum==n)
    {
        ans=add(ans,calc());
        return;
    }
    for(int i=lst;i<=n-sum;i++)
    {
        a[++cnt]=i;
        dfs(sum+i,i);
        cnt--;
    }
}

int main()
{
    n=read();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            gcd[i][j]=ggcd(i,j);
    p2[0]=1;
    for(int i=1;i<=n*n*2;i++) p2[i]=add(p2[i-1],p2[i-1]);
    fac[0]=1;
    for(int i=1;i<=n;i++) fac[i]=mul(fac[i-1],i);
    ifac[n]=poww(fac[n],mod-2);
    for(int i=n;i>=1;i--) ifac[i-1]=mul(ifac[i],i);
    for(int i=1;i<=n;i++) inv[i]=poww(i,mod-2);
    dfs(0,1);
    printf("%d\n",ans);
    return 0;
}