codeforces 1062E Company dfs序+线段树+lca

思路代码来自：ccsu_cat

题目链接

题意：给你一棵树，让你删除区间【l,r】中的某个节点，使得他们的lca的深度最大，输出删除的节点及lca的深度。

思路：一段区间中对lca有影响的就是dfs序（深度优先搜素的顺序）最大的那个节点和dfs序最小的那个节点，那么假设输出最小的那个节点，那么此时的lca就是dfs序次小的与最大的的lca，假设删除最大的，那lca就是最小的和次大的，那么，我们比较一下lca，删除小的那个就可以了。

需要注意的是在预处理倍增父亲节点的时候，内层循环应该是n,外层应该是20，因为如果不这样处理，他的倍增父亲的值可能为0，有一些值没有存进去。

#include<bits/stdc++.h>
#define ok(x) cout<<#x<<'='<<x<<'\n';
using namespace std;
const int maxn=1e5+10;
vector<int>G[maxn];
int dep[maxn],f[maxn][20],id[maxn],cnt,s[maxn],n;
struct node{
	int mn,mx;
	node operator + (const node& t)const
	{
		node tmp;
		tmp.mn=min(mn,t.mn);
		tmp.mx=max(mx,t.mx);
		return tmp;
	}
}tr[maxn*4];
void dfs(int u,int fa,int deep)
{
	dep[u]=deep;
	f[u][0]=fa;
	s[u]=++cnt;	//dfs序
	id[cnt]=u;	//dfs序对应的节点的编号
	for(int i=0;i<G[u].size();i++)
		dfs(G[u][i],u,deep+1);
}

void init()
{
	for(int j=1;j<=18;j++)
	for(int i=1;i<=n;i++)
	{
		f[i][j]=f[f[i][j-1]][j-1];
	}
	
}
int Lca(int u,int v)
{
	if(dep[u]<dep[v])swap(u,v);
	for(int i=18;i>=0;i--)
	{
		if(dep[f[u][i]]>=dep[v])u=f[u][i];
	}
	if(u==v)return dep[u]-1;
	for(int i=18;i>=0;i--)
	{
		if(f[u][i]!=f[v][i])v=f[v][i],u=f[u][i]; 	
	}
	return dep[f[u][0]]-1;
}
#define ls o<<1
#define rs o<<1|1
#define mid (l+r)/2
void up(int o,int l,int r,int p,int v)
{
	if(l==r)
	{
		tr[o].mn=tr[o].mx=v;
		return ;
	}
	if(p<=mid)up(ls,l,mid,p,v);
	else up(rs,mid+1,r,p,v);
	tr[o]=tr[ls]+tr[rs];
}
node qu(int o,int l,int r,int L,int R)
{
	if(l>=L&&r<=R)return tr[o];
	node ans;
	ans.mn=1e9;
	ans.mx=0;
	if(L<=mid)ans=ans+qu(ls,l,mid,L,R);
	if(R>mid)ans=ans+qu(rs,mid+1,r,L,R);
	return ans;
}
int main()
{
	int q,l,r;
	scanf("%d%d",&n,&q);
	
	for(int i=2;i<=n;i++)
	{
		int x;
		scanf("%d",&x);
		G[x].push_back(i);
	}
	dfs(1,0,1);
	init();
	for(int i=1;i<=n;i++)
	{
		up(1,1,n,i,s[i]);
	}
	while(q--)
	{
		scanf("%d%d",&l,&r);
		node tmp=qu(1,1,n,l,r);
		int ti=tmp.mn;
		int tx=tmp.mx;
		
		up(1,1,n,id[ti],1e8);
		tmp=qu(1,1,n,l,r);
		int ii=tmp.mn;
		up(1,1,n,id[ti],ti);
		
		up(1,1,n,id[tx],0);
		tmp=qu(1,1,n,l,r);
		int xx=tmp.mx;
		up(1,1,n,id[tx],tx);	
	
		int t=Lca(id[ti],id[xx]);
		int t1=Lca(id[tx],id[ii]);
	
		if(t>=t1)printf("%d %d\n",id[tx],t);
		else printf("%d %d\n",id[ti],t1);
	}
	return 0;
}