首先可以求出从某点做$2^k$次车能到的最浅的点，这个只要dfs一下，把它的孩子能到的最浅的点更新过来就可以
然后倍增地往上跳，不能跳到lca的上面，记录坐车的次数ans
此时有三种情况（设最远能跳到x,y）：
1.再跳也跳不到lca的上面，就是-1
2.路径(x,y)被某趟车覆盖，答案是ans+1
3.并没有被覆盖，答案是ans+2
那么怎么看有没有覆盖呢
首先，如果这两个点是直上直下的（有一个是lca），只要看下面那个点能不能跳到上面去就行
对于剩下的(x,y)，只要有车的两端点分别在x和y的子树中就可以
所以做一遍dfs，进某个点x的时候记下来y的子树中车端点的个数，然后把从x发的车的终点++，出来的时候再次统计那个个数，如果不相同，就说明(x,y)这条路径被覆盖了。这个在dfs序上用一个树状数组就可以
不要把N写成M然后WA一页，很难看。
#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=2e5+; inline ll rd(){
ll x=;char c=getchar();int neg=;
while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
while(c>=''&&c<='') x=x*+c-'',c=getchar();
return x*neg;
} int eg[maxn][],egh[maxn],ect;
int N,M,Q,dep[maxn],fa[maxn][],nxt[maxn][];
int dfn[maxn][],tot,ans[maxn];
int pt[maxn*][],pth[maxn],tr[maxn];
int lq[maxn*][],lqh[maxn],tmp[maxn*];
bool covered[maxn]; inline void adeg(int a,int b){
eg[++ect][]=b,eg[ect][]=egh[a],egh[a]=ect;
} void dfs(int x){
dfn[x][]=++tot;
for(int i=;fa[x][i]&&fa[fa[x][i]][i];i++)
fa[x][i+]=fa[fa[x][i]][i];
for(int i=egh[x];i;i=eg[i][]){
int b=eg[i][];
dep[b]=dep[x]+;
dfs(b);
}dfn[x][]=tot;
} int getlca(int x,int y){
if(dep[x]<dep[y]) swap(x,y);
for(int i=;i>=&&dep[x]!=dep[y];i--){
if(fa[x][i]&&dep[fa[x][i]]>=dep[y])
x=fa[x][i];
}
if(x==y) return x;
for(int i=;i>=;i--){
if(fa[x][i]!=fa[y][i])
x=fa[x][i],y=fa[y][i];
}
return fa[x][];
} void dfs2(int x){
for(int i=egh[x];i;i=eg[i][]){
int b=eg[i][];
dfs2(b);
if(nxt[b][]&&(dep[nxt[b][]]<dep[nxt[x][]]||!nxt[x][]))
nxt[x][]=nxt[b][];
}
if(dep[nxt[x][]]>=dep[x]) nxt[x][]=;
} inline int lowbit(int x){return x&(-x);}
inline void add(int x,int y){
for(;x<=N;x+=lowbit(x)) tr[x]+=y;
}
inline int query(int x){
int re=;for(;x;x-=lowbit(x)) re+=tr[x];return re;
} void dfs3(int x){
for(int i=lqh[x];i;i=lq[i][]){
int b=lq[i][];
tmp[i]=query(dfn[b][])-query(dfn[b][]-);
}
for(int i=pth[x];i;i=pt[i][]){
int b=pt[i][];
add(dfn[b][],);
}
for(int i=egh[x];i;i=eg[i][]){
int b=eg[i][];
dfs3(b);
}
for(int i=lqh[x];i;i=lq[i][]){
int b=lq[i][];
if(tmp[i]!=query(dfn[b][])-query(dfn[b][]-))
covered[i>>]=;
}
} int main(){
// freopen("983E.in","r",stdin);
int i,j,k;
N=rd();
for(i=;i<=N;i++){
fa[i][]=rd();
adeg(fa[i][],i);
}
dep[]=;dfs();
M=rd();
for(i=;i<=M;i++){
int a=rd(),b=rd();
pt[i<<][]=a,pt[i<<][]=pth[b],pth[b]=i<<;
pt[i<<|][]=b,pt[i<<|][]=pth[a],pth[a]=i<<|;
int x=getlca(a,b);
if(dep[x]<dep[nxt[a][]]||!nxt[a][]) nxt[a][]=x;
if(dep[x]<dep[nxt[b][]]||!nxt[b][]) nxt[b][]=x;
}
dfs2();
for(i=;i<;i++){
for(j=;j<=N;j++){
// if(!nxt[j][i]||!nxt[nxt[j][i]][i]) continue;
nxt[j][i+]=nxt[nxt[j][i]][i];
}
}
Q=rd();
for(i=;i<=Q;i++){
int a=rd(),b=rd(),lca=getlca(a,b);
for(j=;j>=;j--){
if(nxt[a][j]&&dep[nxt[a][j]]>dep[lca])
ans[i]+=<<j,a=nxt[a][j];
}
for(j=;j>=;j--){
if(nxt[b][j]&&dep[nxt[b][j]]>dep[lca])
ans[i]+=<<j,b=nxt[b][j];
}
if((a==lca&&nxt[b][]&&dep[nxt[b][]]<=dep[lca])||(b==lca&&nxt[a][]&&dep[nxt[a][]]<=dep[lca]))
ans[i]+=;
else if(nxt[a][]&&nxt[b][]&&dep[nxt[a][]]<=dep[lca]&&dep[nxt[b][]]<=dep[lca]){
ans[i]+=;
lq[i<<][]=a,lq[i<<][]=lqh[b],lqh[b]=i<<;
lq[i<<|][]=b,lq[i<<|][]=lqh[a],lqh[a]=i<<|;
}else{
ans[i]=-;
}
}
dfs3();
for(i=;i<=Q;i++)
printf("%d\n",ans[i]-covered[i]);
return ;
}