#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e6 + 9;
long long gcd(long long a, long long b) { return b ? gcd(b, a % b) : a; }
long long s[maxn];
int f[maxn], F[maxn], p[maxn];
int main() {
  int n;
  scanf("%d", &n);
  for (int i = 1; i <= n; ++i) {
    scanf("%lld", &s[i]);
  }
  for (int i = 2; i <= n; ++i) {
    scanf("%d", &p[i]);
  }
  for (int i = n; i > 1; --i) {
    s[p[i]] += s[i];
  }
  for (int i = n; i; --i) {
    long long t = s[1] / gcd(s[1], s[i]);
    if (t <= n) {
      ++f[t];
    }
  }
  int i, j, ans = 0;
  for (i = n; i; --i)
    for (j = i; (j += i) <= n;) f[j] += f[i];
  for (F[1] = i = 1; i <= n; ++i) {
    if (f[i] == i) {
      ans += F[i];
      j = i * 2;
      for (; j <= n; j += i) {
        F[j] += F[i];
      }
    }
  }
  printf("%d", ans);
  return 0;
}
