#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
using namespace std;
long long max(long long a, long long b) {
  if (a > b) {
    return a;
  } else {
    return b;
  }
}
long long min(long long a, long long b) {
  if (a < b) {
    return a;
  } else {
    return b;
  }
}
long long power(long long b, long long e) {
  if (e == 0) return 1;
  if (e % 2)
    return ((b * power((b) * (b), (e - 1) / 2)));
  else
    return power((b) * (b), e / 2);
}
long long modpower(long long b, long long e, long long q) {
  long long MOD = q;
  if (e == 0) return 1;
  if (e % 2)
    return ((b % MOD) * modpower((b % MOD) * (b % MOD), (e - 1) / 2, q)) % MOD;
  else
    return modpower((b % MOD) * (b % MOD), e / 2, q) % MOD;
}
void dpv(vector<long long> v) {
  for (long long i = 0; i < v.size(); i++) {
    cout << v[i] << " ";
  }
  cout << '\n';
}
void dpv(vector<pair<long long, long long> > v) {
  for (long long i = 0; i < v.size(); i++) {
    cout << v[i].first << " " << v[i].second << '\n';
  }
}
void dpv(set<long long> v) {
  for (auto i : v) {
    cout << i << " ";
  }
  cout << '\n';
}
long long md = 1e9 + 7;
long long dp[100005][2];
vector<vector<long long> > adj(100005);
long long color[100005];
void dfs(long long x, long long p) {
  if (color[x] == 1) {
    dp[x][1] = 1;
    dp[x][0] = 0;
  } else {
    dp[x][0] = 1;
    dp[x][1] = 0;
  }
  for (auto u : adj[x]) {
    if (u != p) {
      dfs(u, x);
      long long one = dp[x][1];
      long long zero = dp[x][0];
      dp[x][1] = 0;
      dp[x][0] = 0;
      dp[x][0] += dp[u][0] * zero;
      dp[x][0] %= md;
      dp[x][1] += dp[u][1] * zero + one * dp[u][0];
      dp[x][1] %= md;
      dp[x][1] += one * dp[u][1];
      dp[x][1] %= md;
      dp[x][0] += zero * dp[u][1];
      dp[x][0] %= md;
    }
  }
}
void oblivious() {
  memset(dp, 0, sizeof(dp));
  long long n;
  cin >> n;
  for (long long i = 0; i < n - 1; i++) {
    long long x;
    cin >> x;
    adj[i + 1].push_back(x);
    adj[x].push_back(i + 1);
  }
  for (long long i = 0; i < n; i++) {
    cin >> color[i];
  }
  dfs(0, -1);
  cout << dp[0][1] << '\n';
}
signed main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  long long t = 1;
  while (t--) {
    oblivious();
  }
  return 0;
}
