#include<bits/stdc++.h>
using namespace std;
#define pb push_back
#define int long long
const int M=3e5+10;
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-')f=-1;ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		x=(x<<3)+(x<<1)+(ch^48);ch=getchar();
	}
	return x*f;
} 
inline void write(int x){
	if(x<0){
		x=-x;putchar('-');
	}
	if(x>9)write(x/10);
	putchar(x%10+'0');
}
int res;
int n;
vector<int>to[M];
int siz[M],son[M],dfn[M],tot,idfn[M];
inline void dfs(const int x){
	dfn[x]=++tot;
	idfn[tot]=x;siz[x]=1;
	for(int v:to[x]){
//		cout<<v;
		dfs(v);
		siz[x]+=siz[v];
		if(siz[son[x]]<=siz[v])son[x]=v;
	}
}
int fa[M];
inline int get(const int x){
	if(fa[x]==x)return x;
	return fa[x]=get(fa[x]);
}
int fs[M];int cnt,gt;
int tree[M<<2];
inline void build(const int x,const int l,const int r){
	if(l==r){
		if(l==0||l==n+1){
			tree[x]=1;
		}else{
			tree[x]=0;
		}
		return ;
	}
	int mid=(l+r)>>1;
	build(x<<1,l,mid);build(x<<1|1,mid+1,r);
	tree[x]=tree[x<<1]+tree[x<<1|1];
}

inline void update(const int x,const int l,const int r,const int pos,const int val){
	if(l==r){
//		cout<<0;
		tree[x]=val;
		return ;
	}
	int mid=(l+r)>>1;
	if(pos<=mid)update(x<<1,l,mid,pos,val);
	else update(x<<1|1,mid+1,r,pos,val);
	tree[x]=tree[x<<1]+tree[x<<1|1];
	return ;
}
inline int q1(const int x,const int l,const int r,const int pos){
	if(tree[x]==0||l>=pos)return -1;
	if(l==r)return l;
	int mid=(l+r)>>1;
	int ans=q1(x<<1|1,mid+1,r,pos);
	if(ans==-1)ans=q1(x<<1,l,mid,pos);
	return ans;
}
inline int q2(const int x,const int l,const int r,const int pos){
	if(tree[x]==0||r<=pos)return -1;
	if(l==r)return l;
	int mid=(l+r)>>1;
	int ans=q2(x<<1,l,mid,pos);
	if(ans==-1)ans=q2(x<<1|1,mid+1,r,pos);
	return ans;
}
inline void add(const int x){
	fa[x]=x,fs[x]=1;cnt++;
	int f=get(x-1);
	if(fs[f]){
		cnt+=fs[x]*fs[f];
		fa[f]=x;
		fs[x]+=fs[f];
	}
	int ff=get(x+1);
	if(fs[ff]){
		cnt+=fs[x]*fs[ff];
		fa[ff]=x;
		fs[x]+=fs[ff];
	}
}
inline void ds(const int x){
	int L=q1(1,0,n+1,x),R=q2(1,0,n+1,x);
	gt+=(R-L-1)*(R-L)/2;
	gt-=(x-L-1)*(x-L)/2+(R-1-x)*(R-x)/2;
	update(1,0,n+1,x,0);
}
inline void del(const int x){
	int L=q1(1,0,n+1,x),R=q2(1,0,n+1,x);
	gt-=(R-L-1)*(R-L)/2;
	gt+=(x-L-1)*(x-L)/2+(R-1-x)*(R-x)/2;
	update(1,0,n+1,x,1);
}
inline void solve(const int x){
	for(int v:to[x]){
		if(v!=son[x]){
			solve(v);
			cnt=0;
			for(int i=dfn[v];i<dfn[v]+siz[v];i++){
				int j=idfn[i];
				fa[j]=j;
				fs[j]=0;
				ds(j);
			}
		}
	}
	if(son[x])solve(son[x]);
	add(x),del(x);
	for(int v:to[x]){
		if(v!=son[x]){
			for(int i=dfn[v];i<dfn[v]+siz[v];i++){
				int j=idfn[i];
				add(j),del(j);
			}
		}
	}
	if(x!=1){
//		cout<<res;
		res-=gt;res-=cnt;
	}
}
signed main(){
	n=read();
	res=n*(n+1)/2*n;
	for(int i=2;i<=n;i++){
		int u=read();
		to[u].pb(i);
	}dfs(1);
	for(int i=1;i<=n+1;i++)fa[i]=i,fs[i]=0;
	cnt=0,gt=n*(n+1)/2;
//	cout<<0;
	build(1,0,n+1);
	solve(1);
	write(res);
	return 0;
}