Cod sursa(job #1235192)

Utilizator Vally77FMI Calinescu Valentin Gelu Vally77 Data 28 septembrie 2014 22:13:21
Problema Lowest Common Ancestor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.37 kb
#include <iostream>
#include <fstream>
#include <vector>
using namespace std;
ifstream ka("lca.in");
ofstream ki("lca.out");

const int N_MAX = 100005;
const int LOG_MAX = 20;
vector< vector<int> >graf(N_MAX + 1);

int poz[N_MAX + 1], grad[N_MAX + 1];
int repr[2 * N_MAX + 1], log[2 * N_MAX + 1];
int minim[2 * N_MAX + 1][LOG_MAX];

int n, m, k;

void rep_euler(int t)
{
    repr[++k] = t;
    poz[t] = k;
    for(int i = 0; i < graf[t].size(); i++)
    {
        grad[graf[t][i]] = grad[t] + 1;
        rep_euler(graf[t][i]);
        repr[++k]= t;
    }
}

int minul(int a, int b)
{
    if(grad[a] <= grad[b])
        return a;
    else
        return b;
}

int main()
{
    ka >> n >> m;
    int a;
    for(int i = 2; i <= n; i++)
    {
        ka >> a;
        graf[a].push_back(i);
    }
    rep_euler(1);
    minim[1][0] = 1;
    for(int i = 2; i <= k; i++)
    {
        log[i] = log[i >> 1] + 1;
        minim[i][0] = repr[i];
    }
    for(int i = 1; i <= log[k]; i++)
    {
        for(int j = 1; j <= k - (1 << i); j++)
        {
            minim[j][i] = minul(minim[j][i-1], minim[j + (1 << (i-1))][i-1]);
        }
    }
    int x, y;
    for(int i = 1; i <= m; i++)
    {
        ka >> x >> y;
        x = poz[x];
        y = poz[y];
        if(x > y)
            swap(x, y);
        int dif = log[y - x];
        ki << minul(minim[x][dif], minim[y - (1 << dif) +1][dif]) << '\n';
    }
}