Pagini recente » Cod sursa (job #937574) | Cod sursa (job #2176574) | Cod sursa (job #1160932) | Cod sursa (job #444664) | Cod sursa (job #1160558)
#include <fstream>
#include <vector>
using namespace std;
#define NMax 100005
ifstream f("lca.in");
ofstream g("lca.out");
int n,m,k;
int gr[NMax];
int M[NMax][20],P[NMax],poz;
int put[20];
vector<int> v[NMax];
void DFS(int s)
{
M[++poz][0]=s;
P[s]=poz;
for(int i=0;i<(int)v[s].size();i++) gr[v[s][i]]=gr[s]+1,DFS(v[s][i]),M[++poz][0]=s;
}
void Preproc()
{
int i,j,nr=1;
while(nr<poz) put[k++]=nr,nr<<=1;
put[k]=nr;
int a,b;
for(j=1;j<=k;j++) for(i=1;i<poz;i++)
{
if(i+put[j-1]<=poz)
{
a=M[i][j-1];
b=M[i+put[j-1]][j-1];
if(gr[a]<=gr[b]) M[i][j]=a;
else M[i][j]=b;
}
else M[i][j]=M[i][j-1];
}
}
int bs(int x)
{
int step,j;
for(step=1;step<=k;step<<=1);
for(j=0;step;step>>=1)
if(j+step<=k && put[j+step]<=x) j+=step;
return j;
}
int main()
{
int i,x,y,p;
f>>n>>m;
for(i=2;i<=n;i++)
{
f>>x;
v[x].push_back(i);
}
DFS(1);
Preproc();
for(i=1;i<=m;i++);
for(i=1;i<=m;i++)
{
f>>x>>y;
x=P[x]; y=P[y];
p=bs(y-x+1);
if(gr[M[x][p]]<=gr[M[y-put[p]+1][p]]) g<<M[x][p]<<"\n";
else g<<M[y-put[p]+1][p]<<"\n";
}
f.close();
g.close();
return 0;
}