Pagini recente » Cod sursa (job #2028967) | Cod sursa (job #390217) | Cod sursa (job #2731180) | Cod sursa (job #1286841) | Cod sursa (job #2569169)
#include <bits/stdc++.h>
#define dim 100005
using namespace std;
ifstream fin("lca.in");
ofstream fout("lca.out");
int prim[dim],nivel[dim],niv,euler[2*dim],k,Log[2*dim],d[30][2*dim],n,m,i,x,y;
vector<int> L[dim];
void e(int nod)
{
euler[++k]=nod;
prim[nod]=k;
//cout<<nod<<" "<<k<<endl;
nivel[nod]=niv;
for(auto it:L[nod])
{
niv++;
e(it);
niv--;
euler[++k]=nod;
}
}
int query(int st,int dr){
// cout<<st<<" "<<dr<<endl;;
int l=dr-st+1;
int lg=Log[l];
int p1=d[lg][st];
int minn=2*n,sol;
//cout<<st<<" "<<dr<<endl;
int p2=d[lg][dr-(1<<lg)];
if(nivel[euler[p1]]<nivel[euler[p2]])
return euler[p1];
else
return euler[p2];
}
int main()
{
fin>>n>>m;
for(i=2; i<=n; i++)
{
int x;
fin>>x;
L[x].push_back(i);
}
// niv=1;
e(1);
/// RMQ
for(i=1; i<=k; i++)
d[0][i]=euler[i];
// for(i=1;i<=k;i++)
// cout<<d[0][i]<<" ";
// cout<<endl;
int p=1;
for(int lun=1; lun<=22; lun++)
{
p*=2;
///
for(i=1; i<=k-p+1; i++)
{
/// inceput secv
int p1=d[lun-1][i];
int p2=d[lun-1][i+p];
if(nivel[p1]<nivel[p2])
{
d[lun][i]=p1;
}
else
d[lun][i]=p2;
}
}
Log[1]=0;
for(i=2; i<=k; i++)
Log[i]=1+Log[i/2];
// for(i=1;i<=k;i++)
// cout<<i<<" "<<Log[i]<<"\n";
// cout<<endl;
// for(i=1;i<=k;i++)
// cout<<euler[i]<<" ";
// cout<<endl;
// for(i=1;i<=k;i++)
// cout<<nivel[euler[i]]<<" ";
for(i=1;i<=m;i++){
fin>>x>>y;
x=prim[x];
y=prim[y];
if(x>y)
swap(x,y);
fout<<query(x,y)<<"\n";
}
}