Pagini recente » Cod sursa (job #1780415) | Cod sursa (job #576504) | Cod sursa (job #783980) | Cod sursa (job #755159) | Cod sursa (job #1696402)
#include <fstream>
#include <vector>
using namespace std;
const int NMAX=100005;
const int Roxi=2000005;
vector <int> G[NMAX];
int Euler[Roxi],First[NMAX],Rmq[Roxi][21],Lev[Roxi],Lg[Roxi],cnt;
void dfs(int x,int niv)
{
cnt++;
Euler[cnt]=x;
Lev[cnt]=niv;
First[x]=cnt;
for(int i=0;i<(int)G[x].size();i++)
{
int y=G[x][i];
dfs(y,niv+1);
cnt++;
Euler[cnt]=x;
Lev[cnt]=niv;
}
}
int main()
{
ifstream in("lca.in");
ofstream out("lca.out");
int n,m,i,j,x,y,k,l,lg,a,b;
in>>n>>m;
for(i=2;i<=n;i++)
{
in>>x;
G[x].push_back(i);
}
dfs(1,0);
for(i=2;i<=cnt;i++)
Lg[i]=Lg[i/2]+1;
for(i=1;i<=cnt;i++)
{
for(j=0;(1<<j)<=i;j++)
{
if(j==0)
{
Rmq[i][0]=i;
continue;
}
k=i-(1<<(j-1));
if(Lev[Rmq[i][j-1]]<=Lev[Rmq[k][j-1]]) Rmq[i][j]=Rmq[i][j-1];
else Rmq[i][j]=Rmq[k][j-1];
}
}
/*for(i=1;i<=cnt;i++)
{
for(j=0;(1<<j)<=i;j++)
{
printf("%d ",Rmq[i][j]);
}
printf("\n");
}*/
//for(i=1;i<=n;i++) printf("%d ",First[i]);
for(i=1;i<=m;i++)
{
in>>x>>y;
a=First[x];
b=First[y];
if(a>b) swap(a,b);
l=b-a+1;
lg=Lg[l];
//printf("%d %d\n",lg,l);
if(Lev[Rmq[b][lg]]<Lev[Rmq[a+(1<<lg)-1][lg]]) out<<Euler[Rmq[b][lg]]<<"\n";
else out<<Euler[Rmq[a+(1<<lg)-1][lg]]<<"\n";
}
return 0;
}