#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5;
struct node{
    int p,v;
    node(int p,int v):p(p),v(v){}
    bool operator < (const node& o)const{
        return v>o.v;
    }
};
int dep[N],fa[N],n,a[N];
int cnt[N],b[N];
bool cmp(int x,int y){
    return a[x]+x<a[y]+y;
}
void solve(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cnt[i]=0;
        b[i]=i;
        dep[i]=0;
    }
    for(int i=2;i<=n;i++){
        cin>>fa[i];
        dep[i]=dep[fa[i]]+1;
        cnt[fa[i]]++;
    }
    for(int i=1;i<=n;i++){
        cin>>a[i];
    }
    sort(b+1,b+1+n,cmp);
    priority_queue<node> pq;
    for(int i=1;i<=n;i++){
        if(cnt[i]==0){
            pq.push(node(i,dep[i]));
        }
    }
    int ans=0;
    for(int i=1;i<=n;i++){
        int x=b[i];
        node t=pq.top();
        pq.pop();
        if(fa[t.p]){
            cnt[fa[t.p]]--;
            if(cnt[fa[t.p]]==0){
                pq.push(node(fa[t.p],dep[fa[t.p]]));
            }
        }
        ans=max(ans,x+a[x]+t.v);
    }
    cout<<ans<<'\n';
}
int main(){
    ios::sync_with_stdio(false);
    int T;
    cin>>T;
    while(T--){
        solve();
    }
    return 0;
}