BZOJ [FJOI2007]轮状病毒

分析：

我其实想到链状怎么做了，没有往下多想。一直想用组合做。

最后看到CLJ的题解，才发现我已经想对了一半。。。

对于一个链来说，设Dp[i]为长度为i的链和一个“中心”的生成树数量，可以Dp之。。
然后考虑环，枚举1节点所在的环的长度，设为Len，那么这个环的位置有Len种情况，这个环与中心连接也有Len种情况，那么这个节点对答案的贡献就是Len^2*Dp[n-Len]。。

转自：http://hi.baidu.com/wjbzbmr/item/3f5e90f02b913ed66225d201

#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <iostream>

#define N 50

using namespace std;

struct BIGN
{
    int a[N];
}dp[110],ans;

int n;

inline BIGN operator +(BIGN a,BIGN b)
{
    BIGN c; memset(c.a,0,sizeof c.a);
    c.a[0]=max(a.a[0],b.a[0]);
    int jin=0;
    for(int i=1;i<=c.a[0];i++)
    {
        jin+=a.a[i]+b.a[i];
        c.a[i]=jin%10000;
        jin/=10000;
    }
    while(jin) c.a[++c.a[0]]=jin%10000,jin/=10000;
    return c;
}

inline BIGN operator *(BIGN a,int b)
{
    BIGN c; memset(c.a,0,sizeof c.a);
    int jin=0;
    for(int i=1;i<=a.a[0];i++)
    {
        jin+=c.a[i]+a.a[i]*b;
        c.a[i]=jin%10000;
        jin/=10000;
    }
    c.a[0]=a.a[0];
    while(jin) c.a[++c.a[0]]=jin%10000,jin/=10000;
    return c;
}

void prev()
{
    dp[0].a[0]=dp[0].a[1]=1;
    dp[1].a[1]=dp[1].a[0]=1;
    for(int i=2;i<=100;i++)
    {
        dp[i].a[1]=0;dp[i].a[0]=1;
        for(int j=1;j<=i;j++)
            dp[i]=dp[i]+dp[i-j]*j;
    }
    
}

void prt(const BIGN &a)
{
    printf("%d",a.a[a.a[0]]);
    for(int i=a.a[0]-1;i>=1;i--) printf("%.4d",a.a[i]);
    puts("");
}

int main()
{
    prev();
    scanf("%d",&n);
    ans.a[1]=0; ans.a[0]=1;
    for(int i=1;i<=n;i++)
        ans=ans+dp[n-i]*(i*i);
    prt(ans);
    return 0;
}