Pagini recente » Cod sursa (job #3355347) | Statistici Borsos Zalan (borsoszalan) | Profil donea | Cod sursa (job #3353353) | Cod sursa (job #3353344)
#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;
cout<<euler[rmq(poz[x],poz[y])]<<endl;;
}
return 0;
}