poj1737_带标号连通图计数

看到这题总有一种不好的联想， 感觉数学考试出这个题绝对不会做......(有厉害的同桌也木有用，因为这个题目是他当数学题问我的)

如果f(n)就是答案的话， 那么考虑另一个g(n)， 表示n个点的非连通图， 那么显然由完全图的基本性质可得f(n) + g(n) = 2 ^ (n*(n-1)/2) = h(n)； 故f(n)有了新的表现形式；

对于g(n)的计算， 考虑用组合数学方法： 在n个点中选定一个点， 讨论它所属的连通分量中元素的个数， 当有i个元素时， 这样的元素组合的个数共有C(n-1, i-1)个， 同时又有f(i)种连接方式使其联通。 在原图中去掉当前分量， 剩下的即为h(n-i)， 顿时感觉豁然开朗——递推正在招手！ 

能得到式子： g(n) = Σ(C(n-1,i-1)*f(i)*h(n-i))；

呵呵了， 记忆化搜索， 问题圆满解决。

不爽： 高精度占的版面简直过量！！！

#include <cstdio>
#include <cstring>
#define N 50 + 3
#define mod 10000

using namespace std;

struct node
{
    int len, l[1000];
    void clear() { memset(l, 0, sizeof l); len = 1; }
    void print()
    {
        printf("%d", l[len]);
        for (int i = len - 1; i; i--)
        printf("%04d", l[i]);
        putchar('\n');
    }
}f[N], g[N], h[N];
int n;
bool flag[3][N];
inline int max(int x, int y)
{
    return x > y ? x : y;
}
node mul(node a, node b)
{
    node c;
    c.clear();
    c.len = max(a.len, b.len);
    for (int i = 1; i <= a.len; ++i)
    for (int j = 1; j <= b.len; ++j)
    c.l[i+j-1] += a.l[i] * b.l[j];
    for (int i = 1; i <= c.len; ++i)
    {
        c.l[i+1] += c.l[i] / mod;
        c.l[i] %= mod;
    }
    while(c.l[c.len+1])
    {
        ++c.len;
        c.l[c.len+1] += c.l[c.len] / mod;
        c.l[c.len] %= mod;
    }
    return c;
}
node minus(node a, node b)
{
    for (int i = 1; i <= a.len; ++i)
    {
        a.l[i] -= b.l[i];
        while(a.l[i] < 0)
        {
            a.l[i] += mod;
            a.l[i+1]--;
        }
    }
    while(!a.l[a.len] && a.len != 1) a.len--;
    return a;
}
node plus(node a, node b)
{
    a.len = max(a.len, b.len);
    for (int i = 1; i <= a.len; ++i)
    {
        a.l[i] += b.l[i];
        a.l[i+1] += a.l[i] / mod;
        a.l[i] %= mod;
    }
    if (a.l[a.len+1]) ++a.len;
    return a;
}
node pow(int x)
{
    node ans, a;
    ans.clear(); a.clear();
    ans.l[1] = 1; a.l[1] = 2;
    while(x)
    {
        if (x & 1) ans = mul(ans, a);
        a = mul(a, a);
        x >>= 1;
    }
    return ans;
}
void mul(node &a, int b)
{
    for (int i = 1; i <= a.len; ++i)
    a.l[i] *= b;
    for (int i = 1; i <= a.len; ++i)
    {
        a.l[i+1] += a.l[i] / mod;
        a.l[i] %= mod;
    }
    while(a.l[a.len+1])
    {
        ++a.len;
        a.l[a.len+1] += a.l[a.len] / mod;
        a.l[a.len] %= mod;
    }
}
void div(node &a, int b)
{
    int now = 0, p = 0;
    for (int i = a.len; i; i--)
    {
        (now *= mod) += a.l[i];
        a.l[i] = 0;
        if (now < b) continue;
        if (!p)
        {
            p = 1;
            a.len = i;
        }
        a.l[i] = now / b;
        now %= b;
    }
}
node C(int x, int y)
{
    node ans;
    ans.clear();
    ans.l[1] = 1;
    for (int i = 1; i <= y; ++i)
    {
        mul(ans, x-i+1);
        div(ans, i);
    }
    return ans;
}
void geth(int x)
{
    if (flag[2][x]) return;
    flag[2][x] = 1;
    h[x] = pow(x*(x-1)/2);
}
void getg(int x)
{
    void getf(int x);
    if (flag[1][x]) return;
    flag[1][x] = 1;
    for (int i = 1; i < x; ++i)
    {
        getf(i);
        geth(x-i);
        g[x] = plus(g[x], mul(C(x-1, i-1), mul(f[i], h[x-i])));
    }
}
void getf(int x)
{
    if (flag[0][x]) return;
    flag[0][x] = 1;
    getg(x);
    geth(x);
    f[x] = minus(h[x], g[x]);
}
void deal()
{
    getf(n);
    f[n].print();
}
int main()
{
    while(scanf("%d", &n) != EOF)
    {
        if (!n) break;
        deal();
    }
    return 0;
}