#include<bits/stdc++.h>
//#pragma GCC optimize(3)
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
using namespace std;
using namespace __gnu_pbds;
#define int long long
#define ll long long
#define ull unsigned ll
#define db double
#define ldb long db
#define mp make_pair
#define fi first
#define se second
#define pii pair<int,int>
#define pbI push_back
#define inf 1e9
#define mdI int mid=(l+r)>>1
#define lowbit(x) (x&-x)
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define MAXN 400010
#define gint __int128
#define MOD 998244353
int qread(){
	char o=getchar();int f=1,x=0;
	for(;!isdigit(o);o=getchar())if(o=='-')f*=-1;
	for(;isdigit(o);o=getchar())x=x*10+o-'0';
	return f*x;
}
void chmx(int &x,int y){if(y>x)x=y;}
void chmi(int &x,int y){if(y<x)x=y;}
bool FLA;
namespace MathPart{
	#define rep(a,b,c) for(int a=b;a<=c;a++)
	#define vep(a,b,c) for(int a=b;a>=c;a--)
	void exgcd(ll a,ll b,ll &x,ll &y){if(b==0){x=1,y=0;return;}exgcd(b,a%b,y,x);y-=(a/b)*x;}
	ll inv(ll a,ll p=MOD){ll x,y;exgcd(a,p,x,y);return (x+p)%p;}
	gint qp(ll a,ll b,ll p=MOD){gint rs=1,bs=a;while(b){if(b&1){rs=rs*bs%p;}bs=bs*bs%p;b>>=1;}return rs;}
	ll fac[MAXN],ifac[MAXN];
	void initF(ll n,ll p=MOD){
		fac[0]=1;rep(i,1,n){fac[i]=fac[i-1]*i%p;}
		ifac[n]=inv(fac[n],p);vep(i,n-1,0){ifac[i]=ifac[i+1]*(i+1)%p;}
	}
	ll C(ll n,ll m){if(n<m){return 0;}return fac[n]*ifac[m]%MOD*ifac[n-m]%MOD;}
	ll CC(ll n,ll m,ll p=MOD){ll rs=1;rep(i,0,m-1){rs=rs*(n-i)%p*inv(i+1,MOD)%p;}return rs;}
	ll S2OL(ll n,ll m){
		ll rs=0,d=1,f=1;
		for(ll k=0;k<=m;k++){
			rs=(rs+f*C(m,k)*qp(m-k,n,MOD)%MOD+MOD)%MOD;
			d=d*max(k,1ll)%MOD;f*=-1;
		}return rs*inv(d,MOD)%MOD;
	}
	ll phi(ll a){ll rs=a,m=sqrt(a);
		rep(i,2,m){if(a%i==0){rs=rs/i*(i-1);while(a%i==0){a/=i;}}}
		if(a>1){rs=rs/a*(a-1);}return rs;
	}
	bool MR(ll pr){if(pr<2)return 0;if(pr==3)return 1;if(pr==2)return 1;
		ll a=pr-1,b=0;while(!(a&1))a>>=1,b++;
		for(int i=1;i<=8;++i){
		    ll x=rand()%(pr-2)+2,v=qp(x,a,pr),j;if(v==1||v==pr-1)continue;
		    for(j=0;j<b-1;j++){v=(gint)v*v%pr;if(v==pr-1)break;}
			if(v!=pr-1)return 0;
		}return 1;
	}
}
using namespace MathPart;
int n,fa[MAXN],sz[MAXN],dep[MAXN];
vector<int>vec[MAXN];
int dp[MAXN][110],te[110];
void dfs(int u){
    rep(i,0,dep[u])dp[u][i]=1;
    for(auto v:vec[u]){
        dep[v]=dep[u]+1;
        dfs(v);
        memset(te,0,sizeof(te));
        rep(j,0,dep[u]){
            rep(k,0,dep[u]-j){
                te[j+k]=(te[j+k]+dp[u][j]*dp[v][k]%MOD*C(sz[u]+sz[v]+j+k,sz[v]+k)%MOD)%MOD;
            }
        }
        rep(j,0,dep[u])dp[u][j]=te[j];
        sz[u]+=sz[v];
    }
    sz[u]++;
    rep(j,0,dep[u])dp[u][j]=dp[u][j+1];
}
bool FLB;
signed main(){
	cerr<<((&FLB)-(&FLA))/1024.0/1024<<endl;
	initF(4e5);
	dep[1]=1;
	cin>>n;rep(i,2,n){fa[i]=qread();vec[fa[i]].pbI(i);dep[i]=dep[fa[i]]+1;}
    dfs(1);
    if(dep[n]>=100){cout<<214465651;return 0;}
    cout<<dp[1][0];
}
/*
15
1 2 2 2 2 3 3 2 7 7 3 10 3 13
各子树不互相影响
*/