1827: [Usaco2010 Mar]gather 奶牛大集会

Description

Bessie正在计划一年一度的奶牛大集会,来自全国各地的奶牛将来参加这一次集会。当然,她会选择最方便的地点来举办这次集会。每个奶牛居住在 N(1<=N<=100,000) 个农场中的一个,这些农场由N-1条道路连接,并且从任意一个农场都能够到达另外一个农场。道路i连接农场A_i和B_i(1 <= A_i <=N; 1 <= B_i <= N),长度为L_i(1 <= L_i <= 1,000)。集会可以在N个农场中的任意一个举行。另外,每个牛棚中居住者C_i(0 <= C_i <= 1,000)只奶牛。在选择集会的地点的时候,Bessie希望最大化方便的程度(也就是最小化不方便程度)。比如选择第X个农场作为集会地点,它的不方便程度是其它牛棚中每只奶牛去参加集会所走的路程之和,(比如,农场i到达农场X的距离是20,那么总路程就是C_i*20)。帮助Bessie找出最方便的地点来举行大集会。 考虑一个由五个农场组成的国家,分别由长度各异的道路连接起来。在所有农场中,3号和4号没有奶牛居住。 

Input

第一行:一个整数N * 第二到N+1行:第i+1行有一个整数C_i * 第N+2行到2*N行,第i+N+1行为3个整数:A_i,B_i和L_i。

Output

* 第一行:一个值,表示最小的不方便值。

Sample Input

5
1
1
0
0
2
1 3 1
2 3 2
3 4 3
4 5 3

Sample Output

15

HINT

 

Source

 
 
树形DP写起来也是很有趣的啊,2333333
其实算是一道裸的树形DP了。。。
先dp一次求出每个点为集会点时它的儿子节点的花费
然后再dp一次求出每个点为集会点时他的父亲节点的花费。。
然后取min。。。
注意会爆int
 4 #include<iostream>
 5 #include<cstdlib>
 6 #include<cmath>
 7 #include<cstring>
 8 #include<cstdio>
 9 #include<algorithm>
10 #include<string>
11 #include<map>
12 #include<queue>
13 #include<vector>
14 #include<set>
15 #define inf 1LL<<50
16 #define maxn 100005
17 #define maxm 100005
18 #define eps 1e-10
19 #define ll long long
20 #define for0(i,n) for(int i=0;i<=(n);i++)
21 #define for1(i,n) for(int i=1;i<=(n);i++)
22 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
23 #define for3(i,x,y) for(int i=(x);i>=(y);i--)
24 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
25 using namespace std;
26 struct edge{
27     int go,next,w;
28 }e[maxn*2];
29 int tot,a[maxn],head[maxn],f[maxn];
30 ll sum[maxn],SUM,num[maxn],ans[maxn];
31 void insert(int x,int y,int z){
32     e[++tot].go=y;e[tot].next=head[x];head[x]=tot;e[tot].w=z;
33 }
34 void ins(int x,int y,int z){
35     insert(x,y,z);insert(y,x,z);
36 }
37 void dp1(int x,int fa){
38     num[x]=a[x];
39     sum[x]=0;
40     for(int i=head[x];i;i=e[i].next){
41         if(e[i].go==fa)continue;
42         dp1(e[i].go,x);
43         num[x]+=num[e[i].go];
44         sum[x]+=sum[e[i].go]+num[e[i].go]*e[i].w; 
45     }
46 }
47 void dp2(int x,int fa,int k){
48     if(!fa)ans[x]=sum[x];
49     else ans[x]=ans[fa]-num[x]*k+(SUM-num[x])*k;
50     for(int i=head[x];i;i=e[i].next){
51         if(fa==e[i].go)continue;
52         dp2(e[i].go,x,e[i].w);
53     }
54 }
55 int read(){
56     int x=0,f=1;char ch=getchar();
57     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
58     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
59     return x*f;
60 }
61 int main(){
62     //freopen("input.txt","r",stdin);
63     //freopen("output.txt","w",stdout);
64     int n=read();
65     for1(i,n)a[i]=read(),SUM+=a[i];
66     for1(i,n-1){
67         int u=read(),v=read(),z=read();
68         ins(u,v,z);
69     }
70     dp1(1,0);
71     dp2(1,0,0);
72     ll mi=inf;
73     for1(i,n)mi=min(mi,ans[i]);
74     cout<<mi<<endl; 
75     return 0;
76 }
View Code

 

posted @ 2016-03-12 16:35  HTWX  阅读(166)  评论(0编辑  收藏  举报