Pagini recente » Cod sursa (job #1907309) | Cod sursa (job #2627199) | Cod sursa (job #637907) | Cod sursa (job #3210381) | Cod sursa (job #3248633)
#include <fstream>
#define maxlog 18
#define dim 100001
#include <vector>
using namespace std;
ifstream fin("lca.in");
ofstream fout("lca.out");
int pstart[dim], pend[dim], stramos[maxlog][dim], nivel[dim], pcrt, n;
vector <int> l[dim];
int tata[dim];
void dfs(int nod, int niv)
{
nivel[nod]=niv;
pstart[nod]=++pcrt;
for(int i=0;i<l[nod].size();i++)
{
dfs(l[nod][i], niv+1);
}
pend[nod]=++pcrt;
}
void calculeaza_stramos()
{
for(int i=1;i<=n;i++)
{
stramos[0][i]=tata[i];
}
for(int p=1;p<maxlog;p++)
{
for(int i=1;i<=n;i++)
{
stramos[p][i]=stramos[p-1][stramos[p-1][i]];
}
}
}
int este_stramos(int x, int y)
{
if(pstart[x]<=pstart[y]&&pend[y]<=pend[x])
{
return 1;
}
return 0;
}
int lca(int x, int y)
{
if(este_stramos(x, y))
{
return x;
}
if(este_stramos(y, x))
{
return y;
}
for(int p=maxlog-1;p>=0;p--)
{
int z=stramos[p][x];
if(z!=0&&este_stramos(z, y)==0)
{
x=z;
}
}
return stramos[0][x];
}
int main()
{
int m, i, x, y;
fin>>n>>m;
tata[1]=0;
for(i=2;i<=n;i++)
{
fin>>tata[i];
l[tata[i]].push_back(i);
}
pcrt=0;
dfs(1, 0);
calculeaza_stramos();
for(i=1;i<=m;i++)
{
fin>>x>>y;
fout<<lca(x, y)<<"\n";
}
}