Pagini recente » Cod sursa (job #1398740) | Cod sursa (job #1200849) | Cod sursa (job #1361389) | Cod sursa (job #6160) | Cod sursa (job #2249569)
#include <fstream>
#include <vector>
#include <algorithm>
using namespace std;
const int NMAX = 1e5 + 5;
int n, q;
vector <int> graph[NMAX];
int euler[2 * NMAX], len, h[NMAX];
int pos[NMAX], lg[2 * NMAX];
int rmq[2 * NMAX][20];
void DFS(int node, int father)
{
h[node] = h[father] + 1;
euler[++len] = node;
pos[node] = len;
for (auto i : graph[node])
{
DFS(i, node);
euler[++len] = node;
}
}
void BuildRMQ()
{
for (int i = 2;i <= len;++i)
lg[i] = lg[i / 2] + 1;
for (int i = 1;i <= len;++i)
rmq[i][0] = i;
int x, y;
for (int p = 1;(1 << p) <= len;++p)
for (int i = 1;i + (1 << p) - 1 <= len;++i)
{
x = rmq[i][p - 1];
y = rmq[i + (1 << (p - 1))][p - 1];
rmq[i][p] = (h[euler[x]] < h[euler[y]] ? x : y);
}
}
inline int Query(int x, int y)
{
x = pos[x];
y = pos[y];
if (x > y)
swap(x, y);
int p = lg[y - x + 1];
int a = euler[rmq[x][p]];
int b = euler[rmq[y - (1 << p) + 1][p]];
return (h[a] < h[b] ? a : b);
}
int main()
{
ifstream fin("lca.in");
ofstream fout("lca.out");
fin >> n >> q;
int x, y, z;
for (int i = 2;i <= n;++i)
{
fin >> x;
graph[x].push_back(i);
}
DFS(1, 0);
BuildRMQ();
for (int i = 1;i <= q;++i)
{
fin >> x >> y;
fout << Query(x, y) << "\n";
}
fin.close();
fout.close();
return 0;
}