算法竞赛——进阶指南——acwing351. 树网的核 解法四: 单调队列优化到 On + 解法五:利用性质遍历On

解法一(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]。

综上 得:其偏心距等于:(根据前几篇博客的分析)

max(max_{i<=k<=j}(d[k]),dist(1,i),dist(j,t))

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

 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值