【树形DP/搜索】BZOJ 1827: [Usaco2010 Mar]gather 奶牛大集会

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

Time Limit: 1 Sec Memory Limit: 64 MB
Submit: 793 Solved: 354
[Submit][Status][Discuss]

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

　　先求每个点的子树点数

　　然后再求设根节点为集合点的大小

　　然后一路搜下去，更新权值

　　感觉又可以解决一类问题了。。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 #include<cmath>
 5  
 6 #define maxn 100001
 7  
 8 using namespace std;
 9  
10 long long tot=0,w[maxn],f[maxn],ans;
11  
12 inline int in()
13 {
14     int x=0;char ch=getchar();
15     while(ch<'0'||ch>'9')ch=getchar();
16     while(ch<='9'&&ch>='0')x=x*10+ch-'0',ch=getchar();
17     return x;
18 }
19  
20 struct ed{
21     int to,c,last;
22 }edge[maxn*2];
23  
24 int last[maxn],cnt=0,a[maxn];
25  
26 bool vis[maxn];
27  
28 void add(int u,int v,int c)
29 {
30     edge[++cnt].to=v,edge[cnt].c=c,edge[cnt].last=last[u],last[u]=cnt;
31     edge[++cnt].to=u,edge[cnt].c=c,edge[cnt].last=last[v],last[v]=cnt;
32 }
33  
34 void DP(int pos)
35 {
36     w[pos]=a[pos];vis[pos]=1;
37     for(int i=last[pos];i;i=edge[i].last)
38     {
39         int u=edge[i].to;
40         if(vis[u])continue;
41         DP(u);
42         w[pos]+=w[u];
43         f[pos]+=w[u]*edge[i].c+f[u];
44     }
45 }
46  
47 void dfs(int poi,long long val)
48 {
49     vis[poi]=1;ans=min(ans,val);
50     for(int i=last[poi];i;i=edge[i].last)
51         if(!vis[edge[i].to])
52             dfs(edge[i].to,val+(tot-2*w[edge[i].to])*edge[i].c);
53 }
54  
55 int main()
56 {
57     int n,u,v,c;
58     n=in();
59     for(int i=1;i<=n;i++)a[i]=in(),tot+=a[i];
60     for(int i=1;i<n;i++)
61     {
62         u=in();v=in();c=in();
63         add(u,v,c);
64     }
65     DP(1);ans=f[1];
66     memset(vis,0,sizeof(vis));
67     dfs(1,f[1]);
68     printf("%lld",ans);
69     return 0;
70 }

View Code