#include<bits/stdc++.h>
#define int long long
#define ffor(i,a,b) for(int i=(a);i<=(b);i++)
#define roff(i,a,b) for(int i=(a);i>=(b);i--)
using namespace std;
const int MAXN=5000+10,MOD=1e9+7;
int n,sze[MAXN],b[MAXN],frac[MAXN],inv[MAXN],tcnt[MAXN];
vector<int> G[MAXN];
int qpow(int base,int p) {
	int ans=1;
	while(p) {
		if(p&1) ans=ans*base%MOD;
		base=base*base%MOD,p>>=1;
	}
	return ans;
}
void dfs1(int u) {
	tcnt[u]=1,sze[u]=1;
	for(auto v:G[u]) dfs1(v),sze[u]+=sze[v],tcnt[u]=tcnt[u]*tcnt[v]%MOD*inv[sze[v]]%MOD;
	tcnt[u]=tcnt[u]*frac[sze[u]-1]%MOD;
	return ;
}
int dp[MAXN][MAXN];
int C(int u,int d) {
	return frac[d]*inv[u]%MOD*inv[d-u]%MOD;	
}
void dfs2(int u) {
	int mul=1;
	for(auto v:G[u]) mul=mul*tcnt[v]%MOD*inv[sze[v]]%MOD;
	for(auto v:G[u]) {
		int insv=qpow(tcnt[v]*inv[sze[v]]%MOD,MOD-2);
		ffor(r,1,n-sze[u]+1) dp[v][r+1]=(dp[v][r+1]+dp[u][r]*mul%MOD*insv%MOD*frac[sze[u]-sze[v]-1]%MOD*C(sze[u]-sze[v]-1,n-sze[u]+1-r+sze[u]-sze[v]-1))%MOD;	
	}
	for(auto v:G[u]) ffor(r,1,n-sze[v]+1) dp[v][r]=(dp[v][r]+dp[v][r-1])%MOD;
	for(auto v:G[u]) dfs2(v);
	return ;
}
signed main() {
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	cin>>n,frac[0]=1;
	ffor(i,1,n) frac[i]=frac[i-1]*i%MOD;
	inv[n]=qpow(frac[n],MOD-2);
	roff(i,n-1,0) inv[i]=inv[i+1]*(i+1)%MOD;
	ffor(i,2,n) {int f;cin>>f,G[f].push_back(i);}
	dfs1(1);
	dp[1][1]=1;
	dfs2(1);
	ffor(i,1,n) {
		int ans=0;
		ffor(j,1,n-sze[i]+1) {
			int cnt=dp[i][j]*tcnt[i]%MOD*C(sze[i]-1,n-j)%MOD;
			if(i==j) ans=(ans+cnt)%MOD;
		}
		cout<<ans<<' ';
	}
	return 0;
}