BZOJ3456 : 城市规划

设

$a[i]=2^{C_i^2}$

$b[i]=i!$

$c[i]=(i-1)!$

f[i]为i个点的简单无向连通图数目，则有

\[\begin{eqnarray*}f[i]&=&a[i]-\sum_{j=1}^{i-1}C_{i-1}^{j-1}f[j]a[i-j]\\
&=&a[i]-c[i]\sum_{j=1}^{i-1}\frac{f[j]a[i-j]}{c[j]b[i-j]}\end{eqnarray*}\]

设

$A[i]=\frac{f[i]}{c[i]}$

$B[i]=\frac{a[i]}{b[i]}$

则

\[\begin{eqnarray*}f[i]=a[i]-c[i]\sum_{j=1}^{i-1}A[j]B[i-j]\end{eqnarray*}\]

右边是个卷积的形式，可以通过分治+NTT在$O(n\log^2n)$复杂度内解决。

 

#include<cstdio>

typedef long long ll;

const int N=131072,K=16;

int n,i,a[N],b[N],c[N],nb[N],nc[N],f[N],A[N],B[N],P=1004535809,g[K+1],ng[K+1],inv[N+1];

inline int pow(int a,ll b){int t=1;for(;b;b>>=1,a=(ll)a*a%P)if(b&1)t=(ll)t*a%P;return t;}

inline void NTT(int*a,int n,int t){

  for(int i=1,j=0;i<n-1;i++){

    for(int s=n;j^=s>>=1,~j&s;);

    if(i<j){int k=a[i];a[i]=a[j];a[j]=k;}

  }

  for(int d=0;(1<<d)<n;d++){

    int m=1<<d,m2=m<<1,_w=t==1?g[d]:ng[d];

    for(int i=0;i<n;i+=m2)for(int w=1,j=0;j<m;j++){

      int&A=a[i+j+m],&B=a[i+j],t=(ll)w*A%P;

      A=(B-t+P)%P;B=(B+t)%P;w=(ll)w*_w%P;

    }

  }

  if(t==-1)for(int i=0,j=inv[n];i<n;i++)a[i]=(ll)a[i]*j%P;

}

void solve(int l,int r){

  if(l==r){f[l]=((a[l]-(ll)c[l]*f[l]%P)%P+P)%P;return;}

  int mid=(l+r)>>1,i,n=r-l+1,k;

  solve(l,mid);

  for(k=1;k<n;k<<=1);

  for(i=0;i<k;i++)A[i]=i+l<=mid?(ll)f[i+l]*nc[i+l]%P:0,B[i]=(ll)a[i]*nb[i]%P;

  NTT(A,k,1),NTT(B,k,1);

  for(i=0;i<k;i++)A[i]=(ll)A[i]*B[i]%P;

  NTT(A,k,-1);

  for(i=r;i>mid;i--)f[i]=(f[i]+A[i-l])%P;

  solve(mid+1,r);

}

int main(){

  for(g[K]=pow(3,(P-1)/N),ng[K]=pow(g[K],P-2),i=K-1;~i;i--)g[i]=(ll)g[i+1]*g[i+1]%P,ng[i]=(ll)ng[i+1]*ng[i+1]%P;

  for(inv[1]=1,i=2;i<=N;i++)inv[i]=(ll)(P-inv[P%i])*(P/i)%P;

  scanf("%d",&n);

  for(b[0]=c[0]=nb[0]=nc[0]=i=1;i<=n;i++){

    a[i]=pow(2,(ll)i*(i-1)/2),b[i]=(ll)b[i-1]*i%P,c[i]=b[i-1];

    nb[i]=pow(b[i],P-2),nc[i]=nb[i-1];

  }

  solve(1,n);

  printf("%d",f[n]);

  return 0;

}