HDU 1561 The more, The Better

Problem Description

ACboy很喜欢玩一种战略游戏，在一个地图上，有N座城堡，每座城堡都有一定的宝物，在每次游戏中ACboy允许攻克M个城堡并获得里面的宝物。但由于地理位置原因，有些城堡不能直接攻克，要攻克这些城堡必须先攻克其他某一个特定的城堡。你能帮ACboy算出要获得尽量多的宝物应该攻克哪M个城堡吗？

Input

每个测试实例首先包括2个整数，N,M.(1 <= M <= N <= 200);在接下来的N行里，每行包括2个整数，a,b. 在第 i 行，a 代表要攻克第 i 个城堡必须先攻克第 a 个城堡，如果 a = 0 则代表可以直接攻克第 i 个城堡。b 代表第 i 个城堡的宝物数量, b >= 0。当N = 0, M = 0输入结束。

Output

对于每个测试实例，输出一个整数，代表ACboy攻克M个城堡所获得的最多宝物的数量。

Sample Input

3 2
0 1
0 2
0 3
7 4
2 2
0 1
0 4
2 1
7 1
7 6
2 2
0 0

Sample Output

5
13

解题思路：经典的树形ＤＰ入门题。要取子节点，必须先取父节点。同样为ＤＰ求最优解问题。

＄ｄｐ［ｉ］［ｊ］＄：父节点为ｉ，选取其及其下节点数为ｊ的最优解（即最大值）。

状态转移方程：＄ｄｐ［ｉ］［ｊ］　＝　＼ｍａｘ＼｛ｄｐ［ｉ］［ｊ］，ｄｐ［ｉ］［ｊ－ｋ］＋ｄｐ［ｓｏｎ［ｉ］］［ｋ］＄

状态初始化：＄ｄｐ［ｉ］［１］　＝　ｖａｌ［ｉ］＄；（即以ｉ为父节点，只选取１个节点时，其值为选取ｉ点的值）

注意几个细节问题，在代码中均有注释。

#pragma comment(linker, "/STACK:1024000000,1024000000") //黑科技
#include <stdio.h>
#include <map>
#include <vector>
#include <queue>
#include <iostream>
#include <string.h>
#include <algorithm>
#include <string>
#include <math.h>
using namespace std;
typedef long long LL;
const double eqs=1e-9;
//freopen("data.in", "r", stdin);
#define maxd 210
int n, m;
vector <int > G[maxd]; 
int dp[maxd][maxd];
void dfs(int u)
{
    int len = G[u].size();
    for(int i = 0; i < len; i++) dfs(G[u][i]);//预处理所有子节点的值
    for(int i = 0; i < len; i++)
        for(int j = m + 1; j >= 0; j--) //逆推，类似于01背包，所有子节点只能更新一次。 
            for(int k = 0; k <= j - 1; k++) //k + 1 <= j， 保证父节点被取 
                dp[u][j] = max(dp[u][j], dp[u][j - k] + dp[G[u][i]][k]);
}
int main()
{
    #ifndef ONLINE_JUDGE
    //freopen("data.in", "r", stdin);
    #endif
    while(scanf("%d %d", &n, &m) != EOF && n)
    {
        for(int i = 0; i < maxd; i++) G[i].clear();//init
        memset(dp, 0, sizeof(dp));//标记 未求值
        for(int i = 1; i <= n; i++)
        {
            int a, b;
            scanf("%d %d", &a, &b);
            G[a].push_back(i);
            dp[i][1] = b;
        }
        dfs(0);
        printf("%d\n", dp[0][m + 1]);
    }
    return 0;
}