#include<bits/stdc++.h>
using namespace std ;
typedef long long ll ;
const int mod=1000000007 ;
ll n,ans ;
vector<int> t[200005] ;
ll f[200005] ;
inline void dfs(int s) {
	cout<<s<<" " ;
	if(!sizeof(t[s]))return ;
	for(int i=0;i<sizeof(t[s]);++i) {
		dfs(t[s][i]) ;
		f[s]=(f[s]+f[t[s][i]])%mod ;
	}
}
int main() {
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0) ;
	cin>>n ;
	f[1]=1 ;
	for(int i=2,x;i<=n;++i) {
		cin>>x ;
		t[x].push_back(i) ;
		f[i]=1 ;
		cout<<x<<" "<<i<<" "<<n<<endl ;
	}
	dfs(1) ;//这里直接re了
	for(int i=1;i<=n;++i) {
		ans=(ans+f[i])%mod ;
	}
	cout<<ans ;
	return 0 ;
}

求说明