#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const long long N=2e3+10,INF=1e15;
ll n,m,k,a[N],dp[N][N][2],siz[N],sum[N];
bool st[N];
vector<int>p[N];
void dfs(int u)
{
	dp[u][0][0]=0;
	dp[u][1][1]=a[u];
	siz[u]=1;
	sum[u]=a[u];
	for(int x:p[u])dp[u][1][1]+=a[x];
	for(int x:p[u])
	{
		dfs(x);
		//siz[u]+=siz[x];
		sum[u]+=sum[x]+a[x];
		for(int i=siz[u];i>=0;i--)
		{
			for(int j=siz[x];j>=0;j--)
			{
				dp[u][i+j][0]=max(dp[u][i+j][0],dp[u][i][0]+max(dp[x][j][1]+a[x],dp[x][j][0]));
				dp[u][i+j][1]=max(dp[u][i+j][1],dp[u][i][1]+max(dp[x][j][1],dp[x][j][0]));
			}
		}
		siz[u]+=siz[x];
	}
}
void solve()
{
	cin>>n;
	for(int i=1;i<=n;i++)p[i].clear();
	for(int i=2;i<=n;i++)
	{
		int x;cin>>x;
		p[x].push_back(i);
	}
	for(int i=1;i<=n;i++)cin>>a[i];
	for(int i=0;i<=n;i++)
	for(int j=0;j<=n;j++)dp[i][j][0]=dp[i][j][1]=-INF;
	dfs(1);
	for(int i=0;i<=n-1;i++)cout<<sum[1]-max(dp[1][i][0],dp[1][i][1])<<' ';
	cout<<sum[1]-max(dp[1][n][1],dp[1][n][0])<<endl;
}
int main()
{
    cin.tie(0), cout.tie(0), ios::sync_with_stdio(false);
    int T;cin>>T;
    while(T--)solve();
    return 0;
}