#include <bits/stdc++.h>
using namespace std;
const int M = 2e5 + 7;
inline int ri() {
  int x = 0;
  bool f = 1;
  char c = getchar();
  for (; !isdigit(c); c = getchar())
    if (c == '-') f = 0;
  for (; isdigit(c); c = getchar()) x = x * 10 + c - 48;
  return f ? x : -x;
}
int n, m, sn;
int sx[M], f[M], val[M];
long long ans[M];
struct vec {
  int g[M], te;
  struct edge {
    int y, nxt;
  } e[M << 1];
  vec() {
    memset(g, 0, sizeof g);
    te = 1;
  }
  inline void push(int x, int y) {
    e[++te] = (edge){y, g[x]};
    g[x] = te;
  }
  inline void push2(int x, int y) {
    push(x, y);
    push(y, x);
  }
  inline int& operator()(int x) { return g[x]; }
  inline edge& operator[](int x) { return e[x]; }
} e;
struct node {
  int l, r, id, x;
  inline bool operator<(const node& y) const {
    if (l / sn != y.l / sn) return l < y.l;
    return r < y.r;
  }
} qs[M];
int tq = 0;
int top[M], pre[M];
int dep[M], sz[M], son[M];
int st[M], ed[M], pid[M], tdfn;
void dfs1(int x) {
  int p, y;
  sz[x] = 1, son[x] = 0;
  for (p = e(x); p; p = e[p].nxt)
    if ((y = e[p].y) != pre[x]) {
      pre[y] = x;
      dep[y] = dep[x] + 1;
      dfs1(y);
      sz[x] += sz[y];
      if (sz[y] > sz[son[x]]) son[x] = y;
    }
}
void dfs2(int x) {
  int p, y;
  pid[st[x] = ++tdfn] = x;
  if (y = son[x]) top[y] = top[x], dfs2(y);
  for (p = e(x); p; p = e[p].nxt)
    if ((y = e[p].y) != pre[x] && y != son[x]) top[y] = y, dfs2(y);
  pid[ed[x] = ++tdfn] = x;
}
inline int LCA(int x, int y) {
  while (top[x] != top[y]) {
    if (dep[top[x]] < dep[top[y]]) swap(x, y);
    x = pre[top[x]];
  }
  return dep[x] < dep[y] ? x : y;
}
int v[M][2];
long long nw = 0;
bool use[M];
inline void upd(int x) {
  if (!use[x]) {
    nw += v[f[x]][sx[x] ^ 1];
    v[f[x]][sx[x]]++;
  } else {
    nw -= v[f[x]][sx[x] ^ 1];
    v[f[x]][sx[x]]--;
  }
  use[x] ^= 1;
}
int main() {
  n = ri();
  sn = 1 + (int)sqrt(n);
  for (int i = (1); i <= (n); ++i) sx[i] = ri();
  for (int i = (1); i <= (n); ++i) f[i] = val[i] = ri();
  sort(val + 1, val + n + 1);
  int tv = unique(val + 1, val + n + 1) - (val + 1);
  for (int i = (1); i <= (n); ++i)
    f[i] = lower_bound(val + 1, val + tv + 1, f[i]) - val;
  for (int i = (1); i < (n); ++i) e.push2(ri(), ri());
  pre[1] = 0, dep[1] = 1, dfs1(1);
  top[1] = 1, tdfn = 0, dfs2(1);
  m = ri();
  for (int i = (1); i <= (m); ++i) {
    int x = ri(), y = ri(), z = LCA(x, y);
    if (x == y) {
      ans[i] = 0;
      continue;
    }
    if (st[x] > st[y]) swap(x, y);
    if (ed[y] < ed[x])
      qs[++tq] = (node){st[x], st[y], i, 0};
    else
      qs[++tq] = (node){ed[x], st[y], i, z};
  }
  sort(qs + 1, qs + tq + 1);
  int l = 1, r = 0;
  for (int i = (1); i <= (tq); ++i) {
    for (int j = (qs[i].r); j > (r); --j) upd(pid[j]);
    for (int j = (qs[i].l); j < (l); ++j) upd(pid[j]);
    for (int j = (r); j > (qs[i].r); --j) upd(pid[j]);
    for (int j = (l); j < (qs[i].l); ++j) upd(pid[j]);
    if (qs[i].x) upd(qs[i].x);
    ans[qs[i].id] = nw;
    if (qs[i].x) upd(qs[i].x);
    l = qs[i].l, r = qs[i].r;
  }
  for (int i = (1); i <= (m); ++i) printf("%I64d\n", ans[i]);
  return 0;
}
