4557: [JLoi2016]侦察守卫

4557: [JLoi2016]侦察守卫

Time Limit: 20 Sec   Memory Limit: 256 MB
Submit: 232   Solved: 161
[ Submit][ Status][ Discuss]

Description

小R和B神正在玩一款游戏。这款游戏的地图由N个点和N-1条无向边组成，每条无向边连接两个点，且地图是连通的

。换句话说，游戏的地图是一棵有N个节点的树。游戏中有一种道具叫做侦查守卫，当一名玩家在一个点上放置侦

查守卫后，它可以监视这个点以及与这个点的距离在D以内的所有点。这里两个点之间的距离定义为它们在树上的

距离，也就是两个点之间唯一的简单路径上所经过边的条数。在一个点上放置侦查守卫需要付出一定的代价，在不

同点放置守卫的代价可能不同。现在小R知道了所有B神可能会出现的位置，请你计算监视所有这些位置的最小代价

。

Input

第一行包含两个正整数N和D，分别表示地图上的点数和侦查守卫的视野范围。约定地图上的点用1到N的整数编号。

第二行N个正整数，第i个正整数表示在编号为i的点放置侦查守卫的代价Wi。保证Wi≤1000。第三行一个正整数M，

表示B神可能出现的点的数量。保证M≤N。第四行M个正整数，分别表示每个B神可能出现的点的编号，从小到大不

重复地给出。接下来N–1行，每行包含两个正整数U,V，表示在编号为U的点和编号为V的点之间有一条无向边。N<=

500000,D<=20

Output

 仅一行一个整数，表示监视所有B神可能出现的点所需要的最小代价

Sample Input

12 2
8 9 12 6 1 1 5 1 4 8 10 6
10
1 2 3 5 6 7 8 9 10 11
1 3
2 3
3 4
4 5
4 6
4 7
7 8
8 9
9 10
10 11
11 12

Sample Output

10

HINT

Source

[ Submit][ Status][ Discuss]
定义状态f[i][j]:以i为根的子树中状态为j的最优方案 其中，j为负数，则可以向外覆盖j单位距离的B神 j为0则刚好完成所有i子树中B神的覆盖 j大于零则该子树中离根最远的未被覆盖B神的距离为j 转移的话，每个儿子的信息逐个合并到根里 O(n*D^2)，，细节就不提了。。之前写过一篇类似的blog，cf的某题树形dp

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<bitset>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

const int maxn = 5E5 + 50;
const int INF = ~0U>>1;

int n,d,m,f[maxn][50],w[maxn],g[50];
bool Mark[maxn];

vector <int> v[maxn];

int F(const int &x) {return x + d;}

void Merge(int x,int y)
{
	for (int a = -d; a <= d; a++)
	{
		if (f[x][F(a)] == INF) continue;
		for (int b = -d; b <= d; b++)
		{
       		if (f[y][F(b)] == INF) continue;
			int nex,now = f[x][F(a)] + f[y][F(b)];
			if (a <= 0 && b <= 0) nex = min(a,b + 1);
			else if (a > 0 && b > 0) 
			{
				if (b == d) continue;
				nex = max(a,b + 1);
			}
			else if (a <= 0 && b > 0) 
			{
				if (-a < b && b == d) continue;
				nex = -a >= b?a:b + 1;
			}
			else nex = -b >= a?b + 1:a;
			g[F(nex)] = min(g[F(nex)],now);
		}                          		
	}
	memcpy(f[x],g,sizeof(g));
}

void Dfs(int x,int from)
{
	f[x][F(-d)] = w[x];
	if (Mark[x]) f[x][F(1)] = 0;
	else f[x][F(0)] = 0;
	for (int i = 0; i < v[x].size(); i++)
	{
		int y = v[x][i];
		if (y == from) continue;
		Dfs(y,x); 
		for (int A = -d; A <= d; A++) g[F(A)] = INF;
		Merge(x,y);
	}
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	cin >> n >> d;
	for (int i = 1; i <= n; i++)
		for (int j = -d; j <= d + 1; j++)
			f[i][F(j)] = INF;
	for (int i = 1; i <= n; i++) scanf("%d",&w[i]);
	cin >> m;
	for (int i = 1; i <= m; i++)
	{
		int x; scanf("%d",&x);
		Mark[x] = 1;
	}
	for (int i = 1; i < n; i++)
	{
		int x,y; scanf("%d%d",&x,&y);
		v[x].push_back(y);
		v[y].push_back(x);
	}
	
	int ans = ~0U>>1; Dfs(1,0);
	for (int i = -d; i <= 0; i++) 
		ans = min(ans,f[1][F(i)]);
	cout << ans;
	return 0;
}