Pagini recente » Cod sursa (job #2571009) | Cod sursa (job #670547) | Cod sursa (job #2325897) | Cod sursa (job #1672746) | Cod sursa (job #2776540)
#include <bits/stdc++.h>
using namespace std;
ifstream f ("lca.in");
ofstream g ("lca.out");
int n,q;
vector <int> v[100005];
int a[500005];
int lung;
int viz[500005];
int dp[500005];
int rmq[20][500005];
int prim[500005];
void dfs(int nod, int tata)
{
viz[nod]=1;
a[++lung]=nod;
dp[nod]=dp[tata]+1;
int ajuns=0;
for(int i=0; i<v[nod].size(); ++i)
{
int fiu=v[nod][i];
if(viz[fiu]==0)
{
dfs(fiu, nod);
ajuns=1;
a[++lung]=nod;
}
}
}
int main()
{
f>>n>>q;
for(int i=2; i<=n; ++i)
{
int tata;
f>>tata;
v[tata].push_back(i);
v[i].push_back(tata);
}
dfs(1, 0);
int exp=log2(lung)+1;
for(int j=0; j<=exp; ++j)
for(int i=lung+1 ; i<=pow(2, exp); ++i)
rmq[j][i]=9999999;
for(int i=1; i<=lung; ++i)
{
int nod=a[i];
rmq[0][i]=nod;
}
for(int put=1; put<=exp; ++put)
{
for(int i=1; i<=lung-pow(2,put); ++i)
{
int nod=a[i];
int x=pow(2, put-1);
int nod2=i+x;
if(dp[rmq[put-1][i]] < dp[rmq[put-1][nod2]])
rmq[put][i]=rmq[put-1][i];
else
rmq[put][i]=rmq[put-1][nod2];
}
}
for(int i=1; i<=lung; ++i)
{
if(prim[a[i]]==0)
prim[a[i]]=i;
}
for(int i=1; i<=q; ++i)
{
int x,y;
f>>x>>y;
int ii=prim[x];
int jj=prim[y];
if(ii>jj)
swap(ii, jj);
int put=log2(jj-ii+1);
int rez1=rmq[put][ii];
int dif=jj-pow(2,put)+1;
int rez2=rmq[put][dif];
if(dp[rez1]<dp[rez2])
g<<rez1<<"\n";
else
g<<rez2<<"\n";
}
return 0;
}