bzoj 4711 小奇挖矿 —— 树形DP

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4711

就是树形DP，然而也想了半天才把转移想清楚；

f[x][j][0] 表示 x 去上面 j 步的仓库，f[x][j][1] 表示 x 去子树内下去 j 步的仓库而且整个子树都算好了；

f[x][j][0] 就是 d[j] 加上儿子们的 f[u][...][1] 与 f[u][j+1][0] 的最小值；

f[x][j][1] 就是选一个儿子加它的 f[u][j-1][1]，其它儿子是 f[u][...][1] 与 f[u][j+1][0] 的最小值；

f[x][0][1] 要特殊处理一下；

然后一交秒WA，又改了改，似乎也没改什么，就A了...总之注意边界。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int const xn=,inf=0x3f3f3f3f;
int n,K,d[xn],ans,f[xn][xn][],mn[xn];
int hd[xn],ct,nxt[xn<<],to[xn<<];
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
  return f?ret:-ret;
}
void add(int x,int y){to[++ct]=y; nxt[ct]=hd[x]; hd[x]=ct;}
void dfs(int x,int fa)
{
  for(int i=hd[x],u;i;i=nxt[i])
    {
      if((u=to[i])==fa)continue;
      dfs(u,x);
    }
  for(int j=;j<n;j++)
    {
      if(!j)f[x][j][]=K;
      else f[x][j][]=d[j];
      for(int i=hd[x],u;i;i=nxt[i])
        {
          if((u=to[i])==fa)continue;
//          for(int k=0;k<n;k++)ret=min(ret,f[u][k][1]);
          int ret=min(mn[u],f[u][j+][]);
          f[x][j][]+=ret;
        }
      if(!j)
        {
          f[x][j][]=K;
          for(int i=hd[x],u;i;i=nxt[i])
            {
              if((u=to[i])==fa)continue;
              int ret=min(f[u][][],mn[u]);
//              for(int k=0;k<n;k++)ret=min(ret,f[u][k][1]);
              f[x][j][]+=ret;
            }
        }
      else
        {
          f[x][j][]=inf;
          for(int i=hd[x],v;i;i=nxt[i])
            {
              if((v=to[i])==fa)continue;
              int r=f[v][j-][]+d[j];
              for(int l=hd[x],u;l;l=nxt[l])
                {
                  if((u=to[l])==fa||u==v)continue;
                  int ret=min(f[u][j+][],mn[u]);
//                  for(int k=0;k<n;k++)ret=min(ret,f[u][k][1]);
                  r+=ret;
                }
              f[x][j][]=min(f[x][j][],r);
            }
        }
        mn[x]=min(mn[x],f[x][j][]);
    }
}
int main()
{
  n=rd(); K=rd();
  for(int i=;i<n;i++)d[i]=rd();
  for(int i=,x,y;i<n;i++)
    {
      x=rd(); y=rd();
      add(x,y); add(y,x);
    }
  memset(f,0x3f,sizeof f);
  memset(mn,0x3f,sizeof mn);
  dfs(,); ans=f[][][];
  for(int j=;j<n;j++)ans=min(ans,f[][j][]);
  printf("%d\n",ans);
  return ;
}