Pagini recente » Cod sursa (job #91033) | Cod sursa (job #2662548) | Cod sursa (job #694279) | Cod sursa (job #743285) | Cod sursa (job #1888806)
#include <bits/stdc++.h>
#define LGMax 21
#define NMax 100005
using namespace std;
ifstream fin("lca.in");
ofstream fout("lca.out");
vector<int> Graf[NMax];
int Euler[4*NMax],Level[4*NMax],First[4*NMax],lg[4*NMax],x,N,M,k;
int RMQ[LGMax][4*NMax];
void Read()
{
fin>>N>>M;
for(int i=2;i<=N;i++)
{
fin>>x;
Graf[x].push_back(i);
}
}
void DFS(int nod,int level)
{
Euler[++k]=nod;
Level[k]=level;
First[nod]=k;
for(vector<int>::iterator it=Graf[nod].begin();it!=Graf[nod].end();it++)
{
DFS(*it,level+1);
Euler[++k]=nod;
Level[k]=level;
}
}
void rmq()
{
for(int i=2;i<=k;i++)
lg[i]=lg[i/2]+1;
for(int i=1;i<=k;i++)
RMQ[0][i]=i;
for(int i=1;(1<<i)<=k;i++)
for(int j=1;j<=k-(1<<i);j++)
{
int l=1<<(i-1);
RMQ[i][j]=RMQ[i-1][j];
if(Level[RMQ[i][j]]>Level[RMQ[i-1][j+l]])
RMQ[i][j]=RMQ[i-1][j+l];
}
}
int LCA(int x,int y)
{
x=First[x];
y=First[y];
if(x>y)
swap(x,y);
int l=lg[y-x+1];
int sol=RMQ[l][x];
if(Level[sol]>Level[RMQ[l][y-(1<<l)+1]])
sol=RMQ[l][y-(1<<l)+1];
return Euler[sol];
}
int main()
{
Read();
DFS(1,0);
rmq();
for(int i=1;i<=M;i++)
{
int x,y;
fin>>x>>y;
fout<<LCA(x,y)<<"\n";
}
return 0;
}