BZOJ 1827: [Usaco2010 Mar]gather 奶牛大集会 树形DP + 带权重心

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

* 第一行：一个值，表示最小的不方便值。

题解：不难发现选定点必须是树的带权重心 (任意儿子大小不能超过总和一半). 

直接 DFS 求重心并算一下贡献即可. 

#include <bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin) 
#define maxn 1000000 
#define inf 100000000000000 
#define ll long long 
using namespace std; 
ll f[maxn], siz[maxn], sumv[maxn]; 
int C[maxn], hd[maxn], to[maxn << 1], nex[maxn << 1], val[maxn << 1]; 
int n, edges, root; 
ll tot; 
void add(int u, int v, int c)
{
    nex[++edges] = hd[u], hd[u] = edges, to[edges] = v, val[edges] = c; 
}
void dfs(int u, int ff)
{
    siz[u] = C[u], f[u] = 0; 
    for(int i = hd[u]; i ; i = nex[i]) 
    {
        int v = to[i]; 
        if(v == ff) continue; 
        dfs(v, u);         
        siz[u] += siz[v];
        f[u] = max(f[u], siz[v]); 
    }               
    f[u] = max(f[u], tot - siz[u]); 
    if(f[u] < f[root]) root = u; 
}      
void calc(int u,int ff)
{
    siz[u] = C[u], sumv[u] = 0; 
    for(int i = hd[u]; i ; i = nex[i])
    {
        int v = to[i]; 
        if(v == ff) continue; 
        calc(v, u); 
        siz[u] += siz[v]; 
        sumv[u] += sumv[v] + val[i] * siz[v]; 
    }
}
int main()
{
    // setIO("input"); 
    scanf("%d",&n); 
    for(int i = 1; i <= n; ++i) scanf("%d",&C[i]), tot += C[i]; 
    for(int i = 1, u, v, c; i < n; ++i) 
    {
        scanf("%d%d%d",&u,&v,&c); 
        add(u, v, c); 
        add(v, u, c); 
    }
    f[0] = inf, root = 0, dfs(1, 0);
    memset(siz, 0, sizeof(siz)), calc(root, 0);  
    printf("%lld\n",sumv[root]); 
    return 0; 
}