#include <bits/stdc++.h>
using namespace std;
const int N=100100;
vector<int>G[N];
vector<long long>c,w;
int v,h[N];
long long f[N][2][2][2],g[2][2][2],o;
void dfs(int x){
    if(G[x].empty()){
        f[x][0][0][0]=f[x][1][1][1]=0;
        return h[x]=++v,void();
    }
    for(int y : G[x]){
        dfs(y);
        if(!h[x])
            for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) for(int k=0;k<=1;k++)
                f[x][i][j][k]=max(f[y][i^1][j][k]+c[y-1],f[y][i][j][k]);
        else{
            for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) for(int k=0;k<=1;k++)
                g[i][j][k]=-1e18;
            for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) for(int k=0;k<=1;k++)
            for(int a=0;a<=1;a++) for(int b=0;b<=1;b++) for(int z=0;z<=1;z++)
                g[i][j][k]=max(g[i][j][k],f[x][i][j][a]+f[y][z][b][k]+c[y-1]*(i^z)+w[h[x]-1]*(a^b));
            for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) for(int k=0;k<=1;k++)
                f[x][i][j][k]=g[i][j][k];
        } h[x]=h[y];
    }
}
long long place_police(vector<int> P, vector<long long> C, vector<long long> W){
    for(int i=0;i<P.size();i++)
        G[P[i]].push_back(i+1);
    memset(f,0xcf,sizeof(f));
    c=C,w=W,dfs(0),o=-1e18;
    for(int i=0;i<=1;i++) for(int j=0;j<=1;j++) for(int k=0;k<=1;k++)
        o=max(o,f[0][i][j][k]+W.back()*(j^k)); o=-o;
    for(int x : C) o+=x; for(int x : W) o+=x;
    return o;
}