【BZOJ1827】奶牛大集会（树形递推）

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

题解：如果题目已经给定了集会的地点，那么做一遍dfs就好了。我们发现：如果从当前的集会地点移动到它的任意一个儿子，只会有两个值发生变化，合起来变化量为(siz[root]-2*siz[son])*len[root][son],因此我们可令任意一个节点为集会地点，先算一遍，再O(N)寻找解最优的地方就行了。
代码如下：

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#define ll long long
#define inf 0x7f7f7f7f
#define N 100005
using namespace std;
int n,val[N],x,y,z;
ll dis[N],siz[N],ans;
struct nod
{
    int v,w,nex;
    nod(){}
    nod(int a,int b,int c):v(a),w(b),nex(c){}
}e[N<<1];
int cnt,hd[N<<1];
void add(int x,int y,int z)
{
    e[++cnt]=nod(y,z,hd[x]);hd[x]=cnt;
    e[++cnt]=nod(x,z,hd[y]);hd[y]=cnt;
}
ll dfs(int u,int fa)
{
    ll ans=dis[u]*val[u];
    siz[u]=val[u];
    for(int i=hd[u];i;i=e[i].nex)
    {
        int nex=e[i].v,len=e[i].w;
        if(nex==fa) continue;
        dis[nex]=dis[u]+len;
        ans+=dfs(nex,u);
        siz[u]+=siz[nex];
    }
    return ans;
}
void work(int u,int fa)
{
    for(int i=hd[u];i;i=e[i].nex)
    {
        int nex=e[i].v,len=e[i].w;
        if(nex==fa) continue;
        if(siz[1]-2*siz[nex]<0)
        {
            ans=ans+(siz[1]-2*siz[nex])*len;
            work(nex,u);
        }
    }
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&val[i]);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z);
    }
    ans=dfs(1,0);
    work(1,0);
    printf("%lld\n",ans);
    return 0;
}