Luogu P4841 城市规划

题目大意

% 　　给定 $n$，求有 $n$ 个点的简单（无重边无自环）无向连通图的数目。
　　输出模 $479\ast 2^{21}+1=1004535809$ 意义下的答案。
　　数据范围　$1⩽n⩽130000$

题解

% 　　题中所给模数的原根为3。
　　定义 $f_n$ 为答案，即简单无向连通图的数量，$g_n$ 为无向图的数量，显然有 $g_0=g_1=1$。可以发现，对于任意 $n⩾2$，图中最多有 $C_n^2$ 条无向边，每条边可有可无，因而有 $$g_n=2^{C_n^2}\text{\hspace{1em}(1)}$$

% 　　若其中一个点所在的连通图大小为 $i$，则
$$g_n=\sum _{i=1}^n{C}_{n-1}^{i-1}{f}_i{g}_{n-i}\text{\hspace{1em}(2)}$$

% 　　把求和拆一下，得
$$g_n=f_n+\sum _{i=1}^{n-1}{C}_{n-1}^{i-1}{f}_i{g}_{n-i}$$

% 　　移项，组合数展开
$$\begin{array}{rl}{f}_n& =g_n-\sum _{i=1}^{n-1}\frac{(n-1)!}{(i-1)!(n-i)!}{f}_i{g}_{n-i}\\ & =g_n-(n-1)!\sum _{i=1}^{n-1}\frac{f_i}{(i-1)!}\frac{g_{n-i}}{(n-i)!}\end{array}$$

% 　　考虑
$$\begin{array}{r}{a}_i=\frac{f_i}{(i-1)!},b_i=\frac{g_i}{i!}\end{array}$$

% 　　因而 $f_n=(n-1)!a_n,g_n=i!b_i$，则有
$$(n-1)!a_n=n!b_n-(n-1)!\sum _{i=1}^{n-1}{a}_i{b}_{n-i}$$

% 　　同除 $(n-1)!$，得
$$a_n=n{b}_n-\sum _{i=1}^{n-1}{a}_i{b}_{n-i}$$

% 　　至此便可以用分治FFT处理了，时间复杂度 $\Theta (n{\log }_2^{2}n)$。

---

% 　　我们继续考虑 $(2)$ 式，将 $(1)$ 式代入 $(2)$ 式，展开，得
$$2^{C_n^2}=\sum _{i=1}^n\frac{(n-1)!}{(i-1)!(n-i)!}{f}_i{2}^{C_{n-i}^2}$$

% 　　同除 $(n-1)!$，得
$$\frac{2^{C_n^2}}{(n-1)!}=\sum _{i=1}^n\frac{f_i}{(i-1)!}\frac{2^{C_{n-i}^2}}{(n-i)!}$$

% 　　可以发现这是卷积的形式，定义
$$\begin{array}{rl}F(x)& =\sum _{i=1}^n\frac{f_i}{(i-1)!}{x}^i\\ G(x)& =\sum _{i=1}^n\frac{2^{C_i^2}}{i!}{x}^i\\ H(x)& =\sum _{i=1}^n\frac{2^{C_i^2}}{(i-1)!}{x}^i\end{array}$$

% 　　则有 $F\ast G=H$，因而有 $F=H\ast G^{-1}$，多项式求逆即可求出 $F$，再对答案乘 $(n-1)!$ 即可。时间复杂度为 $\Theta (n{\log }_{2}n)$。

代码

分治FFT？不存在的。

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define MAXN 460010
const long long mod=1004535809,G=3,inG=334845270;
inline long long pows(long long a,long long b){
	long long ret=1;
	while(b){
		if(b&1) ret=ret*a%mod;
		b>>=1;a=a*a%mod;
	} return ret;
}
#define invt(x) pows(x,mod-2)
class poly{
	public:
	static const int maxn=MAXN;
	long long *t;
	int limit,l,inmit;
	inline long long& operator[](int x){return t[x];}
	inline const long long& operator[](int x)const{return t[x];}
	inline poly(int n=0):t(NULL),limit(1),l(0){init(n);}
	inline void init(int n){
		if(limit>=n) return;
		if(t) free(t);
		while(limit<n) limit<<=1,++l;
		inmit=invt(limit);
		t=(long long*)calloc(limit+5,sizeof(long long));
	}
	void fnt(int type){
		static long long r[MAXN];
		static int now=0;
		if(now!=limit&&(now=limit))
			for(int i=1;i<limit;i++)
				r[i]=((r[i>>1]>>1)|((i&1)<<(l-1)));
		for(int i=1;i<limit;i++)
			if(i<r[i]) swap(t[i],t[r[i]]);
		for(int mid=1;mid<limit;mid<<=1){
			long long wn=pows(type==1?G:inG,(mod-1)/(mid<<1));
			for(int j=0,B=(mid<<1);j<limit;j+=B){
				long long w=1;
				for(int k=0;k<mid;k++,w=w*wn%mod){
					int x=t[j+k],y=w*t[j+k+mid]%mod;
					t[j+k]=(x+y)%mod;t[j+k+mid]=(x-y+mod)%mod;
				}
			}
		} if(type==-1) for(int i=0;i<limit;i++)
			t[i]=t[i]*inmit%mod;
	}
	inline void inv(int n,poly &g)const{
		g.init(n+n);
		if(n==1) return(void)(g[0]=invt(t[0]));
		inv(n+1>>1,g);
		poly F(n+n),G(n+n);
		for(int i=0;i<n;i++) F[i]=t[i];
		for(int i=0;i<n+1>>1;i++) G[i]=g[i];
		F.fnt(1);G.fnt(1);
		for(int i=0;i<F.limit;i++)
			F[i]=G[i]*((2-F[i]*G[i]%mod+mod)%mod)%mod;
		F.fnt(-1);
		for(int i=0;i<n;i++) g[i]=F[i];
	}
};
int n;
long long frac[MAXN],ifrac[MAXN];
int main(){
	scanf("%d",&n);++n;
	frac[0]=1;
	for(int i=1;i<=n;i++)
		frac[i]=frac[i-1]*i%mod;
	ifrac[n]=invt(frac[n]);
	for(int i=n-1;i>=0;i--)
		ifrac[i]=ifrac[i+1]*(i+1)%mod;
	poly G(n+n),H(n+n),inv_G(n+n);
	for(int i=1;i<n;i++)
		G[i]=ifrac[i]*pows(2,((long long)i*(i-1)/2))%mod;
	for(int i=1;i<n;i++)
		H[i]=ifrac[i-1]*pows(2,((long long)i*(i-1)/2))%mod;
	G[0]=1;G.inv(n,inv_G);
	inv_G.fnt(1);H.fnt(1);
	for(int i=0;i<H.limit;i++)
		H[i]=inv_G[i]*H[i]%mod;
	H.fnt(-1);
	printf("%lld\n",H[n-1]*frac[n-2]%mod);
	return 0;
}