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

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

Time Limit: 1 Sec  Memory Limit: 64 MB
Submit: 689  Solved: 295
[Submit][Status]

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

HINT

Source

Gold

题解：

原来居然不会做，是有多弱。。。

今天又来看这题，发现好水，类似与POI 的 STA，记录一个节点的子树有多少头牛，然后任何一个节点的ans可以由它的父节点在O（1）的时间算出

所以一次dfs+一次dp就ok了

代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 100000000000000

#define maxn 150000

#define maxm 500+100

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
struct edge{int go,next;ll w;}e[2*maxn];
int n,tot,head[maxn];
ll sum,s[maxn],f[maxn];
bool v[maxn];
inline void insert(int x,int y,int z)
{
    e[++tot].go=y;e[tot].next=head[x];head[x]=tot;e[tot].w=z;
    e[++tot].go=x;e[tot].next=head[y];head[y]=tot;e[tot].w=z;
}
void dfs(int x)
{
    v[x]=1;
    for(int i=head[x],y;i;i=e[i].next)
    if(!v[y=e[i].go])
    {
        dfs(y);
        s[x]+=s[y];
        f[x]+=f[y]+s[y]*e[i].w;
    }
}
void dp(int x)
{
    v[x]=0;
    for(int i=head[x],y;i;i=e[i].next)
    if(v[y=e[i].go])
    {
        f[y]=f[x]+(sum-s[y]-s[y])*e[i].w;
        dp(y);
    }
}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    n=read();
    for1(i,n)s[i]=read(),sum+=s[i];
    for1(i,n-1){int x=read(),y=read(),z=read();insert(x,y,z);}
    dfs(1);
    dp(1);
    ll ans=inf;
    for1(i,n)ans=min(ans,f[i]);
    printf("%lld\n",ans);

    return 0;

}

inf开到1000亿然后WA了，怒加了几个0就A了T_T