JZOJ.1729. blockenemy

Problem

Description

你在玩电子游戏的时候遇到了麻烦……你玩的游戏是在一个虚拟的城市里进行，这个城市里有n个点，都从0~n-1编了号，每两个点之间有且仅有一条路径。现在，你的敌人到这个城市来踩点了!!!为了阻止他们更好的踩点， 你决定切断他们所有踩点人员的联系，使他们孤军作战，然后在各个击破。但是这就要切断某些街道，而你每切断一条路，市民就会产生相对的不满值，不满值越大，城市的和谐度就越小。所以你现在需要知道为了使踩点人员所在的点两两之间不联通所切断的边产生的最小不满值是多少？

Input

第一行一个数：n n<=50 以下n-1行，每行3个数 a，b，c 表示a点和b点之间有条路，切断这条路的不满值为c 以下若干行 每行一个数，表示踩点人员的位

Output

一个数，最小不满值

Sample Input

5 1 0 1 1 2 2 0 3 3 4 0 4 3 2 4

Sample Output

4

Solution

我8点多钟的时候已经想到了方法。我们设树的重心为根，咱们先来画个图。

那么我们要求出任意一对敌人的LCA（枚举吧！反正N最大才50，1/2*N^2也不大）
那么当求到(A,B)和(A,E)时，如果边#1是(A,B)最短路径的最短边，边#2是(A,E)最短路径的最短边，且边#1大于边#2，那么我们要删掉边#1和边#2，而不能保留#2.
原因：如果去掉边#2，到时候求(B,E)最短路径的最短边时我们又会求到边#2（自己脑补一下），所以越早删掉越好。注意一下判重，一条边只能被删掉一次。
然后直接计算被删掉边的权值的总和就OK了。

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define N 60
using namespace std;
int head[N*2],next[N*2],go[N*2],val[N*2],deep[N],size[N],mx[N],a[N],fa[N][2];
bool b[N][N];
int i,j,k,n,tot,x,y,z,hv,od,ans;
void lb(int x,int y,int z)
{
    go[++tot]=y;
    next[tot]=head[x];
    head[x]=tot;
    val[tot]=z;
}
void zheavy(int set,int x,int y)
{
    int i;
    mx[x]=0;size[x]=1;
    for (i=head[x];i;i=next[i])
    {
        if (go[i]!=y)
        {
            zheavy(set,go[i],x);
            size[x]+=size[go[i]];
            mx[x]=max(mx[x],size[go[i]]);
        }
    }
    mx[x]=max(mx[x],set-size[x]); 
    if (mx[x]<mx[hv]) hv=x; 
    return;
}
void dfs(int x,int y,int z)
{
    int i;
    deep[x]=z;
    for(i=head[x];i;i=next[i])
    {
        if (go[i]!=y)
        {
            dfs(go[i],x,z+1);
            fa[go[i]][0]=x;
            fa[go[i]][1]=val[i];
        }
    }
}
int lca(int x,int y)
{
    int mn=2147483647,w1,w2;
    if (deep[x]<deep[y]) swap(x,y);
    while (deep[x]>deep[y])
    {
        if (fa[x][1]<mn)
        {
            mn=fa[x][1];
            w1=x;
            w2=fa[x][0];
        }
        x=fa[x][0];
    }
    if (x==y) 
    {
        if (b[w1][w2]) 
        {
            b[w1][w2]=b[w2][w1]=0;
            return mn;
        } else return 0;
    }
    while (x!=y)
    {
        if (fa[x][1]<mn)
        {
            mn=fa[x][1];
            w1=x;
            w2=fa[x][0];
        }
        x=fa[x][0];
        if (fa[y][1]<mn)
        {
            mn=fa[y][1];
            w1=y;
            w2=fa[y][0];
        }
        y=fa[y][0];
    }
    if (b[w1][w2]) 
    {
        b[w1][w2]=b[w2][w1]=0;
        return mn;
    } else return 0;
}
int main()
{
    scanf("%d",&n);
    for (i=1;i<=n-1;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        lb(x,y,z);
        lb(y,x,z);
        b[x][y]=b[y][x]=1;
    }
    while (scanf("%d",&a[od+1])!=EOF) od++;
    mx[51]=2147438647;
    zheavy(n,0,51);
    dfs(hv,0,51);
    for (i=1;i<=od-1;i++)
        for (j=i+1;j<=od;j++) ans+=lca(a[i],a[j]);
    printf("%d",ans);
}

其实这道题还有很多方法的（如贪心），欢迎各位大神讨论。

——2016.7.14