Pagini recente » Cod sursa (job #836606) | Cod sursa (job #885075) | Cod sursa (job #692978) | Cod sursa (job #855450) | Cod sursa (job #473842)
Cod sursa(job #473842)
#include<fstream>
#include<vector>
using namespace std;
#define NMAX 100003
ifstream fi("lca.in");
vector<int> f[NMAX],euler;
int poz[NMAX],niv[NMAX],rmq[20][NMAX<<1],log[NMAX<<1],d[20];
int n,m,x,i,a,b;
void Euler(int nod, int nv) {
int i,s;
poz[nod]=euler.size();
euler.push_back(nod);
niv[nod]=nv;
for(i=0,s=f[nod].size();i<s;i++) {
Euler(f[nod][i],nv+1);
euler.push_back(nod);
}
}
void calcD() {
d[0]=1;
for(int i=1;i<20;i++) d[i]=d[i-1]<<1;
}
void calcLog() {
int i,j=1,s=euler.size();
for(i=1;i<s;i++) {
if(i==d[j]) j++;
log[i]=j-1;
}
}
void preRMQ() {
int i,j,s;
for(i=0,s=euler.size();i<s;i++) rmq[0][i]=euler[i];
for(i=1;d[i]<=s;i++)
for(j=0;j<=s-d[i];j++)
if(niv[rmq[i-1][j]]<niv[rmq[i-1][j+d[i-1]]])
rmq[i][j]=rmq[i-1][j];
else rmq[i][j]=rmq[i-1][j+d[i-1]];
}
int RMQ(int a, int b) {
int best=rmq[log[b-a+1]][a];
a+=d[log[b-a+1]];
while(a<=b) {
if(niv[rmq[log[b-a+1]][a]]<niv[best]) best=rmq[log[b-a+1]][a];
a+=d[log[b-a+1]];
}
return best;
}
int main() {
fi>>n>>m;
for(i=2;i<=n;i++) {
fi>>x;
f[x].push_back(i);
}
Euler(1,1);
calcD();
calcLog();
preRMQ();
freopen("lca.out","w",stdout);
while(m--) {
fi>>a>>b;
if(poz[a]>poz[b]) swap(a,b);
printf("%d\n",RMQ(poz[a],poz[b]));
}
return 0;
}