#include<iostream>
#include<cstdio>
#include<vector>
#include<set>
using namespace std;
const int N=100005;
int n;
int c[N];
int fa[N];
vector<int>G[N];
int dis[N],dep[N],son[N],len[N];
void init(int u,int father)
{
    dis[u]=dep[u]=dep[father]+1;
    son[u]=0;
    for(int v:G[u])
    {
        if(v==father) continue;
        init(v,u);
        if(son[u]==0||dis[v]<dis[son[u]]) son[u]=v;
    }
    if(son[u]) dis[u]=dis[son[u]];
    len[u]=dis[u]-dep[u]+1;
    return;
}
int f[N];
vector<int>pos[N],neg[N],zero[N];
int ans[N];
long long sum[N];
void dfs(int u,int father)
{
    if(son[u]) dfs(son[u],u);
    else
    {
        pos[u].clear(),zero[u].clear(),neg[u].clear();
        if(c[u]==0) f[u]=0,sum[u]=0,pos[u].emplace_back(1);
        else f[u]=1,sum[u]=-1,neg[u].emplace_back(-1);
        ans[u]=0;
        return;
    }
    swap(pos[u],pos[son[u]]);
    swap(zero[u],zero[son[u]]);
    swap(neg[u],neg[son[u]]);
    f[u]=f[son[u]];
    sum[u]=sum[son[u]];
    for(int v:G[u])
    {
        if(v==father||v==son[u]) continue;
        dfs(v,u);
        f[u]+=f[v];
        for(int i=1;i<len[u];i++)
        {
            int d=0;
            if(i<=(int)neg[v].size()) d=neg[v][i-1];
            else if(i<=(int)neg[v].size()+(int)zero[v].size()) d=0;
            else d=pos[v][pos[v].size()-(i-(int)neg[v].size()-(int)zero[v].size())];
            if(d==0) continue;
            if(i<=(int)neg[u].size()) neg[u][i-1]+=d,sum[u]+=d;
            else if(i<=(int)neg[u].size()+(int)zero[u].size()) zero[u][i-(int)neg[u].size()-1]+=d;
            else pos[u][pos[u].size()-(i-(int)neg[u].size()-(int)zero[u].size())]+=d;
        }
    }
    int tot=0;
    for(int i=0;i<(int)zero[u].size()&&zero[u][i]<0;i++)
        sum[u]+=zero[u][i],neg[u].emplace_back(zero[u][i]),tot++;
    for(int i=(int)zero[u].size()-1;i>=0&&zero[u][i]>0;i--)
        pos[u].emplace_back(zero[u][i]),tot++;
    zero[u].resize(zero[u].size()-tot);
    while(!neg[u].empty()&&neg[u].back()>=0)
    {
        if(neg[u].back()==0) zero[u].emplace_back(0);
        else pos[u].emplace_back(neg[u].back());
        sum[u]-=neg[u].back();
        neg[u].pop_back();
    }
    while(!pos[u].empty()&&pos[u].back()<=0)
    {
        if(pos[u].back()==0) zero[u].emplace_back(0);
        else sum[u]+=pos[u].back(),neg[u].emplace_back(pos[u].back());
        pos[u].pop_back();
    }
    if(c[u]==0) pos[u].emplace_back(1);
    else f[u]++,neg[u].emplace_back(-1),sum[u]+=-1;
    ans[u]=f[u]+sum[u];
    return;
}
void solve()
{
    cin>>n;
    string str;
    cin>>str;
    for(int i=1;i<=n;i++)
        c[i]=str[i-1]-'0';
    for(int i=1;i<=n;i++)
        G[i].clear();
    for(int i=2;i<=n;i++)
    {
        cin>>fa[i];
        G[fa[i]].emplace_back(i);
        G[i].emplace_back(fa[i]);
    }
    init(1,0);
    dfs(1,0);
    for(int i=1;i<=n;i++)
        cout<<ans[i]<<" ";
    cout<<"\n";
    return;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr),cout.tie(nullptr);
    int T;
    cin>>T;
    while(T--)
        solve();
    return 0;
}

/*
1
18
111001001101000000
1 2 3 4 5 1 7 1 9 2 10 8 12 6 13 11 16
*/