ii个点的连通块个数为F(i)F(i),所有ii个点的图的个数为G(i)G(i)
从反面考虑问题,发现ii个点连通图个数=所有ii个点的图的个数-ii个点非连通图的个数。
考虑所有ii个点的图的个数,每条边可选可不选,共i×(i1)/2i \times (i-1)/2条边,则G(i)=2i×(i1)/2G(i)=2^{i \times (i-1)/2}
考虑ii个点非连通图的个数,设节点11所在的连通块大小为jj,那么我们选出这个连通块的方法数为F(j)×Ci1j1F(j) \times C^{j-1}_{i-1},也就是从剩下i1i-1个节点中选出j1j-1个,与节点11组成连通块。
剩下iji-j个点随便乱连,有G(ij)G(i-j)种搞法。

那么我们可以推出dpdp方程式:

F(i)=G(i)F(j)×Ci1j1×G(ij)(1ji1)F(i)=G(i)-\sum F(j) \times C^{j-1}_{i-1} \times G(i-j) (1 \le j \le i-1)

代码(高精自己打)

bign C[51][51];
inline void Init(){
    for (register int i=1;i<=50;++i){
        C[i][0]=1,C[i][i]=1;
        for (register int j=1;j<i;++j){
            C[i][j]=C[i-1][j-1]+C[i-1][j];
        }
    }
}
inline bign pow2(int k){
    bign ans=1,p=2;
    while (k){
        if (k&1) ans=ans*p;
        p=p*p;
        k>>=1;
    }
    return ans;
}
bign dp[51];
inline bign f(int x){//x个点随便怎么连
    return pow2(x*(x-1)/2);
}
int main() {
    Init();
    dp[1]=1;
    for (register int i=2;i<=50;++i){
        dp[i]=f(i);
        for (register int j=1;j<=i-1;++j){
            dp[i]=dp[i]-dp[j]*C[i-1][j-1]*f(i-j);
        }
    }
    while (true){
        int x;
        scanf("%d",&x);
        if (x==0) break;
        cout<<dp[x]<<endl;
    }
}