Cod sursa(job #2958992)

Utilizator Luka77Anastase Luca George Luka77 Data 29 decembrie 2022 16:29:44
Problema Lowest Common Ancestor Scor 40
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.44 kb
#include <bits/stdc++.h>
using namespace std;

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

const int NMAX = 1e5+5, LOG = 25;
int n, q, k;
int nodes[2*NMAX], depth[2*NMAX], poz[NMAX];
int rmq[LOG][NMAX];
vector<int>tree[NMAX];

inline void dfs(int node, int level)
{
    nodes[++k] = node;
    depth[k] = level;
    poz[node] = k;
    for(auto itr : tree[node])
    {
        dfs(itr, level +  1);
        nodes[++k] = node;
        depth[k] = level;
    }
}


inline void RangeMinimumQuery()
{
    for(int i = 1; i <= k; ++ i)
        rmq[0][i] = i;
    for(int i = 1; (1<<i) <= k; ++ i)
    {
        for(int j = 1; j <= k; ++ j)
        {
            if(depth[rmq[i-1][j]] < depth[rmq[i-1][j + (1<<(i-1))]])
                rmq[i][j] = rmq[i-1][j];
            else
                rmq[i][j] = rmq[i-1][j + (1 << (i-1))];
        }
    }
}

int main()
{
    fin >> n >> q;
    for(int i = 2; i <= n; ++ i)
    {
        int node;
        fin >> node;
        tree[node].push_back(i);
    }

    dfs(1, 0);
    RangeMinimumQuery();

    while(q--)
    {
        int a, b;
        fin >> a >> b;
        a = poz[a];
        b = poz[b];
        if(a > b)
            swap(a, b);
        int l = log2(b - a + 1);
        if(depth[rmq[l][a]] < depth[rmq[l][1 + b - (1<<l)]])
            fout << nodes[rmq[l][a]] << '\n';
        else
            fout << nodes[rmq[l][1 + b - (1<<l)]] << '\n';
    }
}