经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C++ » 查看文章
BZOJ1812: [Ioi2005]riv(树形dp)
来源:cnblogs  作者:自为风月马前卒  时间:2018/12/7 9:33:22  对本文有异议

题意

题目链接

Sol

首先一个很显然的思路是直接用\(f[i][j] / g[i][j]\)表示\(i\)的子树中选了\(j\)个节点,该节点是否选的最小权值。但是直接这样然后按照树形背包的套路转移的话会有一种情况无法处理,就是说该节点不选,儿子节点也不选,这样我们就不清楚儿子节点的子节点的贡献了

一种暴力的做法是钦定该节点选,并重新枚举子树中的所有节点,转移出dp值之后背包合并

最后再把\(0\)号节点的合并一次

  1. #include<bits/stdc++.h>
  2. #define chmin(x, y) (x = x < y ? x : y)
  3. #define chmax(x, y) (x = x > y ? x : y)
  4. using namespace std;
  5. const int MAXN = 1001, INF = 2e9 + 10;
  6. inline int read() {
  7. char c = getchar(); int x = 0, f = 1;
  8. while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
  9. while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
  10. return x * f;
  11. }
  12. int N, K, a[MAXN], dis[MAXN], siz[MAXN], f[MAXN][MAXN], g[MAXN][MAXN], ans;
  13. vector<int> v[MAXN];
  14. void dfs2(int x, int fa, int root) {
  15. f[x][0] = dis[root] * a[x];
  16. for(int i = 0; i < v[x].size(); i++) {
  17. int to = v[x][i]; if(to == fa) continue;
  18. dfs2(to, x, root);
  19. for(int j = min(K, siz[x]); ~j; j--)
  20. for(int k = 0; k <= min(j, siz[to]); k++)
  21. chmax(f[x][j], f[to][k] + f[x][j - k]);
  22. }
  23. for(int i = siz[x]; i; i--) chmax(f[x][i], g[x][i]);
  24. }
  25. void dfs(int x, int fa) {
  26. dis[x] += dis[fa]; siz[x] = 1;
  27. for(int i = 0; i < v[x].size(); i++) {
  28. int to = v[x][i];
  29. if(to == fa) continue;
  30. dfs(to, x); siz[x] += siz[to];
  31. }
  32. g[x][0] = 0;
  33. memset(f, 0, sizeof(f));
  34. for(int i = 0; i < v[x].size(); i++) {
  35. int to = v[x][i]; if(to == fa) continue;
  36. dfs2(to, x, x);
  37. for(int j = min(K, siz[x]); ~j; j--)
  38. for(int k = 0; k <= min(j, siz[to]); k++)
  39. chmax(g[x][j], f[to][k] + g[x][j - k]);
  40. }
  41. for(int i = min(K, siz[x]); i; i--) g[x][i] = g[x][i - 1] + dis[x] * a[x];
  42. }
  43. int main() {
  44. N = read(); K = read();
  45. for(int i = 1; i <= N; i++) {
  46. a[i] = read(); int fa = read(); dis[i] = read();
  47. v[fa].push_back(i);
  48. }
  49. dfs(0, -1);
  50. for(int i = 1; i <= N; i++) ans += dis[i] * a[i];
  51. memset(f, 0, sizeof(f));
  52. int x = 0;
  53. for(int i = 0; i < v[x].size(); i++) {
  54. int to = v[x][i];
  55. dfs2(to, x, x);
  56. for(int j = min(K, siz[x]); ~j; j--)
  57. for(int k = 0; k <= min(j, siz[to]); k++)
  58. chmax(f[x][j], f[to][k] + f[x][j - k]);
  59. }
  60. int out = INF;
  61. for(int i = 0; i <= K; i++) out = min(out, ans - f[0][i]);
  62. printf("%d\n", out);
  63. return 0;
  64. }
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号