#include <cstdio>
using namespace std;
typedef long long ll;
int read(){
	char c=getchar();int x=0;
	while(c<48||c>57) c=getchar();
	do x=(x<<1)+(x<<3)+(c^48),c=getchar();
	while(c>=48&&c<=57);
	return x;
}
const int N=5003,P=998244353;
int qp(int a,int b=P-2){
	int res=1;
	while(b){
		if(b&1) res=(ll)res*a%P;
		a=(ll)a*a%P;b>>=1;
	}
	return res;
}
int n,f[N],s[N],fac[N];
int dp[N],res[N];
int main(){
	n=read();
	fac[0]=1;
	for(int i=1;i<=n;++i) s[i]=1,fac[i]=(ll)fac[i-1]*i%P;
	for(int i=2;i<=n;++i) f[i]=read();
	for(int i=n;i>1;--i) s[f[i]]+=s[i];
	dp[1]=1;
	for(int x=0;x<n;++x){
		int p=read();
		for(int i=n;i;--i){
			res[i]=(res[i]+(ll)dp[i]*fac[n-x-1]%P*p)%P;
			dp[i]=((ll)dp[f[i]]*s[f[i]]+(ll)dp[i]*(n-x-s[i]))%P;
		}
	}
	int prod=1;
	for(int i=1;i<=n;++i) prod=(ll)prod*s[i]%P;
	prod=qp(prod);
	for(int i=1;i<=n;++i){
		res[i]=(ll)res[i]*prod%P*s[i]%P;
		printf("%d ",res[i]);
	}
	putchar('\n');
	return 0;
}