2016vijos 6-1 松鼠聚会（LCA+卡空间)

#include<cstdio>

using namespace std;

#define N 1000001

const int D=576;

struct node
{
    char a,b,c;
}f[N],g[N];

inline int F(int i)
{
    return ((f[i].a<<7)+f[i].b<<7)+f[i].c; 
}

inline int G(int i)
{
    return ((g[i].a<<7)+g[i].b<<7)+g[i].c; 
}

inline int d(int i)
{
    return G(i)<N ? G(i):G(F(i))+1;
}

void read(int &x)
{
    x=0; char c=getchar();
    while(c<‘0‘||c>‘9‘) c=getchar();
    while(c>=‘0‘&&c<=‘9‘) { x=x*10+c-‘0‘; c=getchar(); }
}    

node turn(int x)
{
    node y;
    y.a=x>>14;
    y.b=(x>>7)&127;
    y.c=x&127;
    return y;
}

int lca(int x,int y)
{
    int u;
    if(d(x)<d(y)) u=x,x=y,y=u;
    int dis=d(x)-d(y);
    while(dis>=D)
    {
        if(G(x)<N) x=F(x),dis--;
        else x=G(x)-N,dis-=D;
    }
    while(dis) x=F(x),dis--;
    while(x!=y)
    {
        if(G(x)<N||G(x)==G(y)) x=F(x),y=F(y);
        else x=G(x)-N,y=G(y)-N;
    }
    return x; 
}

int main()
{
    freopen("squirrel.in","r",stdin);
    freopen("squirrel.out","w",stdout);
    int n,m;
    read(n); read(m);
    g[0]=turn(N);
    int x;
    for(int i=2;i<=n;++i) 
    {
        read(x);
        f[i]=turn(x),g[i]=turn(G(F(i))+1);
    }
    for(int i=n;i;i--)
        if(G(i)&1)
        {
            x=i;
            for(int j=1;j<=6;++j) x=F(x);
            g[i]=turn(x+N);
        }
    for(int i=n;i;i--)
        if(G(i)>N)
        {
            x=i;
            for(int j=1;j<=6;++j) x=G(x)-N;
            g[i]=turn(x+N);
        }
    for(int i=n;i;i--)
        if(G(i)>N)
        {
            x=i;
            for(int j=1;j<=4;++j) x=G(x)-N;
            g[i]=turn(x+N);
        }
    for(int i=n;i;i--)
        if(G(i)>N)
        {
            x=i;
            for(int j=1;j<=4;++j) x=G(x)-N;
            g[i]=turn(x+N);
        }
    int y;
    while(m--)
    {
        read(x); read(y);
        printf("%d\n",lca(x,y));
    }
}