#include <bits/stdc++.h>

template <typename T>
void rd(T& x){
    x = 0; char c = getchar(); int f = 1;
    while(!isdigit(c)) f = c == '-' ? -1 : 1, c = getchar();
    while(isdigit(c)) x = x * 10 + c - '0', c = getchar();
    x *= f;
}
template <typename T, typename... Args>
void rd(T& x, Args&... args){rd(x); rd(args...);}

using ll = long long;
const int N = 4.1e5 + 10;
const ll INF = 1e18;

int n, a[N];
ll f[N][3][2];
/*
0/1/2: 初始的权值被父亲吸走/被儿子吸走/自己去吸（初始权值被吸走后还可以操作这个点，这时不影响答案）
相邻点不可能同时选 2 状态
0/1: 当前点是否吸走父亲的初始权值
*/

std::vector <int> G[N];

void dp(int x, int fath){
    for(auto y : G[x]) dp(y, x);
    ll sum1 = 0, sum2 = 0, dim = INF;
    for(auto y : G[x]){
        ll t = std::max({f[y][0][0] - a[y], f[y][1][0], f[y][2][0]});
        sum1 += t;
        sum2 += std::max({f[y][0][0] - a[y], f[y][1][0]});
        dim = std::min(dim, t - std::max(f[y][1][1], f[y][2][1])); // 去掉一个子树的贡献，因为这个子树要吸走 x
    }
    f[x][0][0] = sum1;
    f[x][1][0] = sum1 - dim, f[x][1][1] = sum1 - dim - a[fath];
    f[x][2][0] = sum2 + a[x], f[x][2][1] = sum2 - a[fath];
    sum1 = 0, sum2 = 0, dim = INF;
    for(auto y : G[x]){
        ll t = std::max({f[y][0][0] + a[y], f[y][1][0], f[y][2][0]});
        sum1 += t;
        sum2 += std::max({f[y][0][0] + a[y], f[y][1][0]});
        dim = std::min(dim, t - std::max(f[y][1][1], f[y][2][1])); // 去掉一个子树的贡献，因为这个子树要吸走 x
    }
    f[x][0][0] = std::max(f[x][0][0], sum1);
    f[x][1][0] = std::max(f[x][1][0], sum1 - dim), f[x][1][1] = std::max(f[x][1][1], sum1 - dim + a[fath]);
    f[x][2][0] = std::max(f[x][2][0], sum2 - a[x]), f[x][2][1] = std::max(f[x][2][1], sum2 + a[fath]);
}

int main(){
    rd(n);
    for(int i = 1; i <= n; ++i) rd(a[i]);
    for(int i = 2, x; i <= n; ++i) rd(x), G[x].push_back(i);
    dp(1, 0);
    printf("%lld\n", std::max(f[1][1][0], f[1][2][0]));
    return 0;
}