【BZOJ3456】城市规划

Description

刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了.
刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通. 为了省钱, 每两个城市之间最多只能有一条直接的贸易路径. 对于两个建立路线的方案, 如果存在一个城市对, 在两个方案中是否建立路线不一样, 那么这两个方案就是不同的, 否则就是相同的. 现在你需要求出一共有多少不同的方案.
好了, 这就是困扰阿狸的问题. 换句话说, 你需要求出n个点的简单(无重边无自环)无向连通图数目.
由于这个数字可能非常大, 你只需要输出方案数mod 1004535809(479 * 2 ^ 21 + 1)即可.

Input

仅一行一个整数n(<=130000)

Output

仅一行一个整数, 为方案数 mod 1004535809.

Sample Input

Sample Output

HINT

对于 100%的数据, n <= 130000

题解：一直打算做来着，前几天终于抽空把这题过了~

设f[n]表示n个点的简单无向连通图数目，先列出大家都会的DP方程：

$f[n]=2^{n(n-1)\over 2}-\sum\limits_{i=1}^{n-1}f[i]C_{n-1}^{i-1}2^{(n-i)(n-i-1)\over 2}$

哎呀好多n-1太烦了，还是令f[n]表示n+1个点的简单无向连通图数目吧：

$f[n]=2^{n(n+1)\over 2}-\sum\limits_{i=0}^{n-1}f[i]C_n^i2^{(n-i)(n-i-1)\over 2}$

组合数很不优美是不是？那就拆了好啦。

$f[n]=2^{n(n+1)\over 2}-n!\sum\limits_{i=0}^{n-1}{f[i]2^{(n-i)(n-i-1)\over 2}\over i!(n-i)!}$

下面我们想把只含有含有n的放到一起，i放到一起，n-i放到一起，于是疯狂移项得到：

${f[n]\over n!}={2^{n(n+1)\over 2}\over n!}-\sum\limits_{i=0}^{n-1}{f[i]\over i!}\cdot{2^{(n-i)(n-i-1)\over 2}\over (n-i)!}$

明朗很多了是不是？令$F(x)={f[x]\over x!},G(x)={2^{x(x+1)\over 2}\over x!},H(x)={2^{x(x-1)\over 2}\over x!}$，特别地，H(0)=0。

所以$F(x)=G(x)-F(x)*H(x)$，$F(x)={G(x)\over 1+H(x)}$，多项式求个逆就好了。

#include <cstdio>

#include <cstring>

#include <iostream>

#include <algorithm>

using namespace std;

typedef long long ll;

const ll P=1004535809;

const int maxn=(1<<19)+4;

int n;

ll A[maxn],B[maxn],C[maxn],D[maxn],ine[maxn],jc[maxn],jcc[maxn];

inline ll pm(ll x,ll y)

{

	ll z=1;

	while(y)

	{

		if(y&1)	z=z*x%P;

		x=x*x%P,y>>=1;

	}

	return z;

}

inline void NTT(ll *a,int len,int f)

{

	int i,j,k,h;

	for(i=k=0;i<len;i++)

	{

		if(i>k)	swap(a[i],a[k]);

		for(j=len>>1;(k^=j)<j;j>>=1);

	}

	ll wn,w,t;

	for(h=2;h<=len;h<<=1)

	{

		if(f==1)	wn=pm(3,(P-1)/h);

		else	wn=pm(3,P-1-(P-1)/h);

		for(i=0;i<len;i+=h)

		{

			w=1;

			for(j=i;j<i+h/2;j++)	t=w*a[j+h/2]%P,a[j+h/2]=(a[j]-t)%P,a[j]=(a[j]+t)%P,w=w*wn%P;

		}

	}

	if(f==-1)

	{

		t=pm(len,P-2);

		for(i=0;i<len;i++)	a[i]=a[i]*t%P;

	}

}

void get_inv(ll *F,ll *G,int len)

{

	if(len==1)

	{

		G[0]=pm(F[0],P-2);

		return ;

	}

	int i;

	get_inv(F,G,len>>1);

	memcpy(C,F,sizeof(F[0])*len);

	memset(C+len,0,sizeof(F[0])*len);

	NTT(C,len<<1,1),NTT(G,len<<1,1);

	for(i=0;i<len<<1;i++)	G[i]=(G[i]*2-C[i]*G[i]%P*G[i])%P;

	NTT(G,len<<1,-1);

	memset(G+len,0,sizeof(F[0])*len);

}

int main()

{

	scanf("%d",&n);

	int len=1,i;

	while(len<=n)	len<<=1;

	ine[0]=ine[1]=jc[0]=jc[1]=jcc[0]=jcc[1]=1;

	for(i=2;i<=len;i++)	ine[i]=P-(P/i)*ine[P%i]%P,jc[i]=jc[i-1]*i%P,jcc[i]=jcc[i-1]*ine[i]%P;

	for(i=0;i<len;i++)	D[i]=pm(2,i*(i+1ll)/2)*jcc[i]%P;

	for(i=1;i<len;i++)	A[i]=pm(2,i*(i-1ll)/2)*jcc[i]%P;

	A[0]=1;

	get_inv(A,B,len);

	NTT(D,len<<1,1),NTT(B,len<<1,1);

	for(i=0;i<len<<1;i++)	B[i]=D[i]*B[i]%P;

	NTT(B,len<<1,-1);

	printf("%lld",(B[n-1]*jc[n-1]%P+P)%P);

	return 0;

}