problem
求 n n n 个点带标号的简单无向连通图(无重边无自环)数目。
数据范围: n ≤ 130000 n \le 130000 n≤130000。
solution
由于是带了标号的,所以容易想到用 EGF(即指数生成函数)。
定义 G G G 为简单无向图(不一定连通)的 EGF,则:
G ( x ) = ∑ i = 0 n 2 C i 2 x i i ! G(x)=\sum_{i=0}^{n}2^{C_{i}^{2}}\frac{x^i}{i!} G(x)=i=0∑n2Ci2i!xi
显然 G G G 是可以直接计算的。
定义 F ( x ) F(x) F(x) 为简单无向连通的 EGF,则:
F ( x ) = ∑ i = 0 n f i x i i ! F(x)=\sum_{i=0}^{n}f_i\frac{x^i}{i!} F(x)=i=0∑nfii!xi
不难发现,最后答案是 f n f_n fn,也就是 [ x n ] F ( x ) × n ! [x^n]F(x)\times n! [xn]F(x)×n!,问题是如何计算 F ( x ) F(x) F(x)。
由于简单无向图是由若干简单无向连通图拼接而成的带标号的集合,这时想到生成集合,即:
s e t ( A ) = ∑ i = 0 ∞ A i i ! = e A set(A)=\sum_{i=0}^{\infty}\frac{A^i}{i!}=e^A set(A)=i=0∑∞i!Ai=eA
至于生成集合的原理我就不解释了,这里主要证一下 ∑ i = 0 ∞ A i i ! = e A \sum\limits_{i=0}^{\infty}\frac{A^i}{i!}=e^A i=0∑∞i!Ai=eA。
令 f ( x ) = e x f(x)=e^x f(x)=ex,我们把它在 0 0 0 处展开,得到:
f ( x ) = ∑ i = 0 ∞ f ( i ) ( 0 ) i ! x i f(x)=\sum_{i=0}^{\infty}\frac{f^{(i)}(0)}{i!}x^i f(x)=i=0∑∞i!f(i)(0)xi由于 e x e^x ex 的导函数还是 e x e^x ex,且 e 0 = 1 e^0=1 e0=1,所以:
e x = f ( x ) = ∑ i = 0 ∞ x i i ! e^x=f(x)=\sum_{i=0}^{\infty}\frac{x^i}{i!} ex=f(x)=i=0∑∞i!xi把 x x x 换成 A A A 代回去就可以了。
所以, G = s e t ( F ) = e F G=set(F)=e^F G=set(F)=eF,即 F = ln G F=\ln G F=lnG,那直接多项式 ln \ln ln 就可以解决此题。
code
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
#define N 500005
#define P 1004535809
using namespace std;
typedef vector<int> poly;
const int g=3;
int n,pos[N],fac[N],inv[N];
int add(int x,int y) {return x+y>=P?x+y-P:x+y;}
int dec(int x,int y) {return x-y< 0?x-y+P:x-y;}
int mul(int x,int y) {return 1ll*x*y%P;}
int power(int a,int b,int ans=1){
for(;b;b>>=1,a=mul(a,a))
if(b&1) ans=mul(ans,a);
return ans;
}
int *w[22],C=21;
void prework(){
for(int i=1;i<=C;++i)
w[i]=new int[1<<(i-1)];
int now=power(g,(P-1)/(1<<C));
w[C][0]=1;
for(int i=1;i<(1<<(C-1));++i) w[C][i]=mul(w[C][i-1],now);
for(int i=C-1;i;--i)
for(int j=0;j<(1<<(i-1));++j)
w[i][j]=w[i+1][j<<1];
}
void init(int lim){
for(int i=0;i<lim;++i)
pos[i]=(pos[i>>1]>>1)|((i&1)*(lim>>1));
}
void NTT(poly &f,int lim,int type){
for(int i=0;i<lim;++i)
if(pos[i]>i) swap(f[i],f[pos[i]]);
for(int mid=1,l=1;mid<lim;mid<<=1,++l){
for(int i=0;i<lim;i+=(mid<<1)){
for(int j=0;j<mid;++j){
int p0=f[i+j],p1=mul(f[i+j+mid],w[l][j]);
f[i+j]=add(p0,p1),f[i+j+mid]=dec(p0,p1);
}
}
}
if(type==-1&&(reverse(f.begin()+1,f.begin()+lim),1)){
int inv=power(lim,P-2);
for(int i=0;i<lim;++i) f[i]=mul(f[i],inv);
}
}
poly operator*(poly A,poly B){
int len=A.size()+B.size()-2,lim=1;
while(lim<=len) lim<<=1;init(lim);
A.resize(lim),NTT(A,lim,1);
B.resize(lim),NTT(B,lim,1);
for(int i=0;i<lim;++i) A[i]=mul(A[i],B[i]);
NTT(A,lim,-1),A.resize(len+1);
return A;
}
poly Inv(poly A,int len){
poly C,B(1,power(A[0],P-2));
for(int lim=4;lim<(len<<2);lim<<=1){
init(lim);
C=A,C.resize(lim>>1);
C.resize(lim),NTT(C,lim,1);
B.resize(lim),NTT(B,lim,1);
for(int i=0;i<lim;++i) B[i]=mul(B[i],dec(2,mul(B[i],C[i])));
NTT(B,lim,-1),B.resize(lim>>1);
}
B.resize(len);return B;
}
poly Deriv(poly A){
for(int i=0;i<A.size()-1;++i) A[i]=mul(A[i+1],i+1);
A.pop_back();return A;
}
poly Integ(poly A){
A.push_back(0);
for(int i=A.size()-1;i;--i) A[i]=mul(A[i-1],inv[i]);
A[0]=0;return A;
}
poly Ln(poly A,int len){
A=Integ(Deriv(A)*Inv(A,len)),A.resize(len);
return A;
}
int main(){
poly G(2,1),F;
scanf("%d",&n),prework();
fac[0]=fac[1]=1,inv[1]=1;
for(int i=2;i<N;++i) inv[i]=mul(P-P/i,inv[P%i]);
for(int i=2;i<=n;++i){
fac[i]=mul(fac[i-1],i);
G.push_back(mul(power(2,(1ll*i*(i-1)/2)%(P-1)),power(fac[i],P-2)));
}
int lim=1;
while(lim<(n<<1)) lim<<=1;
F=Ln(G,n+1);
printf("%d",mul(F[n],fac[n]));
return 0;
}