解法一(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
对于直径1,2,3……,t
我们选取i,j为路径
d[k]表示,点k不经过1-t点,可达的最远距离。
如下图:
显然:i到j的偏心距:
对于点i,j,不经过i,j上的点,且离其最远距离的点是1,t。
对于点i+1 -> j-1 中的点k,不经过i,j上的点,且离其最远距离的点是:d[k]。
综上 得:其偏心距等于:(根据前几篇博客的分析)
i,j可以尺取得到(见解法二尺取贪心)
前面的max(d[k]),用单调队列维护即可(单调队列维护区间[i-k,i],的最大值 ,你当然可以用RMQ代替,不过浪费空间没必要)
具体见代码
#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 = 5e5+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,kk=0;
int p,q,f[M],dp[M],dq[M],dx[M],Q[M];
int vs[M];
int n,s;vector<int>v;
void dfs(int x,int fa,int sm)
{
if(sm>mx)mx=sm,p=x;
if(q)f[x]=fa;
if(!kk)dq[x]=sm;else dp[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()
{
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);
q=p;
mx=0;dfs(p,0,0);
int z=p;
while(z)
{
v.pb(z);
z=f[z];
}
kk=1;dfs(p,0,0);
int L=v.size();
for(int i=0;i<L;i++)vs[v[i]]=1;
for(int i=0;i<L;i++)mx=0,dfs2(v[i],0,0),dx[i]=mx;
//cout<<"-----"<<v[i]<<" "<<dx[i]<<endl;;
Q[0]=-1;
int l=1,r=0,j=0,mi=2e9;
for(int i=0;i<L;i++)
{
//单调队列里放的是:直径上的点在vector里的位置。
while(dp[v[j]]-dp[v[i]]<=s&&j<L)
{
while(r>=l&&dx[j]>dx[Q[r]])//维护单调递增,因为范围最大,右边的大的更优
r--;
Q[++r]=j;
j++;
}
int tp=max(dp[v[i]],dq[v[j-1]]);
while(Q[l]<i)l++;
tp=max(tp,dx[Q[l]]);
mi=min(tp,mi);
//cout<<i<<" "<<j<<" "<<l<<" "<<r<<" "<<tp<<" = "<<Q[l]<<endl;
}
cout<<mi<<endl;
return 0;
}
解法五:
上面解法四还要用单调队列,不够优美。。
我们仔细思考max(d[k]) (i<=k<=j)这部分
可以直接转化为max(d[k])(1<=k<=t)即变成了一个常数。
为什么可以这样转化?
看解法4中的图:
1-i中d中的x,的d[x]:显然小于dist(1,i)//考虑直径的性质
j-t中的y,的d[x]:显然小于dist(j,t)
所以我们可以不用单调队列,只用尺取,省空间和码量!
#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 = 5e5+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,kk=0;
int p,q,f[M],dp[M],dq[M];
int vs[M];
int n,s;vector<int>v;
void dfs(int x,int fa,int sm)
{
if(sm>mx)mx=sm,p=x;
if(q)f[x]=fa;
if(!kk)dq[x]=sm;else dp[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()
{
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);
q=p;
mx=0;dfs(p,0,0);
int z=p;
while(z)
{
v.pb(z);
z=f[z];
}
kk=1;dfs(p,0,0);
int L=v.size();
for(int i=0;i<L;i++)vs[v[i]]=1;
int dx=0;
for(int i=0;i<L;i++)mx=0,dfs2(v[i],0,0),dx=max(dx,mx);
int j=0,mi=2e9;
for(int i=0;i<L;i++)
{
while(dp[v[j]]-dp[v[i]]<=s&&j<L)j++;
int tp=max(dp[v[i]],dq[v[j-1]]);
tp=max(tp,dx);
mi=min(tp,mi);
//cout<<i<<" "<<j<<" "<<l<<" "<<r<<" "<<tp<<" = "<<Q[l]<<endl;
}
cout<<mi<<endl;
return 0;
}