HDU 4359 简单dp

题意：定义Bear Tree为一颗二叉树，这种二叉树每个结点有一个权值，范围在2^0~2^n-1，并且每个值只用一次，对于每个结点，如果同时存在左右子树，那么左子树的权值和要小于右子树的权值和。求点数为N，层次为D的Bear Tree的数量。

思路：

2^0 + 2^1 + ... + 2^n < 2^(n+1)

根据这个性质，我们可以得出权值最大节点必须在右子树上,因为假设权值最大的节点在左子树上，那么左边的权值就会大于右边的权值，不符合条件。

并且只要同时存在左右子树，则将权值最大节点放在右子树上就一定符合条件，因为此时右边的和一定会大于左边的和。

所以我们用dp[i][j]表示点数为i且深度不超过j的所有方案数，那么输出结果就是dp[n][d]-dp[n][d-1]。（深度不超过d的减去深度不超过d-1的即为深度为d的。）

而dp[n][d]的构成分下面两种：

1是只有左子树或只有右子树的情况，我们发现，只需要取任意一个节点来做根节点，乘以可能的子树情况（即dp[n-1][d-1]），再区别开是左子树还是右子树，总共有dp[n-1][d-1] (选出一个根节点后即为dp[n-1][d-1])* C(n,1)(从n个点中取一个做节点） * 2（左子树和右子树）种情况。

2是同时有左右子树的情况，我们假设左子树有k个节点，由于要满足左子树最大值<右子树最大值，所以选择的i-1个节点中最大的那个一定分配给右子树，所以左子树只能从i-2个节点中选k个做其节点则有dp[n-k-1][d-1]【右子树】 * dp[k][d-1]【左子树】 * C(n-2,k) 【左子树节点组成】* C(n,1)【根节点选择】

我们把一二相加即可得到转移方程，值得注意的是，由于n ,k<= 360随时可能爆精度，每次操作都尽可能模除10^9+7.

#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
#define maxn 369
__int64 maxmod=1000000000+7;
__int64 dp[maxn][maxn],c[maxn][maxn];//c表示组合数
int n,m,z,T;
int main()
{
    int i,j,k,z;
    memset(dp,0,sizeof(dp));//清空数组dp
    memset(c,0,sizeof(c)); //清空数组c
    
 for(i=1;i<=360;i++)
    {
         c[0][0]=1;
         c[i][0]=1;//C(0，0)=1;C(i,0)=1
    }
    
 for(i=1;i<=360;i++)
    for(j=1;j<=i;j++)
        c[i][j]=(c[i-1][j-1]+c[i-1][j])%maxmod;//组合公式打表 c（i,j)=c(i-1,j-1)+c(i-1,j)
    for(j=1;j<=360;j++)
    dp[1][j]=1;//只有一个节点情况只有一种
    
 for(i=2;i<=360;i++)
    for(j=1;j<=360;j++)
    {
    dp[i][j]=(c[i][i-1]*2*dp[i-1][j-1])%maxmod;  //c[i][i-1] 它其实也是在选出根结点，（只有一棵子树）
    for(k=1;k<=i-2;k++)//（有两棵子树）//一个节点为根节点，一个权值最大的节点分到右边，所以是i-2
    dp[i][j]=(dp[i][j]+((c[i][i-1]*c[i-2][k])%maxmod)*((dp[k][j-1]*dp[i-1-k][j-1])%maxmod))%maxmod;//一个子树加两个子树
    }  //(c[i][i-1]*c[i-2][k])是先选出根结点，把剩下的最大值分给右子树后再从i-2数中选出左子树的k个数
    scanf("%d",&T); //case数
    for(z=1;z<=T;z++)
    {
        scanf("%d%d",&n,&m);
        printf("Case #%d: %I64d\n",z,(dp[n][m]+maxmod-dp[n][m-1])%maxmod);
    }
    return 0;
}