Cod sursa(job #3227911)

Utilizator gabriel.9619Gabriel Stefan Tita gabriel.9619 Data 3 mai 2024 18:51:58
Problema Lowest Common Ancestor Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.47 kb
#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";
    }
}