[vijos P1518] [tyvj P1506] 河流

背景
安徽省芜湖市集训队练习题
IOI 2005 Rivers(riv)
Description:Amber
Data:Official
Program:JackDavid127

---

描述
几乎整个Byteland 王国都被森林和河流所覆盖。小点的河汇聚到一起，形成了稍大点的河。就这样，所有的河水都汇聚并流进了一条大河，最后这条大河流进了大海。这条大河的入海口处有一个村庄——Bytetown。
在Byteland国，有n个伐木的村庄，这些村庄都座落在河边。目前在Bytetown，有一个巨大的伐木场，它处理着全国砍下的所有木料。木料被砍下后，顺着河流而被运到Bytetown的伐木场。Byteland 的国王决定，为了减少运输木料的费用，再额外地建造k个伐木场。这k个伐木场将被建在其他村庄里。这些伐木场建造后，木料就不用都被送到Bytetown了，它们可以在 运输过程中第一个碰到的新伐木场被处理。显然，如果伐木场座落的那个村子就不用再付运送木料的费用了。它们可以直接被本村的伐木场处理。
注：所有的河流都不会分叉，形成一棵树，根结点是Bytetown。
国王的大臣计算出了每个村子每年要产多少木料，你的任务是决定在哪些村子建设伐木场能获得最小的运费。其中运费的计算方法为：每一吨木料每千米1分钱。
编一个程序： 
1．从文件读入村子的个数，另外要建设的伐木场的数目，每年每个村子产的木料的块数以及河流的描述。 
2．计算最小的运费并输出。

---

格式
输入格式

第一行包括两个数n（2<=n<=100），k（1<=k<=50,且k<=n）。n为村庄数，k为要建的伐木场的数目。除了Bytetown 外，每个村子依次被命名为 1，2，3……n，Bytetown被命名为0。
接下来n行，每行3个整数：
wi——每年 i 村子产的木料的块数。（0<=wi<=10000） 
vi——离 i 村子下游最近的村子。（即 i 村子的父结点）（0<=vi<=n） 
di——vi 到 i 的距离(千米)。（1<=di<=10000） 
保证每年所有的木料流到bytetown 的运费不超过2000,000,000分 
50％的数据中n不超过20。
输出格式

输出最小花费，精确到分。
样例1
样例输入1[复制]
4 2 
1 0 1 
1 1 10 
10 2 5 
1 2 3 
样例输出1[复制]
4

提示
树形动态规划
经典问题
来源
安徽省芜湖市集训队练习题
IOI 2005 Rivers(riv)
Description:Amber
Data:Official
Program:JackDavid127

---

首先多叉转二叉，lch为真儿子，rch为兄弟；
记忆化dfs（root，groot，kk，s） 表示以root为根，最近的有伐木场的祖先为groot，且root的子孙们（不算自己）建造kk个伐木场的最优解。s表示的是root的父节点距离groot的距离，而不是root距groot的距离（参照多叉转二叉，因为若s表示的是root到groot的距离的话，其rch即兄弟节点到groot的距离就无法表示。）本程序中cout<

#include <iostream>
#include <vector>
#include <cstring>
using namespace std;
int fa[111],len[111],a[111],lch[111],rch[111];
int f[111][111][111];
int n,k,h,t;
int dfs(int root,int groot,int kk,int s)
{
    if(root==-1||groot==-1) return 0;
    if(f[root][groot][kk]!=2139062143) return f[root][groot][kk];
    int &now=f[root][groot][kk];
    for(int i=0;i<=kk;i++) // root这个点不放伐木场
        now=min(now,dfs(lch[root],groot,i,s+len[root])+dfs(rch[root],groot,kk-i,s)+(s+len[root])*a[root]);
    for(int i=0;i<kk;i++)//放   伐木场
        now=min(now,dfs(lch[root],root,i,0)+dfs(rch[root],groot,kk-1-i,s));
    return f[root][groot][kk];
}

int main()
{
    memset(rch,-1,sizeof(rch));
    memset(lch,-1,sizeof(lch));
    memset(f,127,sizeof(f));
    //cout<<f[1][1][1]<<endl;// 2139062143
    cin>>n>>k;
    for(int i=1;i<=n;i++)
    {
        cin>>a[i]>>fa[i]>>len[i];
        rch[i]=lch[fa[i]];
        lch[fa[i]]=i;
    }
    // cout<<' '<<rch[0]<<endl;
    cout<<dfs(0,0,k,0)<<endl;
    //cout<<dfs(lch[0],0,k,0)<<endl;
}