https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=5015

设dp[cur][i][j]表示当前是第cur个顶点,自身状态是i(0或者1),爸爸的状态是j(0或者1)的时候的最多白色节点数。

0---白色

1---黑色

如果不合法,就是第cur号顶点不能染成白色,则为-1,注意到任意节点都可以染成黑色,所以dp[cur][1][0] = dp[cur][1][1] = 0

转移:

首先把叶子节点特判掉,因为叶子节点的状态很容易判断,并且只有k == 1的时候,才有dp[left][0][1] = 1

然后留给每一个爸爸的转移就是:

1、如果爸爸染黑色,则从  dp[son][1][1]和dp[son][0][1]娶个max过来即可。

2、比较麻烦的是爸爸染了白色。这样就相当于对于所有的儿子(儿子的状态已经全部算出来了),把k个染成黑色(使得爸爸合法),剩下的染成白色,使得白色节点数最大,

也就是给出一个结构体数组,有arr[i].a和arr[i].b表示这个节点染成白色,得到arr[i].a贡献,这个节点染成黑色,得到arr[i].b贡献。

设d[n][k]表示前n个物品,确定选了k个做黑色的最大贡献。二维费用背包转移即可。hack: 需要用到d[i][0]

所以d[0][0] = 0,而d[1][0] = arr[1].a    d[2][0] = arr[1].a + arr[2].a

这题写那个二维费用dp的时候坑队友了,没写出来,转移的时候没考虑d[i][0]

#include <bits/stdc++.h>
#include <algorithm>
#define inf (0x3f3f3f3f)
using namespace std;
typedef long long int LL;
const int maxn = 1e2 + ;
int dp[maxn][][];
char str[ + ];
struct Edge {
int u, v, tonext;
} e[maxn * ];
int first[maxn], num;
void addEdge(int u, int v) {
e[num].u = u, e[num].v = v, e[num].tonext = first[u];
first[u] = num++;
}
int son[maxn];
bool in[maxn];
int n, k; void show() {
for (int i = ; i <= n; ++i) {
printf("node %d: ", i);
for (int j = first[i]; ~j; j = e[j].tonext) {
printf("%d ", e[j].v);
}
printf("\n");
}
printf("***************\n");
}
vector<int> vc[maxn];
int d[maxn][];
struct Node {
int a, b;
Node(int _a, int _b) {
a = _a, b = _b;
}
Node() {}
} arr[maxn];
int getMax(struct Node arr[], int n, int k) {
if (k < ) return -inf;
if (k == ) {
int sum = ;
for (int i = ; i <= n; ++i) sum += arr[i].a;
return sum;
}
if (n < k) return -inf;
memset(d, -0x3f, sizeof d);
d[][] = ;
for (int i = ; i <= n; ++i) {
for (int j = k; j >= ; --j) {
if (d[i - ][j] >= ) d[i][j] = max(d[i][j], d[i - ][j] + arr[i].a);
if (j >= && d[i - ][j - ] >= ) d[i][j] = max(d[i][j], d[i - ][j - ] + arr[i].b);
}
}
return d[n][k];
}
void dfs(int cur) {
if (!son[cur]) return;
for (int i = first[cur]; ~i; i = e[i].tonext) {
int v = e[i].v;
vc[cur].push_back(v);
dfs(v);
dp[cur][][] += max(dp[v][][], dp[v][][]);
dp[cur][][] += max(dp[v][][], dp[v][][]); // 爸爸是黑色
}
int sel = , val = , to = ;
for (int i = ; i < vc[cur].size(); ++i) {
int v = vc[cur][i];
if (dp[v][][] >= ) {
arr[++to] = Node(dp[v][][], dp[v][][]);
} else {
sel++; //这些不能变成白色,也就是固定必须是黑色
val += dp[v][][]; //其贡献
}
}
dp[cur][][] = getMax(arr, to, k - - sel) + val;
dp[cur][][] = getMax(arr, to, k - sel) + val;
if (dp[cur][][] < ) dp[cur][][] = -;
else dp[cur][][]++;
if (dp[cur][][] < ) dp[cur][][] = -;
else dp[cur][][]++;
}
void work() {
num = ;
memset(in ,false, sizeof in);
memset(first, -, sizeof first);
memset(dp, -, sizeof dp);
memset(son, false, sizeof son);
scanf("%d%d", &n, &k);
for (int i = ; i <= maxn - ; ++i) vc[i].clear();
getchar();
for (int i = ; i <= n; ++i) {
gets(str + );
// printf("%s\n", str + 1);
int lenstr = strlen(str + );
for (int j = ; j <= lenstr;) {
if (str[j] >= '' && str[j] <= '') {
int fuck = str[j] - '';
++j;
while (j <= lenstr && str[j] >= '' && str[j] <= '') {
fuck = fuck * + str[j] - '';
++j;
}
if (fuck == ) break;
son[i]++;
addEdge(i, fuck);
in[fuck] = true;
} else j++;
}
}
int root = ;
for (int i = ; i <= n; ++i) {
if (!in[i]) {
root = i;
break;
}
}
// show();
if (n == ) {
if (k != ) {
printf("0\n");
return;
}
}
for (int i = ; i <= n; ++i) {
dp[i][][] = dp[i][][] = ;
}
if (k == ) {
printf("%d\n", n);
return;
}
if (k == ) {
for (int i = ; i <= n; ++i) {
if (!son[i]) {
dp[i][][] = ;
dp[i][][] = -;
dp[i][][] = ;
dp[i][][] = ;
}
}
}
dfs(root);
int ans = dp[root][][];
ans = max(ans, dp[root][][]);
ans = max(ans, dp[root][][]);
// ans = max(ans, dp[root][0][1]);
printf("%d\n", ans);
} int main() {
#ifdef local
freopen("data.txt", "r", stdin);
#endif
int t;
scanf("%d", &t);
while (t--) work();
return ;
}

