码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj 4557 LNOI2016 Day1

时间:2018-01-06 21:06:06      阅读:153      评论:0      收藏:0      [点我收藏+]

标签:static   name   main   gcc   lin   i++   +=   its   idt   

技术分享图片

SOL

树形DP

f[i][j]表示i的子树中,最高覆盖到i向下第j层的最小花费。

g[i][j]表示i的子树全部覆盖,还能向上覆盖j层的最小花费。

#pragma GCC optimize("-O2")
#include<bits/stdc++.h>
#define sight(c) (‘0‘<=c&&c<=‘9‘)
#define eho(x) for(int i=head[x];i;i=net[i]) 
#define min(a,b) (a)<(b)?(a):(b)
#define DEP 23
#define inf (1<<26)
#define N 500007
int f[N][DEP],g[N][DEP],head[N],fall[N<<1],net[N<<1],w[N],d,tot,n,m,X,A,B;
bool u[N];
void dfs(int x,int fa){
   if (u[x]) f[x][0]=g[x][0]=w[x];
   for (int i=1;i<=d;i++) f[x][i]=w[x]; f[x][d+1]=inf;
   eho(x) if (fall[i]^fa) {
     dfs(fall[i],x);
     for (int j=d;~j;j--) 
         f[x][j]=min(f[x][j]+g[fall[i]][j],g[x][j+1]+f[fall[i]][j+1]);
     for (int j=d;~j;j--)
         f[x][j]=min(f[x][j],f[x][j+1]);
     g[x][0]=f[x][0];
     for (int j=1;j<=d;j++)
         g[x][j]+=g[fall[i]][j-1];
     for (int j=1;j<=d;j++) g[x][j]=min(g[x][j],g[x][j-1]);
   }
}
inline void read(int &x){
    static char c;
    for (c=getchar();!sight(c);c=getchar());
    for (x=0;sight(c);c=getchar()) x=x*10+c-48;
}
inline void add(int x,int y) {
    fall[++tot]=y;net[tot]=head[x];head[x]=tot;
}
void write(int x){
    if (x<10) { putchar(0+x);return;} write(x/10),putchar(0+x%10);
}
inline void writeln(int x){
    if (x<0) putchar(-),x*=-1; write(x); putchar(\n);
} 
using namespace std;
int main () {
//  freopen("a.in","r",stdin);
    read(n); read(d);
    for (int i=1;i<=n;i++) read(w[i]);
    read(m);
    for (int i=1;i<=m;i++) read(X),u[X]=1;
    for (int i=1;i< n;i++) read(A),read(B),
     add(A,B),add(B,A);
    dfs(1,0);
    writeln(f[1][0]); return 0;
}

 

bzoj 4557 LNOI2016 Day1

标签:static   name   main   gcc   lin   i++   +=   its   idt   

原文地址:https://www.cnblogs.com/rrsb/p/8215178.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!