#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define db double
#define ldb long double
#define pb push_back
#define pii pair<int, int>
#define FR first
#define SE second
using namespace std;
inline int read() {
  int x = 0; bool op = 0;
  char c = getchar();
  while(!isdigit(c))op |= (c == '-'), c = getchar();
  while(isdigit(c))x = (x << 1) + (x << 3) + (c ^ 48), c = getchar();
  return op ? -x : x;
}
const int N = 2e5 + 10;
int n, S, tim;
ll ans;
ull val[N];
mt19937_64 rnd(114514);
vector<int> G1[N], G2[N];
int nxt[N], vis[N], id[N];
unordered_map<ull, int> cnt[3];
struct Node {int x, y;};
int st[N][25], dep[N], dfn[N];
void dfs(int u, int fa) {
  st[dfn[u] = ++tim][0] = fa; dep[u] = dep[fa] + 1;
  for(int v : G1[u])if(v != fa)dfs(v, u);
  return ;
}
int cmin(int x, int y) {return dep[x] < dep[y] ? x : y;}
void build() {
  for(int j = 1; (1 << j) <= n; j++) {
    for(int i = 1; i + (1 << j) - 1 <= n; i++) {
      st[i][j] = cmin(st[i][j - 1], st[i + (1 << (j - 1))][j - 1]);
    }
  }
  return ;
}
int lca(int u, int v) {
  if(u == v)return u;
  int l = dfn[u], r = dfn[v];
  if(l > r)swap(l, r);
  int t = __lg(r - (l++));
  return cmin(st[l][t], st[r - (1 << t) + 1][t]);
}
int dist(int u, int v) {return dep[u] + dep[v] - 2 * dep[lca(u, v)];}
Node merge(Node x, int u) {
  if(!~x.x && !~x.y)return (Node){-1, -1};
  Node now = x;
  if(dist(x.x, u) > dist(now.x, now.y))now = (Node){x.x, u};
  if(dist(x.y, u) > dist(now.x, now.y))now = (Node){x.y, u};
  vector<int> p = {x.x, x.y, u};
  auto chk = [&](Node x, int u) {return dist(x.x, u) + dist(x.y, u) == dist(x.x, x.y);};
  for(int x : p)if(!chk(now, x))return (Node){-1, -1};
  return now;
}
void dfs2(int u, int fa, ull now, Node p, int op) {
  now ^= val[u]; now ^= val[nxt[u]];
  p = merge(p, u);
  if(~p.x && ~p.y) {
    if(!op) {
      auto add = [&](ull now, int c) {
        for(int i = 0; i + c <= 2; i++)ans += cnt[i][now];
        return ;
      };
      add(now, 0);
      add(now ^ val[p.x], 1);
      // ans += cnt[2][now] + cnt[1][now];
      // ans += cnt[1][now ^ val[p.x]] + cnt[0][now ^ val[p.x]];
      if(p.x != p.y) {
        add(now ^ val[p.y], 1);
        add(now ^ val[p.x] ^ val[p.y], 2);
        // ans += cnt[1][now ^ val[p.y]];
        // ans += cnt[now ^ val[p.x] ^ val[p.y]];
      }
      // printf("qwq:%d %d %d %d\n", u, p.x, p.y, ans);
    }
    else {
      cnt[0][now]++;
      cnt[1][now ^ val[p.x]]++;
      if(p.x != p.y) {
        cnt[1][now ^ val[p.y]]++;
        cnt[2][now ^ val[p.x] ^ val[p.y]]++;
      }
    }
  }
  for(int v : G2[u])if(v != fa && id[v] == tim)dfs2(v, u, now, p, op);
  return ;
}
void calc(int u) {
  for(int i = 0; i <= 2; i++) {
    unordered_map<ull, int>().swap(cnt[i]);
  }
  cnt[0][0] = 1;
  for(int v : G2[u]) {
    // cout << "calc:" << u << ' ' << v << endl;
    if(id[v] != tim)continue;
    dfs2(v, u, 0, (Node){v, v}, 0); 
    dfs2(v, u, 0, (Node){v, v}, 1);
  }
  return ;
}
int root;
int sz[N];
void findroot(int u, int fa) {
  int mx = 0; sz[u] = 1;
  for(int v : G1[u]) {
    if(v == fa || vis[v])continue;
    findroot(v, u); sz[u] += sz[v];
    mx = max(mx, sz[v]);
  }
  if(max(S - sz[u], mx) * 2 <= S)root = u;
  return ;
}
void dfs1(int u, int fa) {
  id[u] = tim; nxt[u] = fa;
  for(int v : G1[u]) {
    if(v == fa || vis[v])continue;
    dfs1(v, u);
  }
  return ;
}
void solve(int u) {
  // if(S > 0)cerr << u << ' ' << S << endl;
  // cerr << "solve:" << u << ' ' << root << ' ' << vis[root] << ' ' << S << endl;
  findroot(u, 0); u = root;

  vis[u] = true;
  tim++; 
  for(int v : G1[u])if(!vis[v])dfs1(v, 0); 
  // printf("solve:%d\n", u);
  // for(int i = 1; i <= n; i++)cout << id[i] << ' '; cout << endl;
  calc(u);
  // printf("%d %d\n", u, ans);
  int nS = S;
  for(int v : G1[u]) {
    if(vis[v])continue;
    S = sz[v] > sz[u] ? nS - sz[u] : sz[v];
    // cerr << u << ' ' << v << endl;
    solve(v);
  }
  return ;
}
int main() { 
  n = read();
  for(int i = 1; i <= n; i++)val[i] = rnd();
  for(int i = 2; i <= n; i++) {
    int u = read();
    G1[u].pb(i); G1[i].pb(u);
  }
  for(int i = 2; i <= n; i++) {
    int u = read();
    // printf("%d %d\n", u, i);
    G2[u].pb(i); G2[i].pb(u);
  }
  // for(int x : G2[2])cout << x << endl;
  dfs(1, 0); build();
  S = ans = n; solve(1);
  printf("%lld\n", ans);
  return 0;
}