Regionals 2014 >> Asia - Taichung 7003 - A Balance Game on Trees 树形DP + 二维费用背包的更多相关文章

  1. 2014 Super Training #9 E Destroy --树的直径+树形DP

    原题: ZOJ 3684 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3684 题意: 给你一棵树,树的根是树的中心(到其 ...

  2. [The Preliminary Contest for ICPC Asia Nanjing 2019] A-The beautiful values of the palace(二维偏序+思维)

    >传送门< 前言 这题比赛的时候觉得能做,硬是怼了一个半小时,最后还是放弃了.开始想到用二维前缀和,结果$n\leq 10^{6}$时间和空间上都爆了,没有办法.赛后看题解用树状数组,一看 ...

  3. The Preliminary Contest for ICPC Asia Xuzhou 2019 J Random Access Iterator (树形DP)

    每次循环向下寻找孩子时,随机选取一个孩子,设dp[u]为从u出发,不能得出正确答案的概率,则从u出发,走一次的情况下不能得出正确答案的概率是 P = (dp[v1]+dp[v2]+dp[v3]+--d ...

  4. The Preliminary Contest for ICPC Asia Shenyang 2019 D. Fish eating fruit(树形dp)

    题意:求一棵树上所有路径和模3分别为0 1 2 的权值的和 思路:树形dp 增加一个记录儿子节点满足条件的个数的数组 不要放在一起dp不然答案跟新会有问题 #include <bits/stdc ...

  5. 2015 UESTC Winter Training #4【Regionals 2008 :: Asia - Tehran】

    2015 UESTC Winter Training #4 Regionals 2008 :: Asia - Tehran 比赛开始时电脑死活也连不上WIFI,导致花了近1个小时才解决_(:зゝ∠)_ ...

  6. C++ 洛谷 2014 选课 from_树形DP

    洛谷 2014 选课 没学树形DP的,看一下. 首先要学会多叉树转二叉树. 树有很多种,二叉树是一种人人喜欢的数据结构,简单而且规则.但一般来说,树形动规的题目很少出现二叉树,因此将多叉树转成二叉树就 ...

  7. hdu5071 2014 Asia AnShan Regional Contest B Chat

    模拟题: add的时候出现过的则不再添加 close的时候会影响到top rotate(Prior.Choose)的时候会影响到top /*============================== ...

  8. UVALive 7138 The Matrix Revolutions(Matrix-Tree + 高斯消元)(2014 Asia Shanghai Regional Contest)

    题目链接:https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&category=6 ...

  9. UVALive 7143 Room Assignment(组合数学+DP)(2014 Asia Shanghai Regional Contest)

    题目链接:https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&category=6 ...

随机推荐

  1. bzoj 2850 巧克力王国 —— K-D树

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2850 只要暴力判断是否全选一个子树或全不选,如果都不是就进入查询: 要注意值有负,所以不是直 ...

  2. bzoj 3994 约数个数和 —— 反演+数论分块

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3994 推导过程和这里一样:https://www.cnblogs.com/MashiroSk ...

  3. git常见使用情境整理

    一.版本回退 回退到某个commit版本的方法如下: 1. 查看commit历史 git log 找到想要回退的版本的号码,eg:f765889 2. 回退到该版本 git reset f765889 ...

  4. 在CentOS 7上安装Node.js的4种方法(包含npm)

    Node.js和Javascript有着千丝万缕的联系,可以说Node.js让Javascript显得从未如此强大.好吧…微魔其实是个门外汉…但是这并不能阻碍微魔学习探索未知的信心~今天在国外闲逛,看 ...

  5. 无密码登录Linux服务器

    1.使用windows上的SecureCRT生成密钥对. Tools->Create Public Key..->RSA->Passphrase(最好输入,也可为空)->长度默 ...

  6. puppet插件fact和hiera(puppet自动化系列3)

    四.Fact插件 4.1 使用pluginsync进行发布 这种方法比较特殊,节点factpath目录里除了编写好的rb文件之外,还需要在puppet模块中引用,运行一次之后才会转换成fact.通常在 ...

  7. Python:模块详解及import本质

    转于:http://www.cnblogs.com/itfat/p/7481972.html 博主:东大网管 一.定义: 模块:用来从逻辑上组织python代码(变量,函数,类,逻辑:实现一个功能), ...

  8. netty中的EventLoop和EventLoopGroup

    Netty框架的主要线程就是I/O线程,线程模型设计的好坏,决定了系统的吞吐量.并发性和安全性等架构质量属性. 一.Netty的线程模型 在讨论Netty线程模型时候,一般首先会想到的是经典的Reac ...

  9. 查看,检查,修复pg的命令

    标签(空格分隔): ceph,ceph运维,pg 如果集群状态是HEALTH_ERR 并且有pgs inconsistent,需要进行如下操作: 1. 通过下面的命令查看哪些pg状态不一致: # ce ...

  10. UDK编辑器 49条小提示

    转自:http://www.cnblogs.com/hmxp8/archive/2012/02/09/2343674.html Very Helpful~ 01. First time using t ...