// 
#include<bits/stdc++.h>
#define ll long long
#define fi first
#define se second
#define pb push_back
#define ull unsigned long long
#define inf 0x3f3f3f3f3f3f3f3f
using namespace std;
const int N=50005;
int n,fa[N];
vector<int>G[N];
struct node{
    ull x;
    int d;
};
struct node0{
    ull g[64];
    node f[64];
    void ins(ull x){
        for(int i=63;i>=0;i--){
            if(x>>i&1){
                if(!g[i]){
                    g[i]=x;
                    return;
                }
                x^=g[i];
            }
        }
    }
    void ins0(node t){
        for(int i=0;i<=63;i++){
            if(t.x>>i&1){
                if(!f[i].x){
                    f[i]=t;
                    return;
                }
                if(f[i].d>t.d)swap(f[i],t);
                t.x^=f[i].x;
            }
        }
    }
    void build(){
        for(int i=63;i>=0;i--){
            for(int j=i-1;j>=0;j--){
                if(g[i]>>j&1){
                    g[i]^=g[j];
                }
            }
        }
        for(int i=0;i<=63;i++){
            if(!g[i]){
                f[i].x|=1ull<<i;
                for(int j=i+1;j<=63;j++){
                    if(g[j]>>i&1){
                        f[i].x|=1ull<<j;
                    }
                }
            }
        }
    }
    ull ask(int d){
        ull x=0;
        for(int i=63;i>=0;i--){
            if(f[i].x&&d>=f[i].d){
                if(__builtin_parityll(x&f[i].x)){
                    x|=1ull<<i;
                }
            }else x|=1ull<<i;
        }
        return x;
    }
}a[N];
void dfs(int u){
    for(auto v:G[u]){
        dfs(v);
        for(int i=63;i>=0;i--){
            if(a[v].f[i].x)a[u].ins0({a[v].f[i].x,a[v].f[i].d+1});
        }
    }
}
vector<int>vec;
ull res[N];
void dfs0(int u){
    vec.clear(),vec.pb(n),vec.pb(0);
    for(int i=0;i<=63;i++){
        if(a[u].f[i].x){
            vec.pb(a[u].f[i].d);
        }
    }
    sort(vec.begin(),vec.end());
    for(int i=0;i<vec.size()-1;i++){
        if(vec[i]==vec[i+1])continue;
        res[u]+=a[u].ask(vec[i])*(vec[i+1]-vec[i]);
    }
    for(auto v:G[u]){
        for(int i=0;i<=63;i++){
            if(a[u].f[i].x)a[v].ins0({a[u].f[i].x,a[u].f[i].d+1});
        }
        dfs0(v);
    }
}
int main(){
    //freopen("data.in","r",stdin);
    // freopen("own.out","w",stdout);
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    cin>>n;
    for(int i=2;i<=n;i++){
        cin>>fa[i],G[fa[i]].pb(i);
    }
    for(int i=1;i<=n;i++){
        int K;cin>>K;
        for(int j=1;j<=K;j++){
            ull x;cin>>x;
            a[i].ins(x);
        }
        a[i].build();
    }
    dfs(1);
    dfs0(1);
    // for(int i=0;i<=63;i++){
    //     if(a[2].f[i].x){
    //         cout<<"find:"<<i<<" "<<a[2].f[i].x<<" "<<a[2].f[i].d<<"\n";
    //     }
    // }
    // cout<<a[1].ask(0);
    for(int i=1;i<=n;i++){
        cout<<res[i]<<"\n";
    }
}