Pagini recente » Cod sursa (job #2137144) | Cod sursa (job #1383056) | Cod sursa (job #2296738) | Cod sursa (job #1571480) | Cod sursa (job #574757)
Cod sursa(job #574757)
#include <cstdio>
#include <bitset>
#include <vector>
#include <cstring>
#include <string>
#define maxn 100010
using namespace std;
bitset <maxn> viz;
int E[maxn << 1], L[maxn << 1], app[maxn];
int help[20][maxn << 1], rmq[20][maxn << 1], lg[maxn << 1];
vector <int> G[maxn];
int N, Q, p;
void DFS (int node, int lev)
{
viz[node] = 1;
E[++p] = node;
L[p] = lev;
app[node] = p;
for (vector <int> :: iterator it = G[node].begin (); it != G[node].end (); it++)
if (viz[*it] == 0) {
DFS (*it, lev + 1);
E[++p] = node;
L[p] = lev;
}
}
int main ()
{
freopen ("lca.in", "r", stdin);
freopen ("lca.out", "w", stdout);
scanf ("%d %d\n", &N, &Q);
int i, x, j, pw;
for (i = 2; i <= N; i++) {
scanf ("%d", &x);
G[x].push_back (i);
}
DFS (1, 0);
/* for (i = 1; i <= p; i++)
printf ("%d ", E[i]);
printf ("\n");
for (i = 1; i <= p; i++)
printf ("%d ", L[i]);
*/
for (i = 1; i <= p; i++) {
rmq[0][i] = L[i];
help[0][i] = E[i];
}
for (j = 1; (1 << j) <= p; j++)
for (i = 1, pw = 1 << (j - 1); i + pw <= p; i++)
if (rmq[j - 1][i] < rmq[j - 1][i + pw]) {
rmq[j][i] = rmq[j - 1][i];
help[j][i] = help[j - 1][i];
} else {
rmq[j][i] = rmq[j - 1][i + pw];
help[j][i] = help[j - 1][i + pw];
}
for (i = 2; i <= p; i++)
lg[i] = lg[i >> 1] + 1;
int y, lgg;
for (; Q != 0; Q--) {
scanf ("%d %d\n", &x, &y);
if (app[x] > app[y]) swap (x, y);
x = app[x];
y = app[y];
lgg = lg[y - x + 1];
// printf ("%d %d %d\n", x, y, lgg);
if (rmq[lgg][x] > rmq[lgg][x + (y - x + 1) - (1 << lgg)])
printf ("%d\n", help[lgg][x + (y - x + 1) - (1 << lgg)]);
else printf ("%d\n", help[lgg][x]);
}
return 0;
}