Cod sursa(job #1366015)

Utilizator CosminRusuCosmin Rusu CosminRusu Data 28 februarie 2015 17:48:50
Problema Lowest Common Ancestor Scor 20
Compilator cpp Status done
Runda Arhiva educationala Marime 1.53 kb
#include <iostream>
#include <vector>
#include <string.h>
#include <fstream>

using namespace std;

ifstream fin("lca.in");
ofstream fout("lca.out");

const int maxn = 100005;
const int maxlg = 20;

int n, m, first[maxn], euler[maxn << 1], rmq[maxlg][maxn << 1], k, lvl[maxn << 1], lg[maxn << 1];
vector <int> g[maxn];

inline void dfs(int node, int height) {
    first[node] = ++ k;
    lvl[k] = height;
    euler[k] = node;
    for(vector <int> :: iterator it = g[node].begin() ; it != g[node].end() ; ++ it) {
        dfs(*it, height + 1);
        euler[++ k] = node;
        lvl[k] = height;
    }
}

inline int query(int x, int y) {
    x = first[x];
    y = first[y];
    if(x > y)
        swap(x, y);
    int l = lg[y - x + 1];
    int ret = rmq[l][x];
    if(lvl[ret] > lvl[rmq[l][y - (1 << l) + 1]])
        ret = rmq[l][y - (1 << l) + 1];
    return euler[ret];
}

int main() {
    fin >> n >> m;
    for(int i = 2 ; i <= n ; ++ i) {
        int x;
        fin >> x;
        g[x].push_back(i);
    }
    dfs(1, 1);
    for(int i = 1 ; i <= k ; ++ i)
        rmq[0][i] = i;
    for(int i = 2 ; i <= n ; ++ i)
        lg[i] = lg[i >> 1] + 1;
    for(int i = 1 ; (1 << i) <= k ; ++ i) {
        int aux = (1 << (i - 1));
        for(int j = 1 ; j + (1 << i) - 1 <= k ; ++ j) {
            rmq[i][j] = rmq[i - 1][j];
            if(lvl[rmq[i][j]] > lvl[rmq[i - 1][j + aux]])
                rmq[i][j] = rmq[i - 1][j + aux];
        }
    }
    for(int i = 1 ; i <= m ; ++ i) {
        int x, y;
        fin >> x >> y;
        fout << query(x, y) << '\n';
    }
}