rivers ioi2005 树形dp

说句实话，写完这道题，很想吐一口血出来，以示我心情的糟糕；

题目很简单，树形dp，正常做30分钟，硬是做了好几个小时，真是伤心。

题解不写了，只是吐个槽，网上没有用背包写的dp，全是左儿子右兄弟写法，没法对照一下，难受啊。

代码：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 using namespace std;
 #define LL long long
 int n,K;
 const int maxn=;
 const LL inf=3000000000LL;
 struct node{
     int y,next,v;
 }e[maxn];
 int linkk[maxn],len=,w[maxn],d[maxn][maxn],siz[maxn],q[maxn],tail=,fa[maxn],ru[maxn],vis[maxn];
 int read(){
     int x=;char ch=getchar();bool flag=;
     while(ch<''||ch>''){if(ch=='-')flag=;ch=getchar();}
     while(ch<=''&&ch>=''){x=x*+ch-'';ch=getchar();}
     return flag?-x:x;
 }
 void insert(int x,int y,int v){
     e[++len].y=y;
     e[len].v=v;
     e[len].next=linkk[x];
     linkk[x]=len;
 }
 void print(int x){printf("%d\n",x);}
 void print(int x,int y){printf("%d %d\n",x,y);}
 void init(){
     int y,v;
     n=read(),K=read();
     memset(d,,sizeof(d));
     for(int i=;i<=n;i++){
         w[i]=read(),y=read(),v=read();
         insert(i,y,v);insert(y,i,v);
         d[i][y]=v,d[y][i]=v;
     }

 }
 void findd(){
     for(int i=;i<=n;i++)d[i][i]=;
     for(int k=;k<=n;k++)
         for(int i=;i<=n;i++)
             for(int j=;j<=n;j++)
                 if(d[i][k]+d[k][j]<d[i][j])d[i][j]=d[i][k]+d[k][j];
 }
 LL f[][][];
 void dfs(int x,int f){
     siz[x]=;
     fa[x]=f;
     for(int i=linkk[x];i;i=e[i].next){
         if(e[i].y==f)continue;
         dfs(e[i].y,x);
         siz[x]+=siz[e[i].y];
         ru[x]++;
     }
     if(siz[x]==)q[++tail]=x;
 }
 void work(){
     memset(f,-,sizeof(f));
     findd();
     dfs(,);
     for(int head=;head<=tail;head++){
         int x=q[head];
         ru[fa[x]]--;
         if(ru[fa[x]]==)q[++tail]=fa[x];
         LL g[];
         for(int k=;k<=siz[x]&&k<=K;k++){
             for(int prev=fa[x];prev!=-;prev=fa[prev]){

                 for(int p=;p<=k;p++)g[p]=;
                 for(int i=linkk[x];i;i=e[i].next){
                     if(e[i].y==fa[x])continue;
                     LL minn=inf;
                     for(int j=k;j>=;j--){
                         minn=inf;
                         for(int l=;l<=siz[e[i].y]&&l<=j;l++)
                             minn=min(minn,g[j-l]+f[e[i].y][l][prev]);
                         g[j]=minn;
                     }
                 }
                 LL ans=g[k]+w[x]*d[prev][x];
                 for(int p=;p<=k;p++)g[p]=;
                 for(int i=linkk[x];i;i=e[i].next){
                     if(e[i].y==fa[x])continue;
                     LL minn=inf;
                     for(int j=k-;j>=;j--){
                         minn=inf;
                         for(int l=;l<=siz[e[i].y]&&l<=j;l++)
                             minn=min(minn,g[j-l]+f[e[i].y][l][x]);
                         g[j]=minn;
                     }
                 }
                 f[x][k][prev]=min(ans,k-<?inf:g[k-]);
                 if(!prev)break;
             }
         }
     }
     cout<<f[][K][]<<endl;
 }
 int main(){
     freopen("1.in","r",stdin);
     freopen("1.out","w",stdout);
     init();
     work();
 }