【BZOJ2286】【SDOI2011】消耗战

【题目链接】

【思路要点】

  • 补档博客,无题解。

【代码】

#include<bits/stdc++.h>
using namespace std;
#define MAXN	250005
#define INF	1e9
#define INFP	1e18
#define MAXLOG	20
struct edge {int dest, len;};
vector <edge> a[MAXN], b[MAXN];
int n, m, timer;
long long cost[MAXN];
bool mark[MAXN], valued[MAXN];
int dfn[MAXN], depth[MAXN], q[MAXN];
int father[MAXN][MAXLOG], minv[MAXN][MAXLOG];
bool cmp(int x, int y) {
	return dfn[x] < dfn[y];
}
void dfs(int pos, int fa, int len) {
	dfn[pos] = ++timer;
	depth[pos] = depth[fa] + 1;
	father[pos][0] = fa;
	minv[pos][0] = len;
	for (int i = 1; i < MAXLOG; i++) {
		father[pos][i] = father[father[pos][i - 1]][i - 1];
		minv[pos][i] = min(minv[pos][i - 1], minv[father[pos][i - 1]][i - 1]);
	}
	for (unsigned i = 0; i < a[pos].size(); i++)
		if (a[pos][i].dest != fa) dfs(a[pos][i].dest, pos, a[pos][i].len);
}
int Lca(int x, int y) {
	if (depth[x] < depth[y]) swap(x, y);
	for (int i = MAXLOG - 1; i >= 0; i--)
		if (depth[father[x][i]] >= depth[y]) x = father[x][i];
	if (x == y) return x;
	for (int i = MAXLOG - 1; i >= 0; i--)
		if (father[x][i] != father[y][i]) {
			x = father[x][i];
			y = father[y][i];
		}
	return father[x][0];
}
void insert(int x, int y) {
	int pos = y, minnum = INF;
	for (int i = MAXLOG - 1; i >= 0; i--)
		if (depth[father[pos][i]] >= depth[x]) {
			minnum = min(minnum, minv[pos][i]);
			pos = father[pos][i];
		}
	b[x].push_back((edge){y, minnum});
}
long long get(int pos, long long len) {
	if (mark[pos]) return len;
	long long ans = 0;
	for (unsigned i = 0; i < b[pos].size(); i++)
		ans += get(b[pos][i].dest, b[pos][i].len);
	return min(ans, len);
}
int main() {
	scanf("%d", &n);
	for (int i = 1; i < n; i++) {
		int x, y, value;
		scanf("%d%d%d", &x, &y, &value);
		a[x].push_back((edge){y, value});
		a[y].push_back((edge){x, value});
	}
	for (int j = 0; j < MAXLOG; j++)
		minv[0][j] = INF;
	dfs(1, 0, INF);
	int T;
	scanf("%d", &T);
	while (T--) {
		scanf("%d", &m);
		for (int i = 1; i <= m; i++) {
			scanf("%d", &q[i]);
			mark[q[i]] = true;
		}
		sort(q + 1, q + m + 1, cmp);
		static int Stack[MAXN], used[MAXN];
		int top = 1, cnt = 0;
		Stack[1] = 1;
		for (int i = 1; i <= m; i++) {
			int lca = Lca(q[i], Stack[top]);
			if (lca == Stack[top]) Stack[++top] = q[i];
			else {
				while (dfn[lca] < dfn[Stack[top - 1]]) {
					insert(Stack[top - 1], Stack[top]);
					used[++cnt] = Stack[top--];
				}
				if (dfn[lca] == dfn[Stack[top - 1]]) {
					insert(Stack[top - 1], Stack[top]);
					used[++cnt] = Stack[top--];
					Stack[++top] = q[i];
				} else {
					insert(lca, Stack[top]);
					used[++cnt] = Stack[top--];
					Stack[++top] = lca;
					Stack[++top] = q[i];
				}
			}
		}
		while (top >= 2) {
			insert(Stack[top - 1], Stack[top]);
			used[++cnt] = Stack[top--];
		}
		used[++cnt] = Stack[top--];
		for (int i = 1; i <= cnt; i++)
			cost[used[i]] = 0;
		printf("%lld\n", get(1, INFP));
		for (int i = 1; i <= cnt; i++)
			b[used[i]].clear();
		for (int i = 1; i <= m; i++)
			mark[q[i]] = false;
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值