bzoj 4557 [JLoi2016]侦察守卫

http://www.elijahqi.win/archives/3698
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
好神 我好菜

设f[x][i]表示 x的子树 由外部覆盖了x子树的前i层的最小代价 x为第一层

那么考虑f[x][i]=sigmaf[y][i-1]

但是考虑不仅是父节点 兄弟节点也会互相影响 于是增加g[x][i]表示x节点的子树完全被覆盖

x是第0层 现在可以额外覆盖子树以外i层的东西

那么转移就看代码了～

#include<cstdio>
#include<cctype>
#include<algorithm>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int N=5e5+10;
const int inf=0x3f3f3f3f;
int n,m,w[N],q,h[N],f[N][22],g[N][22],num;
bool mark[N];
struct node{
    int y,next;
}data[N<<1];
inline void dfs(int x,int fa){
    g[x][0]=f[x][0]=mark[x]?w[x]:0;
    for (int i=1;i<=m;++i) g[x][i]=w[x];g[x][m+1]=inf;
    for (int owo=h[x];owo;owo=data[owo].next){
        int y=data[owo].y;if (y==fa) continue;dfs(y,x);
        for (int i=0;i<=m;++i) g[x][i]=min(g[x][i]+f[y][i],f[x][i+1]+g[y][i+1]);
        for (int i=m;~i;--i) g[x][i]=min(g[x][i],g[x][i+1]);
        f[x][0]=g[x][0];
        for (int i=1;i<=m;++i) f[x][i]+=f[y][i-1];
        for (int i=1;i<=m;++i) f[x][i]=min(f[x][i],f[x][i-1]);
    }
}
int main(){
    freopen("bzoj4557.in","r",stdin);
    n=read();m=read();
    for (int i=1;i<=n;++i) w[i]=read();q=read();
    for (int i=1;i<=q;++i) mark[read()]=1;
    for (int i=1;i<n;++i){
        int x=read(),y=read();
        data[++num].y=y;data[num].next=h[x];h[x]=num;
        data[++num].y=x;data[num].next=h[y];h[y]=num;
    }dfs(1,0);printf("%d\n",f[1][0]);
    return 0;
}