tyvj 2055 [Nescafé29]圣主的考验

http://www.elijahqi.win/archives/1400
背景
随着四叶草魔杖的插入，Nescafe之塔的塔尖附近出现了光漩涡。随着一束光芒投下，探险队员们面前出现了圣主的身影。
“参见圣主！”
“大家好，我叫杜宇飞。我就是Nescafe的圣主applepi。”
“圣主，我们是专程前来拜访您的！”探险队员们万分激动。
“再有七天就是Nescafe镇塔之宝一年一度的封印之期了，你们可以留下来观看封印。不过在这之前，我还是想考考你们……”

描述
若对于二叉树T的每个节点v，其左子树的高度L和右子树的高度R均满足|L – R|≤1，则这个树T有可能来自超自然之界。规定若某节点子树为空，则该子树的高度是0。你的任务是求有N个节点的可能来自超自然之界的树的数目。

输入格式
每个测试点包含若干个测试数据。
每个测试数据占一行，包含一个整数N。
输入文件以0结尾。

输出格式
对于每个测试数据，在单独的一行内输出结果。由于结果可能会很大，你只需要输出答案在十进制表示下的后九位。若答案不足九位，只需输出原答案。

测试样例1
输入
2
3
5
30
0
输出
2
1
6
11307920
备注
对于30% 的测试点，N≤100。
对于70% 的测试点，N≤1000。
对于100% 的测试点，1≤N≤3000。

这题很树形dp 简单的dp方程可以认为是dp[i][j]表示现在有i个节点深度为j的方案数是多少

时间复杂度o(n^3) 显然会炸 怎么办

题解给了个更好的枚举方法 我们每次取枚举我的深度 然后同时可以预处理 若深度为这么大 我最多能放多少节点 反之最少需要多少子节点

我每次枚举i表示 这个深度的结点个数是多少 然后根据朴素dp可以知道 一共只有三种转移可能

左边右边同时j-1 或者左边j-1右边j-2 左边j-2右边j-1

#include<cstdio>
#include<algorithm>
#define N 3300
#define mod 1000000000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while (ch<'0'||ch>'9') ch=gc();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x;
}
long long dp[N][N];
int max1[N],min1[N],ans[N];
int main(){
    //freopen("fin.in","r",stdin);
    min1[1]=1;max1[1]=1;dp[0][0]=dp[1][1]=1;ans[1]=1;
    for (int dep=2;dep;dep++){
        min1[dep]=min1[dep-1]+min1[dep-2]+1;
        max1[dep]=(max1[dep-1]<<1)+1;if (min1[dep]>3000) break;
        for (int i=min1[dep];i<=max1[dep]&&i<=3000;++i){
            for (int j=min1[dep-1];j<=max1[dep-1]&&j<i;++j) (dp[i][dep]+=dp[j][dep-1]*dp[i-j-1][dep-1])%=mod;
            for (int j=min1[dep-2];j<=max1[dep-2]&&j<i;++j){(dp[i][dep]+=dp[j][dep-2]*dp[i-j-1][dep-1]*2%mod)%=mod;}
            (ans[i]+=dp[i][dep])%=mod;
        }
    }
    while (1){
        int n=read();if (!n) break;
        if (n>=38) printf("%09d\n",ans[n]);else printf("%d\n",ans[n]);
    }
    return 0;
}