HDU 5136 Yue Fei's Battle（计数DP）

最近不少题A了就不想写博客接着看别的题了。。这是懒么。。不过花点力气A的还是写写吧。

题意：直径为K的树问直径最大为n的树，有多少种，一节点最多只能有3条路就是说是一颗二叉树，每条分支都是等价的，也就是说如果翻转或旋转后一样，对称，都属于一种。

分奇数和偶数，如果是偶数，假设前n/2个节点加其他点的情况数为tmp，那么答案就是tmp×（tmp-1）/2+tmp，因为上半部分和下半部分如果不一样，换一下是一种情况，那么就会重复2次，上下一样就不会有重复情况。

到下一个奇数，中间点最多能加n/2个点，但是这样就会有3条主链，所以先不考虑这种，中间点加小于等于n/2-1个点的情况个数是那个tmp的前缀和sum[n/2-1]，所以这种情况答案就是sum[k/2-1]*(tmp*(tmp-1)/2+tmp)，如果中间加了n/2个点，那么要分3种情况讨论。

1.3个主链的样子都不一样。那么就是tmp*(tmp-1)*(tmp-2)重复情况有3！种也就是6种。

2.2个主链样子一样第三个不一样，那么就是tmp*(tmp-1)，没有重复。

3.3个主链都一样，那么就是tmp。

所以最后计算式是sum[n/2-1]*(tmp*(tmp-1)/2+tmp)+tmp*(tmp-1)*(tmp-2)/6+tmp*(tmp-1)+tmp。

最后考虑一下tmp该如何递推，容易看出只需要在偶数的时候递推就可以了，因为下一个奇数还是上一个偶数的tmp。式子跟上一个推导很像tmp = sum[n/2-2]*tmp+tmp*(tmp-1)/2+tmp;意思是在第n/2个节点（之前的已经递推过了）分2种，1种是放小于n/2-1个，另一种是放n/2-1个，那么就是2条主链了。

AC代码：

//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<ctype.h>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
#include<cstdlib>
#include<stack>
#include<cmath>
#include<queue>
#include<set>
#include<map>
#include<ctime>
#include<string.h>
#include<string>
#include<sstream>
#include<bitset>
using namespace std;
#define ll __int64
#define ull unsigned long long
#define eps 1e-8
#define NMAX 1000000005
#define MOD 1000000007
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define PI acos(-1)
template<class T>
inline void scan_d(T &ret)
{
    char c;
    int flag = 0;
    ret=0;
    while(((c=getchar())<'0'||c>'9')&&c!='-');
    if(c == '-')
    {
        flag = 1;
        c = getchar();
    }
    while(c>='0'&&c<='9') ret=ret*10+(c-'0'),c=getchar();
    if(flag) ret = -ret;
}
ll dp[100005],sum[100005];
void cal()
{
    ll a = 3,b = MOD-2,ret = 1;
    while(b)
    {
        if(b&1) ret = ret*a%MOD;
        a = a*a%MOD;
        b >>= 1;
    }
    cout<<ret<<endl;
}
int main()
{
#ifdef GLQ
    freopen("input.txt","r",stdin);
//    freopen("o1.txt","w",stdout);
#endif // GLQ
//    cal();
    dp[1] = dp[2] = 1;
    sum[0] = 1;
    sum[1] = 2;
    ll tmp = 1;
    for(int i = 3; i <= 100000; i++)
    {
        if(i%2)
        {
            ll a = ((tmp*(tmp-1LL+MOD)%MOD*(tmp-2LL+MOD)%MOD)*166666668LL%MOD + tmp*(tmp-1LL+MOD)%MOD + tmp)%MOD;
            dp[i] = (sum[i/2-1]*((tmp*(tmp-1LL+MOD)%MOD)*500000004LL%MOD+tmp)%MOD+a)%MOD;
        }
        else
        {
            ll a = tmp*sum[i/2-2]%MOD;
            tmp = (tmp*(tmp-1LL+MOD)%MOD*500000004LL%MOD+tmp+a)%MOD;
            dp[i] = (tmp*(tmp-1LL+MOD)%MOD*500000004LL%MOD+tmp)%MOD;
            sum[i/2] = (sum[i/2-1]+tmp)%MOD;
        }
    }
    int n;
    while(~scanf("%d",&n) && n)
    {
        printf("%I64d\n",dp[n]);
    }
    return 0;
}