1827: [Usaco2010 Mar]gather 奶牛大集会
Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 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
1
1
0
0
2
1 3 1
2 3 2
3 4 3
4 5 3
Sample Output
15
先求每个点的子树点数
然后再求设根节点为集合点的大小
然后一路搜下去,更新权值
感觉又可以解决一类问题了。。
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 }