【NOIP2007提高组】树网的核

题目背景

NOIP2007提高组试题4。

题目描述

设 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 的距离，即：
ECC（F）=max{d(v,F),v∈V}。
任务：对于给定的树网 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。

输入格式

输入包含 n 行：
第 1 行：两个整数 n 和 s，中间用一个空格隔开。其中 n 为树网的核的长度的上界。设结点编号依次为1，2，…，n 。
从第 2 行到第 n 行，每行给出 3 个用空格隔开的正整数，依次表示每一条边的两个端点编号和长度。例如，“2 4 7 ”表示连接结点 2 与 4 的边的长度为 7。
所给的数据都是正确的，不必检验。

输出格式

输出只有一个非负整数，为指定意义下的最小偏心距。
样例数据 1
输入

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

输出

5

样例数据 2
输入

8 6
1 3 2
2 3 2
3 4 6
4 5 3
4 6 4
4 7 2
7 8 3

输出

5

备注

【数据范围】
40% 的数据满足：5<=n<=15
70% 的数据满足：5<=n<=80
100% 的数据满足：5<=n<=300，0<=s<=1000。边长度为不超过 1000 的正整数。

解析：

  暴力DFS。
  $N$只有$300$啊！先求出直径，然后枚举在直径上的两端点作为起点和终点，并把路径上的点标记，然后从路径上的每一个点DFS更新答案即可，时间复杂度$O(N^3)$

代码：

#include <bits/stdc++.h>
using namespace std;

const int Max=305;
int n,m,size,len,ans=1e9,maxx;
int f[Max][Max],dis[Max],first[Max],vis[Max],v[Max],tag[Max];
struct shu{int to,next,len;}edge[Max<<1];

inline int get_int()
{
	int x=0;char c;
	for(c=getchar();(!isdigit(c))&&(c!='-');c=getchar());
	for(;isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+c-'0';
	return x;
}

inline void build(int x,int y,int z)
{
	edge[++size].next=first[x],first[x]=size,edge[size].to=y,edge[size].len=z;
	edge[++size].next=first[y],first[y]=size,edge[size].to=x,edge[size].len=z;
}

inline void init()
{
	n=get_int(),m=get_int();
	memset(f,0x3f,sizeof(f));
	for(int i=1;i<n;i++)
	{
	  int x=get_int(),y=get_int(),z=get_int();
	  build(x,y,z),f[x][y]=f[y][x]=z;
	}
}
inline void pre()
{
	for(int k=1;k<=n;k++)
	  for(int i=1;i<=n;i++)
	    for(int j=1;j<=n;j++)
	      f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
	for(int i=1;i<=n;i++)
	{
	  f[i][i]=0;
	  for(int j=1;j<=n;j++) if(i!=j) len=max(len,f[i][j]);
	}
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++)
	    if(i!=j&&f[i][j]==len)
	    {
	      vis[i]=vis[j]=1;
	      for(int k=1;k<=n;k++)
	        if(j!=k&&i!=k&&f[i][k]+f[k][j]==len) vis[k]=1;
	    }
}
inline void dfs(int p,int d)
{
	tag[p]=1;
	maxx=max(maxx,d);
	for(int u=first[p];u;u=edge[u].next)
	{
	  int to=edge[u].to;
	  if(tag[to]||v[to]) continue;
	  tag[to]=1;
	  dfs(to,d+edge[u].len);
	}
}
inline void solve()
{
	for(int s=1;s<=n;s++)
	  for(int t=1;t<=n;t++)
	  if(vis[s]&&vis[t]&&f[s][t]<=m)
	  {
	  	maxx=0;
	 	memset(v,0,sizeof(v)),memset(tag,0,sizeof(tag));
	  	for(int k=1;k<=n;k++) if(f[s][k]+f[k][t]==f[s][t]&&vis[k]) v[k]=1;
	  	for(int k=1;k<=n;k++) if(v[k]) dfs(k,0);
	  	ans=min(ans,maxx);
	  }
	cout<<ans;
}

int main()
{
	init();
	pre();
	solve();
	return 0;
}