#include<bits/stdc++.h>
#define L(i, j, k) for(int i = (j); i <= (k); ++i)
#define R(i, j, k) for(int i = (j); i >= (k); --i)
#define ll long long 
#define vi vector < int > 
#define sz(a) ((int) (a).size())
#define ll long long 
#define ull unsigned long long
#define me(a, x) memset(a, x, sizeof(a)) 
using namespace std;
const int N = 5e4 + 7;
struct lb {
	ull f[64];
	lb() {
		me(f, 0);
	}
	ull & operator [] (int x) {
		return f[x];
	}
	void insert (ull x) {
		R(i, 63, 0) if(x >> i & 1) {
			if(!f[i]) f[i] = x;
			x ^= f[i];
		}
	}
	bool check(ull x) {
		R(i, 63, 0) if(x >> i & 1) {
			if(!f[i]) return false;
			x ^= f[i];
		}
		return true;
	}
	void clear () {
		me(f, 0);
	}
	int size () {
		int cnt = 0;
		L(i, 0, 63) if(f[i]) cnt += 1;
		return cnt;
	}
	int mxp () {
		R(i, 63, 0) if(f[i]) return i;
		return -1;
	}
	ull mxv() {
		ull ret = 0;
		R(i, 63, 0) if((ret ^ f[i]) > ret) ret ^= f[i];
		return ret;
	} 
} pre[100], now, F[N];
lb operator | (lb a, lb b) {
	R(i, 63, 0) if(a[i]) b.insert(a[i]);
	return b;
}
ull X[66];
lb R[66];
lb operator & (lb a, lb b) {
	int ma = a.mxp(), mb = b.mxp();
	if(ma > mb) swap(a, b), swap(ma, mb);
	if(ma == -1) return a;
	if(ma != mb) return b[mb] = 0, a & b;
	ull xay = a[ma] ^ b[mb], sw = a[ma];
	a[ma] = b[mb] = 0;
	auto ns = a & b;
	int pn = 0;
	R(i, 63, 0) if(a[i]) X[++pn] = a[i]; 
	R[pn + 1] = b;
	R(i, pn, 1) 
		R[i] = R[i + 1], R[i].insert(X[i]);
	if(!R[1].check(xay)) return ns;
	L(i, 1, pn) if(!R[i + 1].check(xay)) xay ^= X[i], sw ^= X[i];
	ns[ma] = sw;
	return ns;
}
lb ans[N];
int dep[N];
vi e[N];
ull ANS[N];
vi fn[N];
void dfs(int x, int f) {
	fn[dep[x]].emplace_back(x);
	for(auto v : e[x]) if(v != f) {
		dep[v] = dep[x] + 1, ans[v] = ans[x] & F[v], dfs(v, x);
	}
}
struct Depu {
	vector < pair < int, lb > > Fun;
	void reduce() {
		int p = 0;
		L(i, 0, sz(Fun) - 1) {
			if(!i || Fun[i].second.size() != Fun[i - 1].second.size()) {
				swap(Fun[p], Fun[i]);
				++p;
			}
		}
		if(p != sz(Fun)) {
			Fun.resize(p); 
		}
	}
	friend Depu operator & (Depu u, Depu v) {
		int x = 0, y = 0;
		Depu ans;
		while(x < sz(u.Fun) && y < sz(v.Fun)) {
			if(u.Fun[x].first <= v.Fun[y].first) {
				ans.Fun.emplace_back(u.Fun[x].first, 
					y > 0 ? u.Fun[x].second & v.Fun[y - 1].second : 
						u.Fun[x].second);
				++x;
			} else {
				ans.Fun.emplace_back(v.Fun[y].first, 
					x > 0 ? u.Fun[x - 1].second & v.Fun[y].second : 
						v.Fun[y].second);
				++y;
			}
		}
		ans.reduce();
		return ans;
	}
	Depu shift(int k) {
		Depu qwq = (*this);
		for(auto&u : qwq.Fun)
			u.first += k;
		return qwq;
	}
};
Depu dp[N];

int n, fa[N];
void dfs1(int x, int fa) {
	for(auto &v : e[x]) if(v != fa) {
		dfs1(v, x), dp[x] = dp[x] & dp[v].shift(1);
	} 
}
void dfs2(int x, int fa) {
	for(auto &v : e[x]) if(v != fa) {
		dp[v] = dp[v] & dp[x].shift(1), dfs2(v, x);
	} 
}
int main() {
	ios :: sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin >> n;
	L(i, 2, n) {
		cin >> fa[i];
		e[fa[i]].emplace_back(i);
		e[i].emplace_back(fa[i]);
	}
	L(i, 1, n) {
		int cnt;
		cin >> cnt;
		while(cnt--) {
			ull w;
			cin >> w;
			F[i].insert(w);
		} 
		dp[i].Fun.emplace_back(0, F[i]);
		dp[i].Fun.emplace_back(n, lb());
	}
	dfs1(1, 0); 
	dfs2(1, 0);
	L(i, 1, n) {
		ull ans = 0;
		L(j, 0, sz(dp[i].Fun) - 2) 
			ans += dp[i].Fun[j].second.mxv() * (dp[i].Fun[j + 1].first - dp[i].Fun[j].first);
		cout << ans << '\n';
	}
	return 0;
}