#include <bits/stdc++.h>

using namespace std;

const int MAXN = 1e6 + 5;
const long long MOD = 998244353;

vector<int> g[MAXN], sz(MAXN);

long long bp(int a, int b) {
    if (b == 0) {
        return 1;
    }
    if (b % 2) {
        return bp(a, b - 1) * a % MOD;
    }
    long long tmp = bp(a, b / 2);
    return tmp * tmp % MOD;
}

void dfs(int v, int p) {
    sz[v] = 1;
    for (auto u : g[v]) {
        if (u == p) {
            continue;
        }
        dfs(u, v);
        sz[v] += sz[u];
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int n;
    cin >> n;
    int k;
    cin >> k;
    cin >> k;
    vector<int> p(n + 1);
    for (int i = 2; i <= n; i++) {
        cin >> p[i];
        g[p[i]].push_back(i);
    }
    dfs(1, 0);
    long long ans = 1;
    int cnt = sz[k];
    for (auto u : g[k]) {
        cnt -= sz[u];
        ans = (ans + (bp(2, sz[u]) - 1) * (bp(2, cnt) - 1) % MOD) % MOD;
    }
    cout << ans;
}

/*
4 1
1
1 1 1
*/
\
