【集训队互测2013】城市规划

Description

求n个点的带标号连通简单图的个数。答案对1004535809取模。
n<=130000

Solution

根据套路，首先我们需要求出n个点的带标号简单图的个数，设为g[n]
显然我们有g[n]=2^(n*(n-1)/2)
考虑设答案为f[n]，要怎么求出？
枚举1号点所在的联通块的大小，我们可以知道

$$g[n]=\sum_{i=1}^{n}f[i]C_{n-1}^{i-1}g[n-i]$$
因为f是保证联通的所以不会算重
把g带入，C拆开
$$2^{n*(n+1)/2}=\sum_{i=1}^{n}f[i]{(n-1)!\over {(i-1)!(n-i)!}}2^{(n-i)*(n-i-1)/2}$$
两边同除(n-1)!
$${2^{n*(n+1)/2}\over {(n-1)!}}=\sum_{i=1}^{n}{f[i]\over (i-1)!}{2^{(n-i)*(n-i-1)/2}\over (n-i)!}$$
可以发现这是一个卷积的形式
那么我们就直接上FFT就好了
需要多项式求逆
第一次写写的很丑QwQ

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;

typedef long long ll;
typedef double db;
ll read() {
    char ch;
    for(ch=getchar();ch<'0'||ch>'9';ch=getchar());
    int x=ch-'0';
    for(ch=getchar();ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
    return x;
}
void write(int x) {
    if (!x) {puts("0");return;}
    char ch[20];int tot=0;
    for(;x;x/=10) ch[++tot]=x%10+'0';
    fd(i,tot,1) putchar(ch[i]);
    puts("");
}

const int Mo=1004535809,D=3,N=5*1e5+5;
int pwr(int x,ll y) {
    int z=1;
    for(;y;y/=2,x=(ll)x*x%Mo)
        if (y&1) z=(ll)z*x%Mo;
    return z;
}

int n,len,lg,inv,t[N],a[N],b[N],c[N],W[N];
int p[N],h[N],g[N],f[N];
void pre(int len) {
    W[0]=1;W[1]=pwr(D,(Mo-1)/len);
    for(int i=2;i<=len;i++) W[i]=(ll)W[i-1]*W[1]%Mo;
    inv=pwr(len,Mo-2);
}
void DFT(int *a,int len,int lg,int flag) {
    for(int i=0;i<len;i++) {
        int p=0;
        for(int j=i,k=0;k<lg;k++,j/=2) p=(p<<1)+(j&1);
        t[p]=a[i];
    }
    for(int m=2;m<=len;m<<=1) {
        int half=m>>1,times=len/m;
        for(int i=0;i<half;i++) {
            int w=(flag==1)?W[i*times]:W[len-i*times];
            for(int j=i;j<len;j+=m) {
                int u=t[j],v=(ll)t[j+half]*w%Mo;
                t[j]=(u+v)%Mo;
                t[j+half]=(u-v+Mo)%Mo;
            }
        }
    }
    for(int i=0;i<len;i++) a[i]=t[i];
    if (flag==-1) for(int i=0;i<len;i++) a[i]=(ll)a[i]*inv%Mo;
}
void NTT(int *A,int *B,int *C,int len,int lg) {
    for(int i=0;i<len;i++) a[i]=A[i],b[i]=B[i];
    DFT(a,len,lg,1);DFT(b,len,lg,1);
    for(int i=0;i<len;i++) c[i]=(ll)a[i]*b[i]%Mo;
    DFT(c,len,lg,-1);
    for(int i=0;i<len;i++) C[i]=c[i];
}
void solve(int *A,int *B) {
    for(len=1,lg=0;len<=n*2;len*=2,lg++) {
        if (len==1) {B[0]=pwr(A[0],Mo-2);continue;}
        for(int i=0;i<len;i++) p[i]=A[i];
        pre(len*2);
        NTT(p,B,h,len*2,lg+1);
        for(int i=0;i<len*2;i++) h[i]=Mo-h[i];
        (h[0]+=2)%=Mo;
        NTT(h,B,B,len*2,lg+1);
        for(int i=len;i<len*2;i++) B[i]=0;
    } 
}
int fact[N],Inv[N],F[N],G[N],Inv_G[N],H[N];
int main() {
    scanf("%d",&n);

    fact[0]=1;fo(i,1,n) fact[i]=(ll)fact[i-1]*i%Mo;
    Inv[n]=pwr(fact[n],Mo-2);fd(i,n-1,0) Inv[i]=(ll)Inv[i+1]*(i+1)%Mo;

    fo(i,1,n) F[i]=(ll)pwr(2,(ll)i*(i-1)/2)*Inv[i-1]%Mo;
    fo(i,0,n) G[i]=(ll)pwr(2,(ll)i*(i-1)/2)*Inv[i]%Mo;
    solve(G,Inv_G);
    pre(len);NTT(F,Inv_G,H,len,lg);

    printf("%d\n",(ll)H[n]*fact[n-1]%Mo);
}

UPD：听说一次倍增只需要3次DFT QwQ
结果我常数直接*2压倒性劣势