#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=200013;
const int P=998244353;

vector<int>mp[N];
int dep[N];
void dfs(int x,int fa) {
    dep[x]=dep[fa]+1;
    for(auto i:mp[x])if(i^fa)dfs(i,x);
}
int a[N];
void solve(){
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)mp[i].clear();
    for(int i=1,x;i<n;i++)cin>>x,mp[x].push_back(i+1),mp[i+1].push_back(x);
    for(int i=1;i<=n;i++)cin>>a[i],a[i]+=i;
    sort(a+1,a+n+1);
    dfs(1,0);
    int ans=0;
    sort(dep+1,dep+n+1,greater<>());
    for(int i=1;i<=n;i++)ans=max(ans,dep[i]+a[i]-1);
    cout<<ans<<'\n';
}
int main(){
    ios::sync_with_stdio(false);
    int T=1;
    cin>>T;
    while(T--)solve();
}