#include <bits/stdc++.h>
using namespace std;
const int mod = 998244353;
const int base = 19260817;
const long long INF = 1e18;
const int maxn = 2e5 + 10;
const double eps = 1e-8;
int add(int x, int y) {
  x += y;
  return x >= mod ? x - mod : x;
}
int sub(int x, int y) {
  x -= y;
  return x < 0 ? x + mod : x;
}
int mul(int x, int y) { return 1LL * x * y % mod; }
int dp[maxn][3];
vector<int> g[maxn];
long long quickpow(long long a, long long x) {
  long long res = 1;
  while (x) {
    if (x & 1) res = res * a % mod;
    a = a * a % mod;
    x >>= 1;
  }
  return res;
}
void dfs(int u) {
  if (g[u].empty()) {
    dp[u][1] = dp[u][2] = 1;
    return;
  }
  dp[u][0] = dp[u][1] = dp[u][2] = 1;
  int temp = 0;
  for (auto v : g[u]) {
    dfs(v);
    dp[u][0] = mul(dp[u][0], add(dp[v][0], dp[v][1]));
    dp[u][1] = mul(dp[u][1], add(dp[v][0], add(dp[v][1], dp[v][2])));
    temp = add(temp, mul(dp[v][2], quickpow(add(dp[v][0], dp[v][1]), mod - 2)));
  }
  dp[u][2] = dp[u][1], temp = mul(dp[u][0], temp);
  dp[u][1] = sub(dp[u][1], add(dp[u][0], temp));
  dp[u][2] = sub(dp[u][2], dp[u][0]);
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n;
  cin >> n;
  for (int i = 2; i <= n; i++) {
    int p;
    cin >> p;
    g[p].push_back(i);
  }
  dfs(1);
  cout << add(dp[1][0], dp[1][1]);
  return 0;
}
