Cod sursa(job #1363554)

Utilizator mirceadinoMircea Popoveniuc mirceadino Data 27 februarie 2015 01:11:01
Problema Lowest Common Ancestor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.55 kb
#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;
}