Cod sursa(job #2231613)

Utilizator DragosBledeaDragos Bledea DragosBledea Data 15 august 2018 10:34:15
Problema Lowest Common Ancestor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.48 kb
#include <fstream>
#include <vector>
#include <algorithm>
using namespace std;

int N,Q;
ifstream fi("lca.in");
ofstream fo("lca.out");
vector <int> D[100001];
int E[200001],k;
int AD[100001];
int PAP[100001];
int T[200001][18],nrc,p;

void euler(int v,int ad)
{
	vector <int> :: iterator it;
	E[++k]=v;
	AD[v]=ad;
	PAP[v]=k;
	for (it=D[v].begin();it!=D[v].end();it++)
	{
		euler(*it,ad+1);
		E[++k]=v;
	}
}

int query(int st, int dr)
/// returneaza varful cu adancime minima din portiunea de indici st->dr
{
	int vfrez;
	vfrez=E[st];
	for (int p=nrc;p>=0;p--)
		if (st+(1<<p)-1<=dr)
		{
			if (AD[T[st][p]]<AD[vfrez])
				vfrez=T[st][p];
			st=st+(1<<p);
		}
	return vfrez;
}

int main()
{
	fi>>N>>Q;
	for (int i=2;i<=N;i++)
	{
		int v;
		fi>>v;
		// v este parintele lui i
		D[v].push_back(i);
	}
	/// se obtine traversarea Euler a arborelui
	k=0;
	euler(1,0);
	/// se construieste T
	nrc=0;
	p=1;
	while ((p<<1)<=k)
	{
		p=(p<<1);
		nrc++;
	}
	/// coloana 0
	for (int i=1;i<=k;i++)
		T[i][0]=E[i];
	/// alte coloane
	for (int j=1;j<=nrc;j++)
		for (int i=1;i<=k-(1<<j)+1;i++)
			if (AD[T[i][j-1]]<=AD[T[i+(1<<(j-1))][j-1]])
				T[i][j]=T[i][j-1];
			else
				T[i][j]=T[i+(1<<(j-1))][j-1];
	/// se citesc interogarile si se raspunde la ele
	for (int i=1;i<=Q;i++)
	{
		int a,b;
		fi>>a>>b;
		int st,dr;
		st=PAP[a];
		dr=PAP[b];
		if (st>dr)
			swap(st,dr);
		fo<<query(st,dr)<<"\n";
	}
	fi.close();
	fo.close();
	return 0;
}