[SDOI2011]消耗战

题目链接

题目描述

在一场战争中，战场由 $n$ 个岛屿和 $n-1$ 个桥梁组成，保证每两个岛屿间有且仅有一条路径可达。现在，我军已经侦查到敌军的总部在编号为 $1$ 的岛屿，而且他们已经没有足够多的能源维系战斗，我军胜利在望。已知在其他 $k$ 个岛屿上有丰富能源，为了防止敌军获取能源，我军的任务是炸毁一些桥梁，使得敌军不能到达任何能源丰富的岛屿。由于不同桥梁的材质和结构不同，所以炸毁不同的桥梁有不同的代价，我军希望在满足目标的同时使得总代价最小。

侦查部门还发现，敌军有一台神秘机器。即使我军切断所有能源之后，他们也可以用那台机器。机器产生的效果不仅仅会修复所有我军炸毁的桥梁，而且会重新随机资源分布（但可以保证的是，资源不会分布到 $1$ 号岛屿上）。不过侦查部门还发现了这台机器只能够使用 $m$ 次，所以我们只需要把每次任务完成即可。

输入格式

第一行一个整数 $n$，表示岛屿数量。

接下来 $n-1$ 行，每行三个整数 $u,v,w$，表示 $u$ 号岛屿和 $v$ 号岛屿由一条代价为 $w$ 的桥梁直接相连。

第 $n+1$ 行，一个整数 $m$，代表敌方机器能使用的次数。

接下来 $m$ 行，第 $i$ 行一个整数 $k_i$，代表第 $i$ 次后，有 $k_i$ 个岛屿资源丰富。接下来 $k_i$ 个整数 $h_1,h_2,...,h_{k_i}$，表示资源丰富岛屿的编号。

输出格式

输出共 $m$ 行，表示每次任务的最小代价。

输入输出样例

输入 #1

10
1 5 13
1 9 6
2 1 19
2 4 8
2 3 91
5 6 8
7 5 4
7 8 31
10 7 9
3
2 10 6
4 5 7 8 3
3 9 4 6

输出 #1

12
32
22

说明/提示

对于 $100\%$ 的数据，$2\le n\le 2.5\times 10^5,1\le m\le 5\times 10^5,\sum k_i\le 5\times 10^5,1\le k_i<n,h_i\ne 1,1\le u,v\le n,1\le w\le 10^5。$

---

Solution

• 考虑到 $m$ 最大可达 $10^5$，但 $\sum _{i=1}^m{k}_i$ 的最大值却不超过 $5\times 10^5$，想到用 虚树 $+$ 树形DP 实现。

• 不难得出，状态转移方程为
  f [ u ] = { w [ u ] , u ∈ { h i } min ⁡ { w [ u ] , ∑ f a [ v ] = u f [ v ] } , u ∉ { h i } f[u]= \begin{cases} w[u], & u\in \{h_i\}\\ \min\{w[u],\sum\limits_{fa[v]=u}f[v]\},& u\notin \{h_i\}\\ \end{cases} f[u]=⎩⎨⎧​w[u],min{w[u],fa[v]=u∑​f[v]},​u∈{hi​}u∈/​{hi​}​

• 其中 $w[u]$ 表示 $u$ 到虚树中 $fa[u]$ 路径中代价最小值，可以通过倍增求解。

• 特别地，$f[1]=\sum _{fa[v]=1}f[v]$，$f[1]$ 即为答案。

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int maxn=500010,maxl=20;
const LL INF=0x3f3f3f3f3f3f3f3f;
int Fir[maxn],Nxt[maxn],Too[maxn],tot=1;
int n,m,dfn[maxn],Times,h[maxn],dep[maxn];
int fa[maxn][maxl];
int stk[maxn],top;
LL f[maxn],w[maxn],lmin[maxn][maxl];
bool ask[maxn];
template<class T>inline T Min(T a,T b){return a<b?a:b;}
inline void add(int a,int b,LL c){
	Too[++tot]=b;w[tot]=c;Nxt[tot]=Fir[a];Fir[a]=tot;
}
inline bool cmp(int a,int b){
	return dfn[a]<dfn[b];
}
inline void dfs(int u){
	dfn[u]=++Times;
	for(int i=1;i<=19;++i){
		fa[u][i]=fa[fa[u][i-1]][i-1];
		lmin[u][i]=Min(lmin[u][i-1],lmin[fa[u][i-1]][i-1]);
	}
	for(int i=Fir[u];i;i=Nxt[i]){
		int v=Too[i];
		if(!dfn[v]){
			fa[v][0]=u;
			lmin[v][0]=w[i];
			dep[v]=dep[u]+1;
			dfs(v);
		}
	}
}
inline int LCA(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	for(int i=19;~i;--i)
		if(dep[fa[x][i]]>=dep[y])
			x=fa[x][i];
	if(x==y)return x;
	for(int i=19;~i;--i)
		if(fa[x][i]!=fa[y][i]){
			x=fa[x][i];
			y=fa[y][i];
		}
	return fa[x][0];
}
inline LL find(int u,int f){
	LL ans=INF;
	for(int i=19;~i;--i)
		if(dep[fa[u][i]]>=dep[f]){
			ans=Min(ans,lmin[u][i]);
			u=fa[u][i];
		}
	return ans;
}
inline void insert(int u){
	if(top==1){stk[++top]=u;return;}
	int lca=LCA(u,stk[top]);
	if(lca==stk[top]){stk[++top]=u;return;}
	while(top>1&&dep[lca]<=dep[stk[top-1]]){
		add(stk[top-1],stk[top],find(stk[top],stk[top-1]));
		--top;
	}
	if(lca!=stk[top])add(lca,stk[top],find(stk[top],lca)),stk[top]=lca;
	stk[++top]=u;
}
inline void DP(int u){
	if(ask[u])f[u]=INF;
	else f[u]=0;
	for(int &i=Fir[u];i;i=Nxt[i]){
		int v=Too[i];
		DP(v);
		f[u]+=Min(f[v],w[i]);
	}
}
int main(){
	memset(lmin,0x3f,sizeof(lmin));
	scanf("%d",&n);
	for(int i=1;i<n;++i){
		int u,v;LL w;
		scanf("%d%d%lld",&u,&v,&w);
		add(u,v,w);add(v,u,w);
	}
	dfs(1);
	memset(Fir,0,sizeof(Fir));tot=0;
	scanf("%d",&m);
	for(int i=1;i<=m;++i){
		int k;
		scanf("%d",&k);
		for(int j=1;j<=k;++j)
			scanf("%d",&h[j]),ask[h[j]]=1;
		sort(h+1,h+k+1,cmp);
		stk[++top]=1;
		for(int j=1;j<=k;++j)
			insert(h[j]);
		while(top){
			add(stk[top-1],stk[top],find(stk[top],stk[top-1]));
			--top;
		}
		DP(1);
		printf("%lld\n",f[1]);
		tot=0;
		for(int j=1;j<=k;++j)
			ask[h[j]]=0;
	}
	return 0;
}