gxx_slide之城市规划

链接：http://oj.tsinsen.com/A1493
题意：求n个点的连通无向图的种类数$(n<=130000)$
分析：n比较小的版本在做递推专题的时候是碰到过的,思考方法就是枚举第一个点所在连通块的点数，点数为n就是我们要求的。
不妨记$G_i为i个点图的种类数，那么显然G_i=2^{C_i^2}$
记F_i为我们要求的答案，那么根据前面所说，我们有
$G_n=\sum_{i=1}^{n}C_{n-1}^{i-1}*F_i*G_{n-i}$
$F_n=G_n-\sum_{i=1}^{n-1}C_{n-1}^{i-1}*F_i*G_{n-i}$
仿照前几题的方式，右边可以转化为一个卷积形式，由于是一项一项求，我用cdq分治，复杂度为$O(n*log^2n)$
当然，论文里给出了一个$O(n*logn)$的解答，需要用到多项式逆的知识（明天去学一下，还想学矩阵乘法的递推优化，但是下个礼拜要去邀请赛了，字符串还要加强一下，还要期末考，不知道来不来得及。。。。）

#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
const int M=1004535809,g=3;
typedef long long Int;
const int Maxn=130010;
Int fac[Maxn],inv[Maxn],G[Maxn],F[Maxn];
Int powmod(Int x,Int y,Int mod)
{
    Int ret=1,t=x;
    while(y){if(y&1)ret=ret*t%mod;t=t*t%mod;y>>=1;}
    return ret;
}
Int a[162144],b[162144];
void rev(Int *a,int N)
{
    int i,j,k;
    for(i=1,j=N>>1;i<N-1;i++)
    {
        if(i<j)swap(a[i],a[j]);
        for(k=N>>1;j>=k;j-=k,k>>=1);j+=k;
    }
}
void dft(Int *a,int N,int flag=1)
{
    rev(a,N);
    for(int m=2;m<=N;m<<=1)
    {
        Int wm=powmod(g,(M-1)/m,M);
        if(flag<0)wm=powmod(wm,M-2,M);
        for(int k=0;k<N;k+=m)
        {
            Int w=1;
            for(int j=k;j<k+(m>>1);j++,w=w*wm%M)
            {
                Int u=a[j],v=a[j+(m>>1)]*w%M;
                a[j]=(u+v)%M;
                a[j+(m>>1)]=(u-v+M)%M;
            }
        }
    }
}
void mul(Int *a,Int *b,int N)
{
    dft(a,N);dft(b,N);
    for(int i=0;i<N;i++)a[i]=a[i]*b[i]%M;
    dft(a,N,-1);
    Int revN=powmod(N,M-2,M);
    for(int i=0;i<N;i++)a[i]=a[i]*revN%M;
}
void solve(int l,int r)
{
    if(l==r)return;
    int mid=(l+r)>>1;
    solve(l,mid);
    int N=1;
    while(N<r-l)N<<=1;
    for(int i=0;i<=mid-l;i++)a[i]=F[i+l]*(i+l)%M*inv[i+l]%M;
    for(int i=mid-l+1;i<N;i++)a[i]=0;
    for(int i=0;i<r-l;i++)b[i]=G[i+1]*inv[i+1]%M;
    for(int i=r-l;i<N;i++)b[i]=0;
    mul(a,b,N);
    for(int i=mid+1;i<=r;i++)
    {
        F[i]=(F[i]-fac[i-1]*a[i-l-1]%M+M)%M;
    }
    solve(mid+1,r);
}
int main()
{
    fac[0]=fac[1]=inv[0]=inv[1]=1;
    for(int i=2;i<Maxn;i++)fac[i]=fac[i-1]*i%M,inv[i]=inv[M%i]*(M-M/i)%M;
    for(int i=2;i<Maxn;i++)inv[i]=inv[i-1]*inv[i]%M;
    for(int i=0;i<Maxn;i++)G[i]=powmod(2,i*(Int)(i-1)/2%(M-1),M);
    for(int i=0;i<Maxn;i++)F[i]=G[i];
    solve(1,130000);
    int n;
    while(scanf("%d",&n)!=EOF)printf("%lld\n",F[n]);
    return 0;
}