Cod sursa(job #3353345)

Utilizator Mirc100Mircea Octavian Mirc100 Data 6 mai 2026 11:56:47
Problema Lowest Common Ancestor Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.43 kb
#include <iostream>
#include <vector>
#include <fstream>
#include<cmath>
using namespace std;

vector<int> nivel;
vector<int> euler;
vector<int> fii[100001];
int mat[200001][20];
int poz[100001];
ifstream fin("lca.in");
ofstream fout("lca.out");

void parc(int x, int niv){
    euler.push_back(x);
    nivel.push_back(niv);
    poz[x]=euler.size()-1;
    for(int y : fii[x])
    {
        parc(y,niv+1);
        euler.push_back(x);
        nivel.push_back(niv);
    }
}
int rmq(int i, int j) {
    if(i>j) swap(i,j);
    int k = log2(j - i + 1);
    return nivel[mat[i][k]]<nivel[ mat[j - (1 << k) + 1][k]] ?mat[i][k]:mat[j - (1 << k) + 1][k];
}
void calc_rmq(int n){
    for(int i = 0; i < n; i ++) {
        mat[i][0]=i;
    }
    for(int j = 1; j <= log2(n); j++) {
        for(int i = 0; i + (1 << j) <= n; i++) {
                if(nivel[mat[i][j-1]]<nivel[mat[i + (1 << (j-1))][j-1]])
                    mat[i][j]=mat[i][j-1];
                else
                    mat[i][j]=mat[i + (1 << (j-1))][j-1];
           // mat[i][j] = min(mat[i][j-1], mat[i + (1 << (j-1))][j-1]);
        }
    }
}
int main(){
    int n,m,x,y;
    fin>>n>>m;
    for(int i=2;i<=n;i++)
       {
        fin>>x;
        fii[x].push_back(i);
       }
    parc(1,0);
      calc_rmq(nivel.size());
    for(int i=0;i<m;i++){
        fin>>x>>y;
        fout<<euler[rmq(poz[x],poz[y])]<<endl;;

    }

    return 0;
}