Cod sursa(job #2776540)

Utilizator marcumihaiMarcu Mihai marcumihai Data 20 septembrie 2021 09:27:14
Problema Lowest Common Ancestor Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.86 kb
#include <bits/stdc++.h>

using namespace std;

ifstream f ("lca.in");
ofstream g ("lca.out");

int n,q;
vector <int> v[100005];
int a[500005];
int lung;
int viz[500005];
int dp[500005];
int rmq[20][500005];
int prim[500005];



void dfs(int nod, int tata)
{
    viz[nod]=1;
    a[++lung]=nod;
    dp[nod]=dp[tata]+1;
    int ajuns=0;
    for(int i=0; i<v[nod].size(); ++i)
    {
        int fiu=v[nod][i];
        if(viz[fiu]==0)
        {
            dfs(fiu, nod);
            ajuns=1;
            a[++lung]=nod;
        }

    }



}

int main()
{
    f>>n>>q;
    for(int i=2; i<=n; ++i)
    {
        int tata;
        f>>tata;
        v[tata].push_back(i);
        v[i].push_back(tata);
    }
    dfs(1, 0);

    int exp=log2(lung)+1;

    for(int j=0; j<=exp; ++j)
        for(int i=lung+1 ; i<=pow(2, exp); ++i)
            rmq[j][i]=9999999;
    for(int i=1; i<=lung; ++i)
    {
        int nod=a[i];
        rmq[0][i]=nod;
    }

    for(int put=1; put<=exp; ++put)
    {
        for(int i=1; i<=lung-pow(2,put); ++i)
        {
            int nod=a[i];
            int x=pow(2, put-1);
            int nod2=i+x;
            if(dp[rmq[put-1][i]] < dp[rmq[put-1][nod2]])
                rmq[put][i]=rmq[put-1][i];
            else
                rmq[put][i]=rmq[put-1][nod2];

        }
    }

    for(int i=1; i<=lung; ++i)
    {
        if(prim[a[i]]==0)
            prim[a[i]]=i;
    }
    for(int i=1; i<=q; ++i)
    {
        int x,y;
        f>>x>>y;
        int ii=prim[x];
        int jj=prim[y];
        if(ii>jj)
            swap(ii, jj);
        int put=log2(jj-ii+1);
        int rez1=rmq[put][ii];
        int dif=jj-pow(2,put)+1;
        int rez2=rmq[put][dif];
        if(dp[rez1]<dp[rez2])
            g<<rez1<<"\n";
        else
            g<<rez2<<"\n";
    }
    return 0;
}