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
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 }