#include <bits/stdc++.h>
using namespace std;
void addmod(int &a, long long b) {
  a = (a + b);
  if (a >= 1000000007) a -= 1000000007;
}
void mulmod(int &a, long long b) { a = (a * b) % 1000000007; }
template <class T>
bool domin(T &a, const T &b) {
  return a > b ? a = b, 1 : 0;
}
template <class T>
bool domax(T &a, const T &b) {
  return a < b ? a = b, 1 : 0;
}
int gi() {
  int a;
  scanf("%d", &a);
  return a;
}
long long gll() {
  long long a;
  scanf("%lld", &a);
  return a;
}
int main() {
  int n = gi(), bc = 0;
  vector<int> a(n);
  for (int i = 0; i < n; i++) {
    a[i] = gi() - 1;
  }
  for (int i = 0; i < n; i++) {
    bc += gi();
  }
  int cnt = 0;
  vector<bool> vis(n, false);
  for (int i = 0; i < n; i++)
    if (!vis[i]) {
      int xi = i;
      cnt++;
      while (!vis[xi]) {
        vis[xi] = true;
        xi = a[xi];
      }
    }
  if (cnt == 1) cnt = 0;
  cout << (cnt) + ((bc & 1) ? 0 : 1) << endl;
  return 0;
}
