博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[HDOJ1561]The more, The Better(树dp,背包)
阅读量:5149 次
发布时间:2019-06-13

本文共 2856 字,大约阅读时间需要 9 分钟。

题目链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=1561

 

思路:树dp,设dp(i,j)为i节点时要攻破k个城堡可获得的最大收益(包括i节点,即i节点也要消耗k个城堡中的一个),我们考虑攻破的这k个城堡应当包括哪些部分:

假如节点i有一个儿子,则这k个城堡一定来自它的儿子已经它儿子的儿子;假如它有p个儿子,则可能来自于p个儿子中的q个儿子,这样子问题就找到了。

设当前节点为u,攻破j个城堡的最大收益为dp(u,j),其中有k个城堡是从u的直接儿子获得的,则有j-k个是从儿子的儿子获得的,则就变成了熟知的01背包,表达为:

dp[u][j] = max(dp[u][j], dp[u][k]+dp[v][j-k]),注意如果u是虚根的话是不可以算入花费的。所以当u=0的时候我们相当于求从k+1个城堡中获得收益。

1 /*  2 ━━━━━┒ギリギリ♂ eye!  3 ┓┏┓┏┓┃キリキリ♂ mind!  4 ┛┗┛┗┛┃\○/  5 ┓┏┓┏┓┃ /  6 ┛┗┛┗┛┃ノ)  7 ┓┏┓┏┓┃  8 ┛┗┛┗┛┃  9 ┓┏┓┏┓┃ 10 ┛┗┛┗┛┃ 11 ┓┏┓┏┓┃ 12 ┛┗┛┗┛┃ 13 ┓┏┓┏┓┃ 14 ┃┃┃┃┃┃ 15 ┻┻┻┻┻┻ 16 */ 17 #include 
18 #include
19 #include
20 #include
21 #include
22 #include
23 #include
24 #include
25 #include
26 #include
27 #include
28 #include
29 #include
30 #include
31 #include
32 #include
33 #include
34 #include
35 using namespace std; 36 #define fr first 37 #define sc second 38 #define cl clear 39 #define BUG puts("here!!!") 40 #define W(a) while(a--) 41 #define pb(a) push_back(a) 42 #define Rint(a) scanf("%d", &a) 43 #define Rs(a) scanf("%s", a) 44 #define Cin(a) cin >> a 45 #define FRead() freopen("in", "r", stdin) 46 #define FWrite() freopen("out", "w", stdout) 47 #define Rep(i, len) for(int i = 0; i < (len); i++) 48 #define For(i, a, len) for(int i = (a); i < (len); i++) 49 #define Cls(a) memset((a), 0, sizeof(a)) 50 #define Clr(a, x) memset((a), (x), sizeof(a)) 51 #define Full(a) memset((a), 0x7f7f7f, sizeof(a)) 52 #define lrt rt << 1 53 #define rrt rt << 1 | 1 54 #define pi 3.14159265359 55 #define RT return 56 #define lowbit(x) x & (-x) 57 #define onecnt(x) __builtin_popcount(x) 58 typedef long long LL; 59 typedef long double LD; 60 typedef unsigned long long ULL; 61 typedef pair
pii; 62 typedef pair
psi; 63 typedef pair
pll; 64 typedef map
msi; 65 typedef vector
vi; 66 typedef vector
vl; 67 typedef vector
vvl; 68 typedef vector
vb; 69 70 const int maxn = 220; 71 vi G[maxn]; 72 int n, m; 73 int w[maxn]; 74 bool vis[maxn]; 75 int dp[maxn][maxn]; 76 77 void dfs(int u) { 78 For(i, 1, m+1) dp[u][i] = w[u]; 79 vis[u] = 1; 80 Rep(i, G[u].size()) { 81 int v = G[u][i]; 82 if(vis[v]) continue; 83 dfs(v); 84 if(u == 0) { 85 for(int j = m; j >= 0; j--) { 86 for(int k = j; k >= 0; k--) { 87 dp[u][j] = max(dp[u][j], dp[u][k]+dp[v][j-k]); 88 } 89 } 90 } 91 else { 92 for(int j = m; j >= 0; j--) { 93 for(int k = j - 1; k >= 1; k--) { 94 dp[u][j] = max(dp[u][j], dp[u][k]+dp[v][j-k]); 95 } 96 } 97 } 98 } 99 }100 101 int main() {102 // FRead();103 int a, b;104 while(~Rint(n) && ~Rint(m) && n + m) {105 Rep(i, n+5) G[i].clear();106 Cls(dp); Cls(vis);107 For(i, 1, n+1) {108 Rint(a); Rint(b);109 G[a].push_back(i);110 w[i] = b;111 }112 dfs(0);113 printf("%d\n", dp[0][m]);114 }115 RT 0;116 }

 

转载于:https://www.cnblogs.com/kirai/p/5815521.html

你可能感兴趣的文章
Opencv保存摄像头视频&&各种编码器下视频文件占用空间对比
查看>>
「图形学」直线扫描——Bresenham算法改进了中点Bresenham算法?
查看>>
jQuery 给div绑定单击事件
查看>>
Exceptionless 生产部署笔记
查看>>
有关快速幂取模
查看>>
转 ObjExporter Unity3d导出场景地图寻路
查看>>
Linux运维必备工具
查看>>
Ubuntu配置ssh及vnc
查看>>
C语言进阶——const 和 volatile 分析09
查看>>
字符串的查找删除
查看>>
NOI2018垫底记
查看>>
快速切题 poj 1002 487-3279 按规则处理 模拟 难度:0
查看>>
Codeforces Round #277 (Div. 2)
查看>>
一步步学Mybatis-搭建最简单的开发环境-开篇(1)
查看>>
微信小程序图片上传
查看>>
【更新】智能手机批量添加联系人
查看>>
NYOJ-128前缀式计算
查看>>
centos6.7 配置外网端口映射
查看>>
淡定,啊。数据唯一性
查看>>
深入理解 JavaScript 事件循环(一)— event loop
查看>>