#include <bits/stdc++.h>
using namespace std;
long long int Pow(long long int a, long long int b, long long int md,
                  long long int ans = 1) {
  for (; b; b >>= 1, a = a * a % md)
    if (b & 1) ans = ans * a % md;
  return ans % md;
}
const long long int MAXN = 1e5 + 10;
const long long int INF = 8e18;
const long long int MOD = 1e9 + 7;
long long int dp[MAXN][2], pd[2][2], n;
vector<long long int> adj[MAXN];
void add(long long int &x, long long int y) {
  y = (y % MOD + MOD) % MOD;
  x = (x + y) % MOD;
}
void DFS(long long int v) {
  if ((long long int)adj[v].size() == 0) {
    dp[v][1] = 1;
    return;
  }
  for (long long int u : adj[v]) {
    DFS(u);
  }
  dp[v][1] = 1;
  for (long long int u : adj[v]) {
    long long int x =
        (dp[v][1] * dp[u][1] % MOD + dp[v][0] * dp[u][0] % MOD) % MOD;
    long long int y =
        (dp[v][0] * dp[u][1] % MOD + dp[v][1] * dp[u][0] % MOD) % MOD;
    add(dp[v][0], x);
    add(dp[v][1], y);
  }
  add(dp[v][0], dp[v][0]), add(dp[v][1], dp[v][1]);
  pd[0][1] = pd[0][0] = 1;
  pd[1][1] = pd[1][0] = 0;
  for (long long int u : adj[v]) {
    long long int x = pd[0][0], y = pd[1][0];
    add(pd[1][0], dp[u][0] * x);
    add(pd[0][0], dp[u][0] * y);
    x = pd[0][1], y = pd[1][1];
    add(pd[0][1], dp[u][1] * y);
    add(pd[1][1], dp[u][1] * x);
  }
  add(dp[v][0], -pd[1][1]);
  add(dp[v][1], -pd[1][0] - pd[0][0]);
}
int main() {
  scanf("%lld", &n);
  for (long long int i = 2; i <= n; i++) {
    long long int p;
    scanf("%lld", &p);
    adj[p].push_back(i);
  }
  DFS(1);
  printf("%lld\n", (dp[1][0] + dp[1][1]) % MOD);
  return 0;
}
