#include <bits/stdc++.h>
using namespace std;
#define For(Ti,Ta,Tb) for(auto Ti=(Ta);Ti<=(Tb);++Ti)
#define Dec(Ti,Ta,Tb) for(auto Ti=(Ta);Ti>=(Tb);--Ti)
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define range(Tx) begin(Tx),end(Tx)
using ll=long long;
const int N=1e6+5;
int n,q,a[N],fa[N],dep[N],pre[N],desc[N][20],g[N][20][2];
ll f[N][20];
vector<int> e[N],deps[N];
ll query(int u,int k){
	ll ans=0;
	Dec(j,19,0){
		if(k>>j&1){
			int v=desc[u][j];
			ans+=f[u][j]-(ll(g[v][j][1])<<j)+(ll(g[v][j][0])<<j);
			u=v;
		}
	}
	For(j,0,19){
		if(k>>j&1){
			ans+=(ll(g[u][j][1])<<j)-(ll(g[u][j][0])<<j);
		}
	}
	return ans;
}
int main(){
	cin.tie(nullptr)->sync_with_stdio(false);
	cin>>n;
	For(i,1,n){
		cin>>a[i];
	}
	For(i,2,n){
		cin>>fa[i];
		e[fa[i]].push_back(i);
		dep[i]=dep[fa[i]]+1;
	}
	{
		queue<int> q;
		q.push(1);
		while(q.size()){
			int u=q.front();
			q.pop();
			pre[u]=(deps[dep[u]].size()?deps[dep[u]].back():0);
			deps[dep[u]].push_back(u);
			for(int v:e[u]){
				q.push(v);
			}
		}
	}
	For(i,1,n){
		for(int u:deps[i]){
			desc[fa[u]][0]=u;
		}
	}
	Dec(i,n,0){
		for(int u:deps[i]){
			if(!desc[u][0]){
				desc[u][0]=desc[pre[u]][0];
			}
			For(j,1,19){
				desc[u][j]=desc[desc[u][j-1]][j-1];
			}
		}
	}
	Dec(i,n,0){
		for(int u:deps[i]){
			For(j,0,19){
				For(x,0,1){
					g[u][j][x]=((a[u]>>j&1)==x);
					for(int v:e[u]){
						g[u][j][x]+=g[v][j][x];
					}
				}
			}
		}
	}
	Dec(i,n,0){
		for(int u:deps[i]){
			For(j,0,19){
				For(x,0,1){
					g[u][j][x]+=g[pre[u]][j][x];
				}
			}
		}
	}
	Dec(i,n,0){
		for(int u:deps[i]){
			f[u][0]=f[pre[u]][0]+a[u];
			For(j,1,19){
				int v=desc[u][j-1],v_=desc[v][j-1];
				f[u][j]=f[u][j-1]+f[v][j-1];
				f[u][j]-=(1LL*g[v][j-1][1]-g[v_][j-1][1])<<(j-1);
				f[u][j]+=(1LL*g[v][j-1][0]-g[v_][j-1][0])<<(j-1);
			}
		}
	}
	cin>>q;
	For($,1,q){
		int u,k;
		cin>>u>>k;
		cout<<query(u,k+1)-query(pre[u],k+1)<<'\n';
	}
	return 0;
}