codeforces 1062E

4 篇文章 0 订阅

传送门:https://codeforces.com/contest/1062/problem/E

题意:给你一棵树,每次询问[l,r]区间删掉一个点后,剩余的点的lca深度最大是多少,输出删除的点和lca深度。

题解:首先一堆点的lca就是dfs序最大的点和最小的点的lca,那么我们只用考虑删掉这两个点后的lca是多少就可以了,维护区间dfs序最大和次大,最小和次小。

#include<bits/stdc++.h>
using namespace std;
int head[100005],cnt,bi[100005],num;
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define drep(i,a,b) for(int i=a;i>=b;i--)
typedef pair<int,int> pi;
pi p[4];
#define fi first
#define se second
#define mp make_pair
struct nodes{
	pi zd,cd,zx,cx;
}tree[400005];
struct node{
	int to,nex;
}edge[200005];
void build(int t,int l,int r){
	if(l==r){
		tree[t].zd=tree[t].zx=mp(bi[l],l);
		tree[t].cx=mp(10000000,0);
		return ;
	}
	int mid=l+r>>1;
	build(t<<1,l,mid);
	build(t<<1|1,mid+1,r);
	p[0]=tree[t<<1].zd;
	p[1]=tree[t<<1].cd;
	p[2]=tree[t<<1|1].zd;
	p[3]=tree[t<<1|1].cd;
	sort(p,p+4);
	tree[t].zd=p[3];
	tree[t].cd=p[2];
	p[0]=tree[t<<1].zx;
	p[1]=tree[t<<1].cx;
	p[2]=tree[t<<1|1].zx;
	p[3]=tree[t<<1|1].cx;
	sort(p,p+4);
	tree[t].zx=p[0];
	tree[t].cx=p[1];
}
void add(int u,int v){
	edge[cnt].to=v;
	edge[cnt].nex=head[u];
	head[u]=cnt++;
}
pi nzd,ncd,nzx,ncx;
void query(int t,int l,int r,int x,int y){
	if(x<=l&&y>=r){
		p[0]=tree[t].zd;
		p[1]=tree[t].cd;
		p[2]=nzd;
		p[3]=ncd;
		sort(p,p+4);
		nzd=p[3];
		ncd=p[2];
		p[0]=tree[t].zx;
		p[1]=tree[t].cx;
		p[2]=nzx;
		p[3]=ncx;
		sort(p,p+4);
		nzx=p[0];
		ncx=p[1];
		return ;
	}
	int mid=l+r>>1;
	if(y<=mid)query(t<<1,l,mid,x,y);
	else if(x>mid)query(t<<1|1,mid+1,r,x,y);
	else{
		query(t<<1,l,mid,x,y);
		query(t<<1|1,mid+1,r,x,y);
	}
}
int f[100005][20],dep[100005];
int lca(int x,int y){
	if(x==0||y==0)return 0;
    if(dep[x]<dep[y])swap(x,y);
//令x为深度较深的点
    drep(i,16,0)if(dep[f[x][i]]>=dep[y])x=f[x][i];
//让x向上走到与y同一深度
    if(x==y)return x; //如果直接是lca,直接返回
    drep(i,16,0)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
//x,y同时向上走,直到父节点相同
    return f[x][0]; //返回父节点
}
void dfs(int t,int fa,int d){
	f[t][0]=fa;
	bi[t]=++num;
	dep[t]=d;
	for(int i=head[t];~i;i=edge[i].nex){
		int v=edge[i].to;
		if(v==fa)continue;
		dfs(v,t,d+1);
	}
}
#define fuck(x) cout<<#x<<' '<<x<<endl
int main(){
	int n,q;
	memset(head,-1,sizeof(head));
	scanf("%d%d",&n,&q);
	for(int i=2;i<=n;i++){
		int x;
		scanf("%d",&x);
		add(x,i);
		add(i,x);
	}
	dfs(1,1,1);
	build(1,1,n);
	rep(k,1,16) //假设10W个点
	rep(i,1,n)f[i][k]=f[f[i][k-1]][k-1]; 
	while(q--){
		int l,r;
		scanf("%d%d",&l,&r);
		nzd=ncd=mp(0,0);
		nzx=ncx=mp(n+1,n+1);
		query(1,1,n,l,r);
		int anslab=nzx.se,ansdep=dep[lca(nzd.se,ncx.se)];
//		fuck(ansdep);
//		fuck(dep[lca(ncd.se,nzx.se)]);
		if(dep[lca(ncd.se,nzx.se)]>ansdep){
			ansdep=dep[lca(ncd.se,nzx.se)];
			anslab=nzd.se;
		}
		printf("%d %d\n",anslab,ansdep-1);
	}
	return 0;
}
/*
11 5
1 1 3 3 3 4 2 7 7 6
4 8
1 11
9 11
8 11

*/

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值