#include <bits/stdc++.h>
using namespace std;
int n, a[100005], t;
vector<int> g[100005];
int d[100005];
map<int, int> c;
void dfs(int x, int dd) {
  d[x] = dd;
  for (int y : g[x]) dfs(y, dd ^ 1);
}
int main() {
  scanf("%d", &n);
  for (int i = 0; i < (n); i++) scanf("%d", a + i);
  for (int i = 1; i < (n); i++) {
    int p;
    scanf("%d", &p);
    p--;
    g[p].push_back(i);
  }
  dfs(0, 0);
  bool asd = false;
  for (int i = 0; i < (n); i++)
    if (!g[i].size() && !d[i]) {
      asd = true;
      break;
    }
  if (asd)
    for (int i = 0; i < (n); i++) d[i] ^= 1;
  for (int i = 0; i < (n); i++)
    if (d[i]) t ^= a[i];
  long long r = 0;
  if (t == 0) {
    int q[2];
    q[0] = q[1] = 0;
    for (int i = 0; i < (n); i++) q[d[i]]++;
    r += 1LL * q[0] * (q[0] - 1) / 2;
    r += 1LL * q[1] * (q[1] - 1) / 2;
  }
  for (int i = 0; i < (n); i++)
    if (!d[i]) c[a[i]]++;
  for (int i = 0; i < (n); i++)
    if (d[i]) r += c[t ^ a[i]];
  printf("%lld\n", r);
  return 0;
}
