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; }