#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=4e5+10;
int n,a[N],fa[N];
vector<int> g[N];
ll f[N][2][2];
void dfs(int x) {
	f[x][1][0]=-a[x], f[x][1][1]=-a[x]+a[fa[x]];
	ll s=0;
	for(auto y:g[x]) {
		dfs(y);
		f[x][0][0]+=f[y][1][0];
		f[x][1][0]+=max(f[y][0][0]+a[y],f[y][1][0]);
		f[x][1][1]+=max(f[y][0][0]+a[y],f[y][1][0]);
		s+=f[y][1][0];
	}
	for(auto y:g[x]) f[x][1][0]=max(f[x][1][0],s-f[y][1][0]+f[y][1][1]);
	f[x][0][0]=abs(f[x][0][0]), f[x][1][0]=abs(f[x][1][0]), f[x][1][1]=abs(f[x][1][1]);
	return;
}
int main() {
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=2;i<=n;i++) scanf("%d",&fa[i]), g[fa[i]].emplace_back(i);
	dfs(1);
	printf("%lld\n",f[1][1][0]);
	return 0;
}