Poj1737 Connected Graph 计数dp、组合数、高精度

又是一道令人胆寒的组合数。。。。

【一句话题意】求n个结点的无向图的个数。（每个节点有标号）n<=50
【分析】对于dp来说，状态的划分和转移是十分重要的。一个连通图是难以进行状态划分的，但一个不连通的无向图是则很容易划分成结点更少的两个部分。所以我们考虑从所有的无向图中减去所有的不连通图，就是所有连通图的数量。
一个无向联通图中至多有$n\ast (n-1)/2$条边，每条边均可有可无，所以共有2^n*(n-1)/2种无向图。

组合数是十分注重“不重不漏”的严谨思维的。所以下文且先划上重点标识符。计算n个点的不连通无向图数量。利用围绕基准点构造一个整体的思想，我们可以枚举标号为1的结点所在的连通块包含的结点个数k，从2~n这n-1个结点中选出k-1个，与1号结点一起构成了大小为k的连通块。显然，我们有$C_{n-1}^{k-1}$种选法。剩余n-k个结点构成任意无向图有2^{(n-k)*(n-k-1)/2}种方法。
根据以上分析不难得出状态转移方程：
$f[i]=2^{i\ast (i-1)/2}-{\Sigma }_{j=1}^{i-1}f[j]\ast C_{i-1}^{j-1}\ast 2^{(i-j)\ast (i-j-1)/2}$
初值f[1]=1;

然而这并没有卵用，因为你还需要高精度乘法、加法。。。。
然后就可以离线做了。
【code】

#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(){
    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;
}