Pagini recente » Cod sursa (job #2803427) | Cod sursa (job #2179498) | Cod sursa (job #2750978) | Cod sursa (job #562066) | Cod sursa (job #3214113)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("lca.in");
ofstream fout("lca.out");
int n, m, Q;
int e[200005], niv[200005], expo[200005], P[100001];
int rmq[20][200005];
vector<int>L[100001];
void DFS(int k, int level)
{
m++;
niv[m] = level;
e[m] = k;
P[k] = m;
for(int i : L[k])
{
DFS(i, level + 1);
m++;
e[m] = k;
niv[m] = level;
}
}
void MakeE()
{
int i;
expo[0] = 1;
for(i = 2; i <= m; i++)
expo[i] = 1 + expo[i / 2];
}
int main()
{
ios_base::sync_with_stdio(0);
fin.tie(0);
fout.tie(0);
int i, j, x, y, l;
int A, B, rez;
fin >> n >> Q;
for(i = 2; i <= n; i++)
{
fin >> j;
L[j].push_back(i);
}
DFS(1, 1);
for(i = 1; i <= m; i++)
rmq[0][i] = i;
MakeE();
for(i = 1; i <= expo[m]; i++)
{
l = (1 << i);
for(j = 1; j <= m - l + 1; j++)
{
A = rmq[i - 1][j];
B = rmq[i - 1][j + l / 2];
if(niv[A] <= niv[B])
rmq[i][j] = A;
else rmq[i][j] = B;
}
}
for(i = 1; i <= Q; i++)
{
fin >> x >> y;
x = P[x];
y = P[y];
if(x > y)
swap(x, y);
int t = expo[y - x + 1];
l = (1 << t);
A = rmq[t][x];
B = rmq[t][y - l + 1];
if(niv[A] <= niv[B])
rez = A;
else rez = B;
fout << e[rez] << "\n";
}
return 0;
}