BZOJ 4557 侦查守卫

好迷的树形dp。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxv 500500
#define maxe 1000500
using namespace std;
int n,m,d,x,y,nume=,g[maxv],dp1[maxv][],dp2[maxv][],w[maxv],fath[maxv];
bool mark[maxv][];
struct edge
{
    int v,nxt;
}e[maxe];
void addedge(int u,int v)
{
    e[++nume].v=v;
    e[nume].nxt=g[u];
    g[u]=nume;
}
void dp(int x)
{
    int flag=;
    for (int i=g[x];i;i=e[i].nxt)
    {
        int v=e[i].v;
        if (v!=fath[x])
        {
            fath[v]=x;
            dp(v);flag=;
            for (int j=;j<=d;j++)
                mark[x][j]|=mark[v][j-];
        }
    }
    if (!flag)
    {
        if (mark[x][]) dp1[x][]=dp2[x][]=w[x];
        for (int i=;i<=d;i++)
        {
            dp2[x][i]=w[x];
            dp1[x][i]=;
        }
        return;
    }
    for (int i=;i<=d;i++)
        for (int j=g[x];j;j=e[j].nxt)
        {
             int v=e[j].v;
             if (v!=fath[x])
                dp1[x][i]+=dp1[v][i-];
        }
    dp2[x][d]+=w[x];
    for (int i=g[x];i;i=e[i].nxt)
    {
        int v=e[i].v;
        if (v!=fath[x])
            dp2[x][d]+=dp1[v][d];
    }
    for (int i=d-;i>=;i--)
    {
        int ret=;
        for (int j=g[x];j;j=e[j].nxt)
        {
             int v=e[j].v;
             if (v!=fath[x])
                ret+=dp1[v][i];
        }
        dp2[x][i]=dp2[x][i+];
        for (int j=g[x];j;j=e[j].nxt)
        {
             int v=e[j].v;
             if (v!=fath[x])
                dp2[x][i]=min(dp2[x][i],dp2[v][i+]+ret-dp1[v][i]);
        }
    }
    dp1[x][]=dp2[x][];
    for (int i=;i<=d;i++)
        dp1[x][i]=min(dp1[x][i],dp1[x][i-]);
    for (int i=d-;i>=;i--)
        if (!mark[x][i]) dp1[x][i]=min(dp1[x][i],dp1[x][i+]);
    dp2[x][]=dp1[x][];
}
int main()
{
    scanf("%d%d",&n,&d);
    for (int i=;i<=n;i++)
        scanf("%d",&w[i]);
    scanf("%d",&m);
    for (int i=;i<=m;i++)
    {
        scanf("%d",&x);
        mark[x][]=true;
    }
    for (int i=;i<=n-;i++)
    {
        scanf("%d%d",&x,&y);
        addedge(x,y);addedge(y,x);
    }
    dp();
    printf("%d\n",dp1[][]);
    return ;
}