Pagini recente » Cod sursa (job #1446834) | Cod sursa (job #680574) | Cod sursa (job #2335454) | Cod sursa (job #915059) | Cod sursa (job #3002680)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("lca.in");
ofstream fout("lca.out");
int n, m, st[500005], nivel[1000005], poz[100005], p2[100005], p, rmq[20][100005];
vector<int>g[100005];
void Euler(int nod, int niv){
st[++p] = nod;
nivel[p] = niv;
poz[nod] = p;
for(auto i : g[nod]){
Euler(i, niv + 1);
st[++p] = nod;
nivel[p] = niv;
}
}
void Rmq(){
for(int i = 1; i <= p; i++)
rmq[0][i] = i;
for(int i = 1; (1 << i) <= p; i++)
for(int j = 1; j + (1 << i) - 1 <= p; j++){
if(nivel[rmq[i-1][j]] < nivel[rmq[i-1][j + (1 << (i-1))]]){
rmq[i][j] = rmq[i - 1][j];
}
else{
rmq[i][j] = rmq[i - 1][j + (1 << (i-1))];
}
}
}
int main()
{
int x, y;
fin >> n >> m;
for(int i = 2; i <= n; i++){
fin >> x;
g[x].push_back(i);
}
Euler(1, 1);
p2[1] = 0;
for (int i = 2; i <= p; i++)
p2[i] = 1 + p2[i / 2];
Rmq();
for(int i = 1; i <= m; i++){
fin >> x >> y;
x = poz[x];
y = poz[y];
if(x > y)
swap(x, y);
p = p2[y - x + 1];
if(nivel[rmq[p][x]] < nivel[rmq[p][y - (1 << p) + 1]])
fout << st[rmq[p][x]] << "\n";
else fout << st[rmq[p][y - (1 << p) + 1]] << "\n";
}
return 0;
}