#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define int long long
#define MAXN 2010
#define N 1000000000000
ll hd[MAXN],ed[MAXN*2],nxt[MAXN*2];
ll dp[MAXN][2][MAXN][2],val[MAXN],siz[MAXN],cnt[MAXN];
ll n,idx,T;
inline ll read(){
    ll x=0,f=1;
    char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-')f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
} 
void build(ll x,ll y){
    nxt[++idx]=hd[x];
    hd[x]=idx;
    ed[idx]=y;
}
void dfs(ll now){
//    cout<<now<<endl;
    for(int i=0;i<=n;i++)dp[now][1][i][0]=N,dp[now][0][i][0]=N;
    dp[now][1][1][0]=0,dp[now][0][0][0]=val[now];
    siz[now]=1;
    cnt[now]=0;
    for(int i=hd[now];i;i=nxt[i]){
        ll y=ed[i];
        dfs(y);
        cnt[now]++;
        siz[now]+=siz[y];
        for(int j=0;j<=siz[now];j++){
            dp[now][1][j][cnt[now]%2]=N,dp[now][0][j][cnt[now]%2]=N;
            for(int u=0;u<=min(siz[y],j);u++){
                dp[now][0][j][cnt[now]%2]=min(dp[now][0][j][cnt[now]%2],dp[now][0][j-u][(cnt[now]-1)%2]+min(dp[y][0][u][cnt[y]%2]+val[y],dp[y][1][u][cnt[y]%2]));
                if(j)dp[now][1][j][cnt[now]%2]=min(dp[now][1][j][cnt[now]%2],dp[now][1][j-u][(cnt[now]-1)%2]+min(dp[y][0][u][cnt[y]%2],dp[y][1][u][cnt[y]%2]));
//                cout<<now<<" "<<j<<" "<<u<<" "<<dp[now][0][j][cnt[now]%2]<<" "<<dp[now][1][j][cnt[now]%2]<<endl;
            }
//            cout<<now<<" "<<j<<" "<<dp[now][0][j][cnt[now]%2]<<" "<<dp[now][1][j][cnt[now]%2]<<endl;
        }
    }
}
void clear(){
//    memset(dp,125,sizeof(dp));
//    memset(siz,0,sizeof(siz));
//    memset(cnt,0,sizeof(cnt));
    memset(hd,0,sizeof(hd));
    memset(ed,0,sizeof(ed));
    memset(nxt,0,sizeof(nxt));
    idx=0;
}
signed main(){
    T=read();
    while(T--){
        n=read();
        clear();
        for(int i=1;i<n;i++){
            ll p=read();
            build(p,i+1);
        } 
        for(int i=1;i<=n;i++)val[i]=read();
        dfs(1);
        for(int i=0;i<=n;i++){
            if(!i)cout<<dp[1][0][0][cnt[1]%2]<<" ";
            else if(i==n)cout<<"0\n";
            else cout<<min(dp[1][1][i][cnt[1]%2],dp[1][0][i][cnt[1]%2])<<" "; 
        }
        cout<<endl;
    }
    return 0;
}