#include<bits/stdc++.h>
using namespace std;
ifstream f("lca.in");
ofstream g("lca.out");
int n, q, m;
int euler[100010], nivel[100011], lvl[100011], lst[100011];
int lg[100011], rmq[25][100011];
bool viz[100011];
vector<int>v[100001];
void dfs(int nod)
{
    m++;
    nivel[m]=lvl[nod];
    euler[m]=nod;
    for(int i=0; i<v[nod].size(); i++)
    {
        lvl[v[nod][i]]=lvl[nod]+1;
        dfs(v[nod][i]);
        m++;
        nivel[m]=lvl[nod];
        euler[m]=nod;
    }
    lst[nod]=m;
}
void build()
{
    for(int i=2; i<=100001; i++)
        lg[i]=lg[i/2]+1;
    for(int i=1; i<=m; i++)
        rmq[0][i]=i;
    for(int i=1; i<=lg[m]; i++)
        for(int j=1; j<=m; j++)
    {
        if(nivel[rmq[i-1][j]]<nivel[rmq[i-1][j+(1<<(i-1))]])
            rmq[i][j]=rmq[i-1][j];
        else
            rmq[i][j]=rmq[i-1][j+(1<<(i-1))];
    }
}
int LCA(int x, int y)
{
    int a=min(lst[x], lst[y]);
    int b=max(lst[x], lst[y]);
    int ln=lg[b-a+1];
    int ans=rmq[ln][a];
    if(nivel[rmq[ln][b-(1<<ln)+1]]<nivel[ans])
        ans=rmq[ln][b-(1<<ln)+1];
    return euler[ans];
}
int main()
{
    f>>n>>q;
    for(int i=2; i<=n; i++)
    {
        int x;
        f>>x;
        v[x].push_back(i);
    }
    dfs(1);
    build();
    for(; q; q--)
    {
        int a, b;
        f>>a>>b;
        g<<LCA(a, b)<<'\n';
    }
    return 0;
}