#define _CRT_SECURE_NO_WARNINGS 1

#include<stdio.h>
#include<vector>
using namespace std;
const int N=200005;

#define ll long long
const ll P=998244353;

inline ll Q(ll x,ll y=P-2,ll z=1){
	for(;y;y>>=1,x=x*x%P)
		(y&1)&&(z=z*x%P);

	return z;
}

int n,sz[N];
ll fc[N],nf[N];
vector<int> e[N];
vector<ll> f,dp[N];

inline ll C(int x,int y){
	if(x>=y&&y>=0)
		return fc[x]*nf[y]%P*nf[x-y]%P;
	else return 0;
}

void DP(int x,int d){
	d++;

	for(int v:e[x]) DP(v,d);

	int i,j;
	dp[x].resize(d);
	f.resize(d);

	for(i=0;i<d;i++) dp[x][i]=1;

	for(int v:e[x]){
		for(i=0;i<d;i++){
			f[i]=dp[x][i];
			dp[x][i]=0;
		}

		sz[x]+=sz[v];

		for(i=0;i<d;i++) for(j=0;i+j<d;j++)
			dp[x][i+j]=(dp[x][i+j]+f[i]*dp[v][j]%P*C(sz[x]+i+j,sz[v]+j))%P;
	}

	sz[x]++;

	for(i=0;i+1<d;i++) dp[x][i]=dp[x][i+1];
}

signed main(){
	int i,u;
	scanf("%d",&n);

	for(i=2;i<=n;i++){
		scanf("%d",&u);
		e[u].push_back(i);
	}

	fc[0]=1;

	for(i=1;i<N;i++) fc[i]=fc[i-1]*i%P;

	nf[N-1]=Q(fc[N-1]);

	for(i=N-1;i>0;i--) nf[i-1]=nf[i]*i%P;

	DP(1,1);
	printf("%lld\n",dp[1][0]);
	return 0;
}