#include<bits/stdc++.h>
using namespace std;
#define int long long
const int MOD=998244353;
const int N=3e3+5; 
vector<int>out[N];
vector<int>outt[N];
void add(int u,int v){
	out[u].push_back(v);
}
int n,k;
int lf[N],lfcnt=0;
void dfs(int start,int father){
	bool flag=0;
	int lenn=out[start].size();
	for(int i=0;i<lenn;i++){
		dfs(out[start][i],start);
		flag=1;
	}
	if(!flag)lf[++lfcnt]=start;
}
int d[(1<<20)+5][22];
bool con[25][25];
void case20(){
	dfs(1,0); 
	for(int i=1;i<=n;i++){
		int lenn=out[i].size();
		for(int j=0;j<lenn;j++)con[i][out[i][j]]=con[out[i][j]][i]=1;
	}
	for(int i=1;i<=lfcnt;i++){
		for(int j=i+1;j<=lfcnt;j++){
			if(j-i<=k||i+lfcnt-j<=k)con[lf[i]][lf[j]]=con[lf[j]][lf[i]]=1;
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)if(con[i][j]&&i!=j)outt[i].push_back(j);
	}
	/*for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++)cout<<con[i][j]<<" ";
		cout<<"\n";
	}*/
	d[1][1]=1;
	for(int i=1;i<(1<<n);i++){
		for(int j=1;j<=n;j++){
			if(i&(1<<(j-1))){
				//cout<<i<<" "<<j<<" "<<d[i][j]<<"\n";
				int lenn=outt[j].size();
				for(int k=0;k<lenn;k++){
					if(i&(1<<(outt[j][k]-1)))continue;
					d[i|(1<<(outt[j][k]-1))][outt[j][k]]+=d[i][j];
					d[i|(1<<(outt[j][k]-1))][outt[j][k]]%=MOD;
				}
			}
		}
	}
	int ans=0;
	for(int i=2;i<=n;i++){
		if(!con[1][i])continue;
		ans+=d[(1<<n)-1][i];
		ans%=MOD;
	}
	cout<<ans/2<<"\n";
}
int f[N][3][3];
int tmp[N];
int tmplen=0;
int pre[N],suf[N];
void dp(int start,int father){
	int son=0;
	int lenn=out[start].size();
	for(int i=0;i<lenn;i++){
		dp(out[start][i],start);
		son++;
	}
	if(!son){
		for(int i=0;i<3;i++){
			for(int j=0;j<3;j++){
				if(i!=j)f[start][i][j]=1;
			}
		}
		lfcnt++;
		return;
	}
	if(son==1){
		int v=out[start][0];
		f[start][0][1]=f[v][0][1];
		f[start][0][2]=f[v][0][2];
		f[start][1][0]=f[v][1][0];
		f[start][2][0]=f[v][2][0];
		return;
	}
	tmplen=0;
	for(int i=0;i<lenn;i++){
		tmp[++tmplen]=out[start][i];
	}
	pre[0]=1;
	for(int i=1;i<=tmplen;i++){
		pre[i]=pre[i-1]*f[tmp[i]][1][2]%MOD;
	}
	suf[tmplen+1]=1;
	for(int i=tmplen;i>=1;i--){
		suf[i]=suf[i+1]*f[tmp[i]][1][2]%MOD;
	}
	int ans=0;
	for(int i=0;i<=tmplen-2;i++){
		ans+=pre[i]*suf[i+3]%MOD*f[tmp[i+1]][1][0]%MOD*f[tmp[i+2]][0][2]%MOD;
		ans%=MOD;
	}
	f[start][1][2]=f[start][2][1]=ans;
	ans=f[tmp[1]][0][2];
	for(int i=2;i<=tmplen;i++){
		ans=ans*f[tmp[i]][1][2]%MOD;
	}
	f[start][0][2]=f[start][2][0]=ans;
	ans=f[tmp[tmplen]][0][1];
	//if(start==2)cout<<ans<<"?\n";
	for(int i=1;i<tmplen;i++){
		ans=ans*f[tmp[i]][2][1]%MOD;
		//if(start==2)cout<<tmp[i]<<"?\n";
	}
	//if(start==2)cout<<ans<<"?\n";
	f[start][0][1]=f[start][1][0]=ans;
}
void casek1(){
	dp(1,0);
	tmplen=0;
	int lenn=out[1].size();
	for(int i=0;i<lenn;i++){
		tmp[++tmplen]=out[1][i];
	}
	if(tmplen==1){
		cout<<0<<"\n";
		return;
	}
	if(lfcnt==2){
		cout<<1<<"\n";
		return;
	}
	if(tmplen==2){
		cout<<(f[tmp[1]][0][2]*f[tmp[2]][1][0]+f[tmp[1]][0][1]*f[tmp[2]][2][0])%MOD<<"\n";
		return;
	}
	/*for(int i=1;i<=n;i++){
		cout<<i<<":";
		for(int j=0;j<=2;j++){
			for(int k=0;k<=2;k++){
				cout<<f[i][j][k]<<" ";
			}
			cout<<"\n";
		}
	}*/
	pre[0]=1;
	for(int i=1;i<=tmplen;i++){
		pre[i]=pre[i-1]*f[tmp[i]][1][2]%MOD;
	}
	suf[tmplen+1]=1;
	for(int i=tmplen;i>=1;i--){
		suf[i]=suf[i+1]*f[tmp[i]][1][2]%MOD;
	}
	int ans=0;
	for(int i=0;i<=tmplen-2;i++){
		ans+=pre[i]*suf[i+3]%MOD*f[tmp[i+1]][1][0]%MOD*f[tmp[i+2]][0][2]%MOD;
		ans%=MOD;
	}
	int sum=f[tmp[1]][0][2]*f[tmp[tmplen]][1][0]%MOD;
	for(int i=2;i<tmplen;i++){
		sum=sum*f[tmp[i]][1][2]%MOD;
	}
	//cout<<ans<<" "<<sum<<"\n";
	cout<<(ans+sum)%MOD<<"\n";
}
signed main(){
	//freopen("6.txt","r",stdin);
	scanf("%lld%lld",&n,&k);
	for(int i=1;i<n;i++){
		int v;
		scanf("%lld",&v);
		add(v,i+1); 
	} 
	for(int i=1;i<=n;i++)sort(out[i].begin(),out[i].end());
	if(n<=20){
		case20();
		return 0;
	}
	if(k==1){
		casek1();
		return 0;
	}
	return 0;
}