显然考虑二分,设当前二分出的值是 \(lim\)

那么就是称一个点能覆盖另一个点当且仅当它被选中且它与那个点的距离不大于 \(lim\),要判断是否能选出不多于 \(m\) 个点覆盖到所有的关键节点。

\(f_i\) 表示节点 \(i\) 的子树中距离 \(i\) 最近的被选中的节点到 \(i\) 的距离,\(g_i\) 表示节点 \(i\) 的子树中距离 \(i\) 最远的还没有被覆盖的关键节点\(i\) 的距离。

定义有点绕,可以在这里停下来读一读。

假设当前在点 \(u\)

初始时 \(f_u=\infty,g_u=-\infty\)

然后 \(f_u=\min\limits_{v\in son_u}(f_u,f_v+1)\)\(g_u=\max\limits_{v\in son_u}(g_u,g_v+1)\)

然后如果点 \(u\) 是关键节点并且 \(f_u\gt lim\),那么 \(g_u=\max(g_u,0)\)

如果 \(f_u+g_u\le lim\),那么 \(g_u=-\infty\)

如果 \(g_u=lim\),那么 \(u\) 必须被选择,所以 \(f_u=0,g_u=-\infty\)

最后如果根节点的 \(g\ge 0\),那么根节点也必须被选择。

Code:

#include <bits/stdc++.h>
using namespace std;
const int N = 300005, inf = 0x3f3f3f3f;
int n, m;
int ty[N];
int head[N], ver[N*2], nxt[N*2], cnt;
int f[N], g[N], tot = 0;

void add(int u, int v) {
	ver[++cnt] = v, nxt[cnt] = head[u], head[u] = cnt;
}

void dfs(int u, int fa, int lim) {
	f[u] = inf, g[u] = -inf;
	for (int i = head[u]; i; i = nxt[i]) {
		int v = ver[i];
		if (v == fa) continue;
		dfs(v, u, lim);
		f[u] = min(f[u], f[v] + 1);
		g[u] = max(g[u], g[v] + 1);
	}
	if (ty[u] && f[u] > lim) g[u] = max(g[u], 0);
	if (g[u] + f[u] <= lim) g[u] = -inf;
	if (g[u] == lim) g[u] = -inf, f[u] = 0, ++tot;
}

bool check(int x) {
	tot = 0, dfs(1, 0, x);
	if (g[1] >= 0) ++tot;
	return tot <= m;
}

int main() {
	scanf("%d%d", &n, &m);
	for (int i = 1; i <= n; ++i) scanf("%d", &ty[i]);
	for (int i = 1, u, v; i < n; ++i) scanf("%d%d", &u, &v), add(u, v), add(v, u);
	int l = 0, r = n - 1;
	while (l < r) {
		int mid = l + r >> 1;
		if (check(mid)) r = mid;
		else l = mid + 1;
	}
	printf("%d", l);
	return 0;
}

原文地址:http://www.cnblogs.com/Kobe303/p/16885931.html

1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长! 2. 分享目的仅供大家学习和交流,请务用于商业用途! 3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入! 4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解! 5. 如有链接无法下载、失效或广告,请联系管理员处理! 6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需! 7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员! 8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载 声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性