#include<bits/stdc++.h>
using namespace std;
void Solve(){   
    int n;
    cin>>n;
    vector<vector<int>> e(n+1);
    vector<int> a(n+1),dep(n+1),f(n+1); 
    for(int i=2;i<=n;i++){
        int x;
        cin>>x;
        e[i].push_back(x),e[x].push_back(i);
    }
    for(int i=1;i<=n;i++)cin>>a[i],a[i]+=i;
    sort(a.begin()+1,a.end());
    auto dfs=[&](auto&& self,int x,int fa)->void{
        for(auto y:e[x]){
            if(y==fa)continue;
            dep[y]=dep[x]+1;
            self(self,y,x);
        }
    };
    dfs(dfs,1,0);
    sort(dep.begin()+1,dep.end());
    int ans=0;
    for(int i=1;i<=n;i++)
        ans=max(ans,a[i]+dep[n-i+1]);
    cout<<ans<<"\n";
}
int main(){
    ios::sync_with_stdio(false),cin.tie(0);
    int T;
    cin>>T;
    while(T--)Solve();
    return 0;
}