#include <bits/stdc++.h>
// #define int long long
#define pii pair<int,int>
#define fi first
#define se second
#define mkp make_pair
#define vi vector<int>
#define pb emplace_back
using namespace std;
constexpr int N=1e6+5;
typedef long long ll;
constexpr ll inf=1e18;
int n,r[N],c[N],ir[N];
struct node{
	int f[N][2],h[N],p[N],val[N];
	node()=default;
	void init(int r[])
	{
		for(int i=1;i<=n;i++)
		{
			f[i][1]=f[i-1][0]+r[i];
			f[i][0]=max(f[i-1][0],f[i-1][1]);
			h[i]=max(f[i][1],f[i][0]);
			val[i]=r[i];
		}
		for(int i=n,fl=0,las=-1;i>=1;i--)
		{
			if(!r[i])
				p[i]=i,fl=0;
			else
			{
				if(!fl)fl=1,las=i;
				p[i]=las;
			}
		}
	}
	ll calc(int l,int r)
	{
		if(r-l+1<=-2)return -inf;
		if(r-l+1<=0)return 0;
		return h[r]-h[min(r,p[l])]+(val[l]?(min(r,p[l])-l)/2+1:0);
	}
}v[2];
int a[N];
int idx;
pii id[N<<1];
vi e[N<<1];
int las[N];
ll f[N<<1][2];
void dfs(int u)
{
	for(auto v:e[u])
	{
		dfs(v);
		ll mx=0;
		mx=::v[c[id[v].se]].calc(id[u].fi+1,id[v].fi-1)+f[v][0];
		mx=max(mx,::v[c[id[v].se]].calc(id[u].fi+1,id[v].fi-2)+f[v][1]);
		f[u][0]+=mx;
		// cout<<"v:"<<v<<' '<<f[v][0]<<' '<<f[v][1]<<'\n';
		// cout<<"mx1:"<<mx<<' ';
		mx=::v[c[id[v].se]].calc(id[u].fi+2,id[v].fi-1)+f[v][0];
		mx=max(mx,::v[c[id[v].se]].calc(id[u].fi+2,id[v].fi-2)+f[v][1]);
		// cout<<"mx2:"<<mx<<'\n';
		f[u][1]+=mx;
	}
	f[u][1]+=r[id[u].fi]^c[id[u].se];
}
signed main()
{
#ifndef ONLINE_JUDGE
	freopen("data.in","r",stdin);
	freopen("wa.out","w",stdout);
#endif
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>r[i],ir[i]=(r[i]^1);
	v[0].init(r);
	v[1].init(ir);
	for(int i=1;i<=n;i++)
		cin>>c[i];
	for(int i=2;i<=n;i++)
		cin>>a[i];
	for(int i=2;i<=n;i++)
	{
		id[++idx]=mkp(i-1,a[i]);
		if(las[a[i]])e[las[a[i]]].pb(idx);
		las[a[i]]=las[i]=idx;
	}
	for(int i=1;i<=n;i++)
		id[++idx]=mkp(n,i),e[las[i]].pb(idx);
	// for(int i=1;i<=idx;i++)
	// {
	// 	for(auto v:e[i])
	// 	{
	// 		cerr<<i<<' '<<v<<'\n';
	// 	}
	// }
	dfs(1);
	cout<<max(f[1][0],f[1][1]);
	return 0;
}