BZOJ 1999: [Noip2007]Core树网的核【最长链】

Description

设T=(V, E, W) 是一个无圈且连通的无向图（也称为无根树），每条边带有正整数的权，我们称T为树网（treenetwork），其中V, E分别表示结点与边的集合，W表示各边长度的集合，并设T有n个结点。 路径：树网中任何两结点a,b都存在唯一的一条简单路径，用d(a,b)表示以a,b为端点的路径的长度，它是该路径上各边长度之和。我们称d(a,b)为a,b两结点间的距离。 一点v到一条路径P的距离为该点与P上的最近的结点的距离： d(v，P)=min{d(v，u)，u为路径P上的结点}。 树网的直径：树网中最长的路径称为树网的直径。对于给定的树网T，直径不一定是唯一的，但可以证明：各直径的中点（不一定恰好是某个结点，可能在某条边的内部）是唯一的，我们称该点为树网的中心。 偏心距ECC(F)：树网T中距路径F最远的结点到路径F的距离，即 。 任务：对于给定的树网T=(V, E,W)和非负整数s，求一个路径F，它是某直径上的一段路径（该路径两端均为树网中的结点），其长度不超过s（可以等于s），使偏心距ECC(F)最小。我们称这个路径为树网T=(V,E,W)的核（Core）。必要时，F可以退化为某个结点。一般来说，在上述定义下，核不一定只有一个，但最小偏心距是唯一的。 下面的图给出了树网的一个实例。图中，A-B与A-C是两条直径，长度均为20。点W是树网的中心，EF边的长度为5。如果指定s=11，则树网的核为路径DEFG（也可以取为路径DEF），偏心距为8。如果指定s=0（或s=1、s=2），则树网的核为结点F，偏心距为12。

题解

显然，任意一条直径求得的答案都是一样的，因为如果有多条直径，那么直径一定经过一个公共点。那么，我们只要求出任意一条直径，然后计算出直径上每一个点对答案的贡献，枚举核的一个端点，从所有答案中求最小值就是最小偏心距了。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 500006
using namespace std;
inline char nc(){
    static char buf[100000],*i=buf,*j=buf;
    return i==j&&(j=(i=buf)+fread(buf,1,100000,stdin),i==j)?EOF:*i++;
}
inline int _read(){
    char ch=nc();int sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
    return sum;
}
int n,nn,tt,K,ans,tot,hed,tal,len,dis[maxn],a[maxn],fa[maxn],f[maxn],lnk[maxn],son[maxn*2],nxt[maxn*2],w[maxn*2],que[maxn];
bool vis[maxn];
void add(int x,int y,int z){
    nxt[++tot]=lnk[x];son[tot]=y;w[tot]=z;lnk[x]=tot;
}
void dfs(int x){
    vis[x]=0;
    for(int j=lnk[x];j;j=nxt[j]) if(vis[son[j]]){
        fa[son[j]]=x;f[son[j]]=max(f[son[j]],f[x]+w[j]);
        dfs(son[j]);
    }
}
void dfs2(int x,int sum){
    vis[x]=0;tt=max(tt,sum);
    for(int j=lnk[x];j;j=nxt[j]) if(vis[son[j]])dfs2(son[j],sum+w[j]);
}
int main(){
    freopen("core.in","r",stdin);
    freopen("core.out","w",stdout);
    n=_read();K=_read();
    for(int i=1,x,y,z;i<n;i++)x=_read(),y=_read(),z=_read(),add(x,y,z),add(y,x,z);
    memset(vis,1,sizeof(vis));dfs(1);
    int t=0;
    for(int i=1;i<=n;i++)if(f[i]>f[t])t=i;
    memset(vis,1,sizeof(vis));memset(f,0,sizeof(f));memset(fa,0,sizeof(fa));
    dfs(t);
    for(int i=1;i<=n;i++)if(f[i]>f[t])t=i;
    len=f[t];int x=t;
    memset(vis,1,sizeof(vis));
    while(x)vis[x]=0,x=fa[x];
    x=t;
    while(x){
        tt=0;dfs2(x,0);
        dis[++nn]=tt;a[nn]=x;
        x=fa[x];
    }
    int i=1,j=1;ans=1e9;
    while(i<=nn){
        int sum=len-f[a[i]];
        while(j<=nn&&f[a[i]]-f[a[j]]<=K)sum=max(sum,dis[j]),j++;
        sum=max(sum,f[a[--j]]);ans=min(sum,ans);
        i=j+1;
    }
    printf("%d\n",ans);
    return 0;
}