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;
}