#include <bits/stdc++.h>
using namespace std;
using i64 = long long;
const i64 inf = 1e18;

void solve(){
    int n;
    cin>>n;
    vector<i64> hp(n+1),edge[n+1];
    for(int i=2;i<=n;i++){
        int u;
        cin>>u;
        edge[u].push_back(i);
    }
    for(int i=1;i<=n;i++){
        cin>>hp[i];
    }

    vector f(n+1,vector(n+1,vector<i64>(2,inf)));

    vector<i64> siz(n+1);
    function<void(int)> dfs=[&](int u){
        f[u][0][1]=0;
        f[u][1][0]=hp[u];
        siz[u]=1;
        for(auto x:edge[u]){
            dfs(x);
            for(int i=siz[u];i>=0;i--){
                for(int j=0;j<=siz[x];j++){
                    f[u][i+j][0]=min(f[u][i+j][0],f[u][i][0]+min(f[x][j][1],f[x][j][0]+hp[x]));
                    f[u][i+j][1]=min(f[u][i+j][1],f[u][i][1]+min(f[x][j][0],f[x][j][1]));
                }
            }
            siz[u]+=siz[x];
        }
    };
    dfs(1);

    for(int i=0;i<=n;i++){
        cout<<min(f[1][n-i][0],f[1][n-i][1])<<" \n"[i==n];
    }
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);

    int t;
    cin>>t;
    while(t--){
        solve();
    }
    return 0;
}