Pagini recente » Cod sursa (job #2499439) | Cod sursa (job #3126197) | Cod sursa (job #1178778) | Cod sursa (job #1794465) | Cod sursa (job #1363554)
#include<cstdio>
#include<string>
#include<vector>
#include<cmath>
using namespace std;
#ifdef HOME
const string inputFile = "input.txt";
const string outputFile = "output.txt";
#else
const string problemName = "lca";
const string inputFile = problemName + ".in";
const string outputFile = problemName + ".out";
#endif
typedef pair<int, int> PII;
const int NMAX = 100000 + 5;
int N, M;
vector<int> V[NMAX];
PII RMQ[19][2 * NMAX];
int euler[2 * NMAX], cnt;
int first[NMAX];
int lvl[NMAX];
void dfs(int x) {
euler[++cnt] = x;
first[x] = cnt;
for(auto y : V[x]) {
lvl[y] = lvl[x] + 1;
dfs(y);
euler[++cnt] = x;
}
}
void build() {
int i, j, p;
for(i = 1; i <= cnt; i++)
RMQ[0][i] = make_pair(lvl[euler[i]], euler[i]);
for(i = 1, p = 2; 1 + p / 2 <= cnt; i++, p *= 2)
for(j = 1; j + p / 2 <= cnt; j++)
RMQ[i][j] = min(RMQ[i - 1][j], RMQ[i - 1][j + p / 2]);
}
int query(int x, int y) {
int d = y - x + 1;
int i = log2(d + 0.5);
int p = (1 << i);
return min(RMQ[i][x], RMQ[i][y - p + 1]).second;
}
int lca(int x, int y) {
x = first[x];
y = first[y];
if(x > y)
swap(x, y);
return query(x, y);
}
int main() {
int i, x, y;
freopen(inputFile.c_str(), "r", stdin);
freopen(outputFile.c_str(), "w", stdout);
scanf("%d%d", &N, &M);
for(i = 2; i <= N; i++) {
scanf("%d", &x);
V[x].push_back(i);
}
dfs(1);
build();
while(M--) {
scanf("%d%d", &x, &y);
printf("%d\n", lca(x, y));
}
return 0;
}