Pagini recente » Cod sursa (job #2011789) | Cod sursa (job #2010486) | Cod sursa (job #2499846) | Cod sursa (job #2004084) | Cod sursa (job #2868800)
#include <bits/stdc++.h>
using namespace std;
ifstream f("lca.in");
ofstream g("lca.out");
int lvl[400005],d[400005][20],n,m,k,euler[400005],x,poz[400005],y,mn,lg;
vector<int>v[100005];
void dfs(int nod,int nivel)//parcurgerea euler
{
k++;
lvl[k]=nivel;
euler[k]=nod;
poz[nod]=k;
for(int i=0;i<v[nod].size();i++)
{
dfs(v[nod][i],nivel+1);
k++;
lvl[k]=nivel;
euler[k]=nod;
}
}
int main()
{
int i,j;
f>>n>>m;
for(i=2; i<=n; i++)
{
f>>x;
v[x].push_back(i);
}
dfs(1,0);
//rmq;
for(i=0;i<=k;i++) d[i][0]=i;
for(j=1;(1<<j)<=k;j++)
{
for(i=1;i<=k;i++)
{
if(lvl[d[i][j-1]]<lvl[d[i+(1<<(j-1))][j-1]]) d[i][j]=d[i][j-1];
else d[i][j]=d[i+(1<<(j-1))][j-1];
}
}
for(i=1; i<=m; i++)
{
//lca
f>>x>>y;
x=poz[x];
y=poz[y];
if(x>y) swap(x,y);
lg=log2(y-x+1);
if(lvl[d[x][lg]]<lvl[d[y-(1<<lg)+1][lg]]) mn=d[x][lg];
else mn=d[y-(1<<lg)+1][lg];
g<<euler[mn]<<'\n';
}
return 0;
}