#include <iostream>
#include <vector>
#include <algorithm>
#include <atcoder/modint>
using namespace std;
using i64 = long long;
using u64 = unsigned long long;
using m32 = atcoder::static_modint<998244353>;
#define rep(i,n) for(int i=0; i<(int)(n); i++)


int main(){
    int N; cin >> N;
    vector<int> P(N, -1);
    vector<vector<int>> E(N);
    for(int i=1; i<N; i++){ int p; cin >> p; P[i] = p-1; }
    for(int i=1; i<N; i++) E[P[i]].push_back(i);
    vector<m32> dp(N, 1);
    for(int i=N-1; i>=1; i--) if(!E[i].empty()){
        for(int e : E[i]) dp[i] *= (dp[e] + 1);
    }
    m32 ans = 1;
    for(auto a : E[0]) ans *= dp[a];
    cout << ans.val() << '\n';
    return 0;
}


struct ios_do_not_sync{
    ios_do_not_sync(){
        std::ios::sync_with_stdio(false);
        std::cin.tie(nullptr);
    }
} ios_do_not_sync_instance;