#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
typedef unsigned long long ull;
mt19937_64 rnd(0x66ccff);
ull Rnd(ull x){
	x^=x<<13;
	x^=x>>7;
	x^=x<<17;
	return x;
}
int gcd(int x,int y){
	return y?gcd(y,x%y):x;
}
int n,m,a[N],sz[N],b[N],ans;
ull f[N],V,va[N];
vector<int>e[N];
void dfs(int x){
	sz[x]=1,f[x]=V;
	for(auto v:e[x]){
		dfs(v);
		sz[x]+=sz[v];
		f[x]+=Rnd(f[v]);
	}
	va[sz[x]]=f[x];
}
bool check(int x,int y){
	if(sz[x]<y)return true;
	if(sz[x]==y)return f[x]==va[y];
	ull o=V;
	for(auto v:e[x]){
		if(!check(v,y))return false;
		if(sz[v]<y)o+=Rnd(f[v]);
	}
	return o==va[y];
}
void solve(){
	V=rnd(),ans=1;
	for(int i=1;i<=n;++i)e[i].clear();
	for(int i=2,x;i<=n;++i){
		scanf("%d",&x);
		e[x].push_back(i);
	}
	dfs(1);
	for(int i=1;i<=n;++i)b[i]=sz[i];
	sort(b+1,b+n+1);
	for(int i=n-1;i;--i)b[i]=gcd(b[i],b[i+1]);
	m=unique(b+1,b+n+1)-b-1;
	for(int i=2;i<m;++i)if(check(1,b[i]))++ans;
	printf("%d\n",ans);
}
int main(){
	while(scanf("%d",&n)!=EOF)solve();
	return 0;
}