poj1737 Connected Graph（计数，组合数学，递推，高精）

求n个点的连通图的个数f[n]，带标号。有两种做法。
1.直接搞出f[n]的递推式。
$f[n]=\sum\limits_{k=1}^{n-1}f[k]*f[n-k]*C_{n-2}^{k-1}*(2^k-1)$
如何理解呢？
对于一张连通图，我们首先拿掉1号点，考虑此时2号点的情况，假设此时有k-1个点仍与2号点相连（方案数为$C_{n-2}^{k-1}$），则我们把原图视作两个连通块，这k个点视作一个连通块（方案数为f[k]），剩余的n-k个点视作连通块（方案数为f[n-k]）。则这两个连通块一定要联通，且只能通过1号点来联通。于是有$2^k-1$种连边方案。所以方案数就是$f[k]*f[n-k]*C_{n-2}^{k-1}*(2^k-1)$
写一下高精就好了。

2.利用补集思想,用n个点的所有图个数h[n]-所有不联通的图个数g[n]

首先显然有$h[n]=2^{C_n^2}$

如何计算g[n]呢？

我们假设1所在连通块有k个点，剩余n-k个点随意，则方案数就是$C_{n-1}^{k-1}*f[k]*h[n-k]$，所以我们有
$g[n]=\sum\limits_{k=1}^{n-1}C_{n-1}^{k-1}*f[k]*h[n-k]$

不过这种写法还得写高精减，比较麻烦。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 400
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n;
ll bin[60],C[60][60];
struct bigint{
    ll a[N];int n;
    bigint(){memset(a,0,sizeof(a));n=0;}
    ll& operator[](int x){return a[x];}
    inline void print(){
        for(int i=n;i>=1;--i) printf("%lld",a[i]);puts("");
    }
    friend bigint operator+(bigint x,bigint y){
        bigint res;res.n=max(x.n,y.n);
        for(int i=1;i<=res.n;++i) res[i]=x[i]+y[i];
        for(int i=1;i<=res.n;++i) res[i+1]+=res[i]/10,res[i]%=10;
        if(res[res.n+1]) res.n++;return res;
    }friend bigint operator*(bigint x,bigint y){
        bigint res;res.n=x.n+y.n-1;
        for(int i=1;i<=x.n;++i)
            for(int j=1;j<=y.n;++j)
                res[i+j-1]+=x[i]*y[j];
        for(int i=1;i<=res.n;++i) res[i+1]+=res[i]/10,res[i]%=10;
        while(res[res.n+1]) res.n++,res[res.n+1]+=res[res.n]/10,res[res.n]%=10;
        return res;
    }friend bigint operator*(bigint x,ll y){
        bigint res;res.n=x.n;
        for(int i=1;i<=x.n;++i) res[i]=y*x[i];
        for(int i=1;i<=res.n;++i) res[i+1]+=res[i]/10,res[i]%=10;
        while(res[res.n+1]) res.n++,res[res.n+1]+=res[res.n]/10,res[res.n]%=10;
        return res;
    }
}f[60];
int main(){
//  freopen("a.in","r",stdin);
    for(int i=0;i<=48;++i) C[i][0]=1;bin[0]=1;
    for(int i=1;i<=50;++i) bin[i]=bin[i-1]<<1;  
    for(int i=1;i<=48;++i)
        for(int j=1;j<=i;++j)
            C[i][j]=C[i-1][j]+C[i-1][j-1];
    f[1].n=1;f[1][1]=1;f[2].n=1;f[2][1]=1;
    for(int i=3;i<=50;++i)
        for(int j=1;j<i;++j)
            f[i]=f[i]+f[j]*f[i-j]*C[i-2][j-1]*(bin[j]-1);
    while(1){
        int x=read();if(!x) break;f[x].print();
    }return 0;
}