Cod sursa(job #1160558)

Utilizator ThomasFMI Suditu Thomas Thomas Data 30 martie 2014 17:13:09
Problema Lowest Common Ancestor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.27 kb
#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;
}