HNOI 2009 图的同构记数 题解

题目传送门

题目大意： 问有多少个不同构的 $n$ 个点的图。

题解

A图的顶点经过一定的重新标号以后，A图的顶点集和边集要完全与B图一一对应。这是题目里告诉我们的同构的条件，这应该算半个提示了，我们可以认为，重新标号是一种置换，这种置换形如：
$$\left(\begin{array}{ccccc}1& 2& ...& n-1& n\\ a_1& a_2& ...& a_{n-1}& a_n\end{array}\right)$$

那么总的置换数量就是 $|G|=n!$，然后用burnside定理来求，也就是要求出每种置换有多少个不动点。

先将一个置换表示成若干个循环的乘积，然后考虑连边，假如循环 $A$ 内第 $x$ 个元素和第 $y$ 个元素之间连了边，因为要在置换后图保持不变，所以第 $x+1$ 和第 $y+1$ 个元素之间也要有连边，同理，距离为 $|x-y|$ 的元素之间都需要连边。

可以发现，对于一个长度 $d$，要么全部连边，要么全都不连，而 $d$ 的取值范围是 $[1,\lfloor \frac{L_i}{2}\rfloor ]$，其中 $L_A$ 是这个循环的长度，则一共有 $2^{\lfloor \frac{L_i}{2}\rfloor }$ 种连边方案。

再考虑两个不同的循环 $A$ 和 $B$，类似的，假如 $A$ 中的第 $i$ 个元素和 $B$ 中的第 $j$ 个元素之间连了边，那么 $i+1$ 和 $j+1$ 之间也要连，以此类推下去要连 $lcm(L_A,L_B)$ 条边，不妨称这样的 $(i,j),(i+1,j+1)$ 为同一组边。

因为一共有 $L_A\times L_B$ 条边，所以一共有 $\frac{L_A\times L_B}{lcm(,L_A),L_B}=\gcd (L_A,L_B)$ 组边，每组边也是要么连要么不连，方案数为 $2^{\gcd (L_A,L_B)}$。

所以答案就是：
$$\frac{1}{|G|}\sum _{L_1,L_2,...,L_k}\left(\prod _{i=1}^k{2}^{\lfloor \frac{L_i}{2}\rfloor }\times \prod _{1\le A<B\le k}{2}^{\gcd (L_A,L_B)}\right)$$

但是这样的话时间复杂度为 $O(n!\times n^2)$，可以愉快地得到 $TLE$ 的好成绩。

考虑优化，观察后可以发现，我们并不关心每个循环内的点是谁，只关心这个置换能拆成几个循环以及每个循环的大小 $L_i$。

于是考虑直接枚举 $L_i$，相当于将 $n$ 进行自然数拆分，经暴力发现，当 $n=60$ 时，拆分方案大概只有 $10^6$ 种，可以接受。

当我们枚举出所有 $L_i$ 之后，就可以用上面的柿子求解，但是还需要求多一个东西：有多少个置换拆成循环之后，循环的长度分别是 $L_1,L_2,...,L_k$。

这个很简单，排列组合一下就能得到：

先给这 $k$ 个循环安排元素进去，方案数为 $\frac{n!}{{\prod }_{i=1}^k{L}_i!}$。

再考虑每个循环内的元素顺序，不考虑重复的话就有 $L_i!$ 排列，但是对于一个循环而言，$(a_1{a}_2{a}_3...a_n)$ 和 $(a_2{a}_3{a}_4...a_n{a}_1)$ 其实是同样的排列，所以要除以 $L_i$，即方案数为 $\frac{L_i!}{L_i}=(L_i-1)!$。

然后对于长度相同的两个循环，比如说 $(12)(34)$，它其实等价于 $(34)(12)$，设长度为 $i$ 的循环有 $S_i$ 个，那么还需要除以 ${\prod }_{i=1}^k{S}_i!$。

所以一共有 $\frac{n!}{{\prod }_{i=1}^k{L}_i\times {\prod }_{i=1}^k{S}_i!}$ 个置换满足要求。

枚举出 $L$ 之后，把这个方案数和上面那个柿子乘起来就是答案了。以及，这里的 $n!$ 和上面的 $\frac{1}{|G|}$ 可以抵消掉。

设对 $n$ 进行自然数拆分的时间复杂度为 $K$，那么最后的时间复杂度为 $O(K\times n^2)$。

讲道理，当 $n=60$ 的时候时间复杂度大概为 $3.6\times 10^9$，虽然说完全跑不满，但是实测也需要 $3s$ 左右，为什么能过呢……

太烧脑了果然还是不管好了qwq，代码如下：

#include <cstdio>
#define mod 997
#define maxn 110

int n;
int ksm(int x,int y)
{
	int re=1;
	while(y)
	{
		if(y&1)re=1ll*re*x%mod;
		x=1ll*x*x%mod;y>>=1;
	}
	return re;
}
int fac[maxn],inv_fac[maxn],inv[maxn],gcd[maxn][maxn];
int gcd_(int x,int y){return y==0?x:gcd_(y,x%y);}
void work()
{
	fac[0]=inv_fac[0]=inv[1]=1;
	for(int i=1;i<=n;i++)fac[i]=1ll*fac[i-1]*i%mod;
	inv_fac[n]=ksm(fac[n],mod-2);
	for(int i=n-1;i>=1;i--)inv_fac[i]=1ll*inv_fac[i+1]*(i+1)%mod;
	for(int i=2;i<=n;i++)inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)gcd[i][j]=gcd_(i,j);
}
int L[maxn],S[maxn],ans=0;
void work(int k)
{
	int prod=1;
	for(int i=1;i<=k;i++)prod=1ll*prod*((1ll<<(L[i]/2))%mod)%mod;
	for(int i=1;i<=k;i++)for(int j=i+1;j<=k;j++)prod=1ll*prod*(1ll<<gcd[L[i]][L[j]])%mod;
	for(int i=1;i<=k;i++)prod=1ll*prod*inv[L[i]]%mod;
	for(int i=1;i<=n;i++)prod=1ll*prod*inv_fac[S[i]]%mod;
	ans=(ans+prod)%mod;
}
int min(int x,int y){return x<y?x:y;}
void dfs(int x,int k)
{
	if(!x)return (void)work(k-1);
	for(int i=min(L[k-1],x);i>=1;i--)
	L[k]=i,S[i]++,dfs(x-i,k+1),S[i]--;
}

int main()
{
	scanf("%d",&n);work();
	L[0]=n;dfs(n,1);
	printf("%d",ans);
}