P1099 [NOIP2007 提高组] 树网的核
题面
设 \(T=(V,E,W)\) 是一个无圈且连通的无向图(也称为无根树),每条边都有正整数的权,我们称 \(T\) 为树网(treenetwork
),其中 \(V\),\(E\) 分别表示结点与边的集合,\(W\) 表示各边长度的集合,并设 \(T\) 有 \(n\) 个结点。
路径:树网中任何两结点 \(a\),\(b\) 都存在唯一的一条简单路径,用 \(d(a, b)\) 表示以 \(a, b\) 为端点的路径的长度,它是该路径上各边长度之和。我们称
\(d(a, b)\) 为 \(a, b\) 两结点间的距离。
\(D(v, P)=\min\{d(v, u)\}\), \(u\) 为路径 \(P\) 上的结点。
树网的直径:树网中最长的路径成为树网的直径。对于给定的树网 \(T\),直径不一定是唯一的,但可以证明:各直径的中点(不一定恰好是某个结点,可能在某条边的内部)是唯一的,我们称该点为树网的中心。
偏心距 \(\mathrm{ECC}(F)\):树网 \(T\) 中距路径 \(F\) 最远的结点到路径 \(F\) 的距离,即
\(\mathrm{ECC}(F)=\max\{D(v, F),v \in V\}\)
任务:对于给定的树网 \(T=(V, E, W)\) 和非负整数 \(s\),求一个路径 \(F\),他是某直径上的一段路径(该路径两端均为树网中的结点),其长度不超过 \(s\)(可以等于 \(s\)),使偏心距 \(\mathrm{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\)。
求指定意义下的最小偏心距。
【数据范围】
- 对于 \(40\%\) 的数据,保证 \(n \le 15\)。
- 对于 \(70\%\) 的数据,保证 \(n \le 80\)。
- 对于 \(100\%\) 的数据,保证 \(n \le 300\),\(0\le s\le10^3\),\(1 \leq u, v \leq n\),\(1 \leq w \leq 10^3\)。
简化题意
给定一棵带边权无根树,在其直径上求出一段长度不超过 \(s\) 的路径 \(F\),使得离路径距离最远的点到路径的距离最短。
思路
首先求出树的直径。可以使用众所周知的两遍DFS大法。
先选择任意一个顶点为源点,然后DFS求出最远的点 \(k\),然后再进行一遍这个操作。得到的就是直径的两个端点。
然后用单调队列的思想,再尺取求答案的最小值。
把直径标记出来,重新计算贡献(偏心距),然后统计答案。
代码
#include <bits/stdc++.h>
using namespace std;
int n,m,k,ans=INT_MAX;
int dis[500005],fa[500005],head[500005];
bool flags[500005];
struct edge{
int next,to,weight;
} e[500005*2];
int ec;
void add(int u,int v,int w){
e[++ec].to=v;
e[ec].weight=w;
e[ec].next=head[u];
head[u]=ec;
}
void dfs(int father,int x){
fa[x]=father;
if(dis[x]>dis[k]){
k=x;
}
for(int i=head[x];i;i=e[i].next){
int y=e[i].to;
if(y==father||flags[y]){
continue;
}
dis[y]=dis[x]+e[i].weight;
dfs(x,y);
}
}
int main(){
cin>>n>>m;
for(int i=1,x,y,z;i<n;i++){
cin>>x>>y>>z;
add(x,y,z);
add(y,x,z);
}
dis[1]=1;
dfs(0,1);
dis[k]=0;
dfs(0,k);
int top=k;
for(int i=top,j=top,l=1,r=0;i;i=fa[i]){
while(dis[j]-dis[i]>m){
j=fa[j];
}
int x=max(dis[top]-dis[j],dis[i]);
ans=min(ans,x);
}
for(int i=top;i;i=fa[i]){
flags[i]=true;
}
for(int i=top;i;i=fa[i]){
k=i;
dis[k]=0;
dfs(fa[i],i);
}
for(int i=1;i<=n;i++){
ans=max(ans,dis[i]);
}
cout<<ans<<endl;
return 0;
}