#include <bits/stdc++.h>  
using namespace std;  

#define int long long  

void solve() {
	int n;
	cin >> n;

	vector<vector<int>> adj(n + 1);
	for (int i = 2; i <= n; i++) {
		int p;
		cin >> p;
		adj[p].push_back(i);
	}

	vector<int> v(n + 1);
	for (int i = 1; i <= n; i++) {
		cin >> v[i];
	}

	vector<int> sz(n + 1);
	vector<vector<int>> ans;
	vector dp(n + 1, vector(n + 1, vector<int>(2, 1e18)));

	function<void(int)> dfs = [&](int x) {
		sz[x] = 1;
		for (auto y : adj[x]) {
			dfs(y);			
			sz[x] += sz[y];
		}
		dp[x][0][0] = v[x];
		dp[x][1][1] = 0;
		int p = 1;
		for (auto y : adj[x]) {
			vector ndp(sz[x] + 1, vector<int>(2, 1e18));
			// for (int i = 0; i <= sz[x]; i++) {
			// 	for (int j = 0; j <= min(i, sz[y]); j++) {
			// 		ndp[i][0] = min(ndp[i][0], dp[x][i - j][0] + min(dp[y][j][1], dp[y][j][0] + v[y]));
			// 		ndp[i][1] = min(ndp[i][1], dp[x][i - j][1] + min(dp[y][j][1], dp[y][j][0]));
			// 	}
			// }
			for (int i = sz[x]; i >= 0; i--) {
				for (int j = 0; j <= sz[y] && i + j <= sz[x]; j++) {
					ndp[i + j][0] = min(ndp[i + j][0], dp[x][i][0] + min(dp[y][j][1], dp[y][j][0] + v[y]));
					ndp[i + j][1] = min(ndp[i + j][1], dp[x][i][1] + min(dp[y][j][1], dp[y][j][0]));
				}
			}
			dp[x] = ndp;
			p += sz[x];
		}
	};

	dfs(1);

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

signed main() {  
    ios::sync_with_stdio(false);
    cin.tie(nullptr);  

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

    return 0;
}