#include <bits/stdc++.h>

using namespace std;

const int NMAX = 1e5, LG = 17;

int n, q, niv[NMAX + 1], first[NMAX + 1], etl,
    sparse[2 * NMAX + 1][LG + 1], lg2[2 * NMAX + 1];
vector<int> adj[NMAX + 1], euler;

void precalcLg2() {
    lg2[1] = 0;
    for(int i = 2; i <= 2 * NMAX; ++i)
        lg2[i] = lg2[i / 2] + 1;
}

void eulerTour(int nod = 1, int par = -1, int nivCrt = 0) {
    niv[nod] = nivCrt;
    first[nod] = euler.size();
    euler.push_back(nod);
    for(const auto &el : adj[nod]) {
        if(el == par) continue;
        eulerTour(el, nod, nivCrt + 1);
        euler.push_back(nod);
    }
}

void constructSparse() {
    for(int i = 1; i <= etl; ++i)
        sparse[i][0] = euler[i];
    for(int sz = 1; (1 << sz) <= etl; ++sz)
        for(int i = 1; i + (1 << sz) - 1 <= etl; ++i) {
            const int v1 = sparse[i][sz - 1],
                      v2 = sparse[i + (1 << (sz - 1))][sz - 1];
            if(niv[v1] < niv[v2]) sparse[i][sz] = v1;
            else sparse[i][sz] = v2;
        }
}

int lca(int x, int y) {
    x = first[x];
    y = first[y];
    if(x > y) swap(x, y);
    int lg2dist = lg2[y - x + 1];
    const int v1 = sparse[x][lg2dist],
              v2 = sparse[y - (1 << lg2dist) + 1][lg2dist];
    if(niv[v1] < niv[v2]) return v1;
    return v2;
}

int main()
{
    precalcLg2();
    freopen("lca.in", "r", stdin);
    freopen("lca.out", "w", stdout);
    scanf("%d%d", &n, &q);
    for(int i = 2, x; i <= n; ++i) {
        scanf("%d", &x);
        adj[x].push_back(i);
        adj[i].push_back(x);
    }

    euler.push_back(0);
    eulerTour();
    etl = euler.size() - 1;
    constructSparse();
    for(int i = 1, x, y; i <= q; ++i) {
        scanf("%d%d", &x, &y);
        printf("%d\n", lca(x, y));
    }
    return 0;
}