P1099 树网的核 （证明为何选一条直径即可） 一道很好的树的直径性质题 4种解法 这一篇是O(n^3) 性质分析+暴力

解法一（On^3）：https://blog.csdn.net/bjfu170203101/article/details/106337107
解法二（On^2）：https://blog.csdn.net/bjfu170203101/article/details/106338373
解法三（Onlogm）：https://blog.csdn.net/bjfu170203101/article/details/106339621
解法四/五（On）：https://blog.csdn.net/bjfu170203101/article/details/106342984

首先题目定义了点到路径的距离为：

点到该路径上任意点最近的距离。

然后给出偏心距的定义：一段路径F的偏心距为：其他点到这条路径距离的最大值。

给出s，找出一条路径（必须在直径上）长度小于等于s的路径T，使其偏心距最小。

题目要求必须在直径上，其实仔细分析发现，一定是在直径上，因为如果有一段不在直径上，离这一段路径最远的点，一定是直径的某个（不在路径上的）端点，但这段距离显然大于端点直接到在直径上路径的距离，不优，白白浪费s。

 

找一条直径很好找，两边DFS即可（不清楚的可以看我之前blog的分析）。但如果有多个直径怎么办？

显然：多个直径的偏心距相同，且一定等于直径中心（点或路径）到直径两端的距离。

灵魂画师启动！

如上图：

S为直径中点（如果直径中点是在路径上直接缩成点一样分析），AB,CD是两条直径。

由于要选取一段路径T。

1.如果T不经过S，假设在直径AB上，那么偏心距一定大于max(CS,SD)（因为还要加上S到路径T的距离）

路径在CD上同理。

也就是说，如果路径不经过S，偏心距始终大于直径中心S到直径两端的距离。

 

2.如果T经过，假设S无限大，取AB为路径。

显然其他点到S最大距离为max（CS,SD）。即偏心距等于：直径中心S到直径两端的距离。

 

由1.2有：我们只需要在一条直径上找路径T，即可求出唯一的最小偏心距。

（如果多条直径，偏心距固定为直径中心S到直径两端的距离，否则为枚举的最小值）

 

然后就有了n^3暴力的做法：

枚举直径上路径长度小于s的路径。标记当前路径的点。

然后从这些点出发dfs。能得到这些点到其他每个点的最小距离（因为树上任意两点的距离唯一，路径唯一）

然后这题就能AC小数据了。 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
const double PI= acos(-1.0);
const int M = 300+7;

int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt,w;}ee[M*2];
void add(int x,int y,int w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}

int mx=0;
int p,q,f[M];
int d[M][M],rt,vs[M];
void dfs(int x,int fa,int sm)
{
	if(sm>mx)mx=sm,p=x;
	if(q)f[x]=fa;
	d[rt][x]=sm;
	for(int i=head[x];i;i=ee[i].nxt)
	{
		int y=ee[i].to,w=ee[i].w;
		if(y==fa)continue;
		dfs(y,x,sm+w);
	}
}
bool flag;
void dfs2(int x,int fa,int sm)
{
	if(flag)return ;
	mx=max(sm,mx);//树上任意两点距离唯一，走非标记点所到达的点，一定是最近的 
	for(int i=head[x];i;i=ee[i].nxt)
	{
		int y=ee[i].to,w=ee[i].w;
		if(y==fa||vs[y])continue;
		dfs2(y,x,sm+w);
	}
}
int main()
{
	int n,s;
	scanf("%d%d",&n,&s);
	for(int i=1;i<n;i++)
	{
		int u,v,w;
		scanf("%d%d%d",&u,&v,&w);
		add(u,v,w);add(v,u,w);
	}
	dfs(1,0,0);
	vector<int>v;q=p;
	mx=0;dfs(p,0,0);
	while(p)
	{
		v.pb(p);
		p=f[p];
	}
	for(int i=1;i<=n;i++)
	{
		rt=i;
		dfs(i,0,0);
	}
	
//	cout<<"okk"<<"   "<<p<<" "<<q<<endl;
	
	int mi=1e9;
	for(int i=0;i<v.size();i++)
	for(int j=i;j<v.size();j++)
	{
		int x=v[i],y=v[j];
//		cout<<" -  "<<x<<" "<<y<<" "<<d[x][y]<<endl;
		if(d[x][y]>s)continue;
		for(int k=i;k<=j;k++)vs[v[k]]=1;
		mx=0;
		for(int k=i;k<=j;k++)dfs2(v[k],0,0);
		for(int k=i;k<=j;k++)vs[v[k]]=0;
	//	cout<<tp<<endl;
		mi=min(mi,mx);
	}
	printf("%d\n",mi);
	return 0;
}