Problem: Earlier, when there was no Internet, each bank had a lot of offices all around Bankopolis, and it caused a lot of problems. Namely, each day the bank had to collect cash from all the offices.

Once Oleg the bank client heard a dialogue of two cash collectors. Each day they traveled through all the departments and offices of the bank following the same route every day. The collectors started from the central department and moved between some departments or between some department and some office using special roads. Finally, they returned to the central department. The total number of departments and offices was n, the total number of roads was n - 1. In other words, the special roads system was a rooted tree in which the root was the central department, the leaves were offices, the internal vertices were departments. The collectors always followed the same route in which the number of roads was minimum possible, that is 2n - 2.

One of the collectors said that the number of offices they visited between their visits to offices a and then b (in the given order) is equal to the number of offices they visited between their visits to offices b and then a (in this order). The other collector said that the number of offices they visited between their visits to offices c and then d (in this order) is equal to the number of offices they visited between their visits to offices d and then c (in this order). The interesting part in this talk was that the shortest path (using special roads only) between any pair of offices among a, b, c and d passed through the central department.

Given the special roads map and the indexes of offices a, b, c and d, determine if the situation described by the collectors was possible, or not.

Input

The first line contains single integer n (5 ≤ n ≤ 5000) — the total number of offices and departments. The departments and offices are numbered from 1 to n, the central office has index 1.

The second line contains four integers a, b, c and d (2 ≤ a, b, c, d ≤ n) — the indexes of the departments mentioned in collector's dialogue. It is guaranteed that these indexes are offices (i.e. leaves of the tree), not departments. It is guaranteed that the shortest path between any pair of these offices passes through the central department.

On the third line n - 1 integers follow: p2, p3, ..., pn (1 ≤ pi < i), where pi denotes that there is a special road between the i-th office or department and the pi-th department.

Please note the joint enumeration of departments and offices.

It is guaranteed that the given graph is a tree. The offices are the leaves, the departments are the internal vertices.

Output

If the situation described by the cash collectors was possible, print "YES". Otherwise, print "NO".

Examples

Input

5
2 3 4 5
1 1 1 1


Output

YES

Input

10
3 8 9 10
1 2 2 2 2 2 1 1 1


Output

NO

Input

13
13 12 9 7
1 1 1 1 5 5 2 2 2 3 3 4


Output

YES

Note

In the first example the following collector's route was possible: <image>. We can note that between their visits to offices a and b the collectors visited the same number of offices as between visits to offices b and a; the same holds for c and d (the collectors' route is infinite as they follow it each day).

In the second example there is no route such that between their visits to offices c and d the collectors visited the same number of offices as between visits to offices d and c. Thus, there situation is impossible. 

In the third example one of the following routes is: <image>.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int MOD = (int)1e9 + 7;
const int INF = (int)1e9;
const long long LINF = (long long)1e18;
const long double PI = acos((long double)-1);
const long double EPS = 1e-9;
inline long long gcd(long long a, long long b) {
  long long r;
  while (b) {
    r = a % b;
    a = b;
    b = r;
  }
  return a;
}
inline long long lcm(long long a, long long b) { return a / gcd(a, b) * b; }
inline long long fpow(long long n, long long k, int p = MOD) {
  long long r = 1;
  for (; k; k >>= 1) {
    if (k & 1) r = r * n % p;
    n = n * n % p;
  }
  return r;
}
template <class T>
inline int chkmin(T& a, const T& val) {
  return val < a ? a = val, 1 : 0;
}
template <class T>
inline int chkmax(T& a, const T& val) {
  return a < val ? a = val, 1 : 0;
}
inline long long isqrt(long long k) {
  long long r = sqrt(k) + 1;
  while (r * r > k) r--;
  return r;
}
inline long long icbrt(long long k) {
  long long r = cbrt(k) + 1;
  while (r * r * r > k) r--;
  return r;
}
inline void addmod(int& a, int val, int p = MOD) {
  if ((a = (a + val)) >= p) a -= p;
}
inline void submod(int& a, int val, int p = MOD) {
  if ((a = (a - val)) < 0) a += p;
}
inline int mult(int a, int b, int p = MOD) { return (long long)a * b % p; }
inline int inv(int a, int p = MOD) { return fpow(a, p - 2, p); }
inline int sign(long double x) { return x < -EPS ? -1 : x > +EPS; }
inline int sign(long double x, long double y) { return sign(x - y); }
const int maxn = 5000 + 5;
int n, a, b, c, d;
int f[maxn];
int p[maxn];
vector<int> adj[maxn];
int size[maxn];
int g[maxn];
bitset<maxn> can[maxn];
int dfs(int u) {
  int res = f[u] ? u + 1 : 0;
  size[u] = !int((adj[u]).size());
  for (int i = (0); i < (int((adj[u]).size())); i++) {
    int v = adj[u][i];
    chkmax(res, dfs(v));
    size[u] += size[v];
  }
  return res;
}
void upd(bitset<maxn>& bt, vector<int> vsz) {
  for (int i = (0); i < (int((vsz).size())); i++) {
    bt |= bt << vsz[i];
  }
}
void solve() {
  cin >> n >> a >> b >> c >> d;
  a--, b--, c--, d--;
  f[a] = f[b] = f[c] = f[d] = 1;
  map<int, int> hs;
  hs[a] = 0, hs[b] = 1, hs[c] = 2, hs[d] = 3;
  for (int i = (1); i < (n); i++) {
    cin >> p[i], p[i]--;
    adj[p[i]].push_back(i);
  }
  int leaves = 0;
  for (int i = (0); i < (n); i++) leaves += !int((adj[i]).size());
  if (leaves & 1) {
    cout << "NO\n";
    return;
  }
  vector<int> vsz;
  for (int i = (0); i < (int((adj[0]).size())); i++) {
    int u = adj[0][i];
    int r = dfs(u);
    if (!r) {
      vsz.push_back(size[u]);
    } else {
      g[hs[r - 1]] = size[u];
    }
  }
  int ar[] = {a, b, c, d};
  for (int i = (0); i < (4); i++) {
    int u = ar[i];
    can[i][0] = 1;
    while (p[u]) {
      int v = p[u];
      vector<int> vsz;
      for (int j = (0); j < (int((adj[v]).size())); j++) {
        int w = adj[v][j];
        if (w != u) {
          vsz.push_back(size[w]);
        }
      }
      upd(can[i], vsz);
      u = v;
    }
  }
  bitset<maxn> exbt;
  exbt[0] = 1, upd(exbt, vsz);
  int tot = 0, dif = leaves / 2 - 1;
  for (int i = (0); i < (maxn); i++)
    for (int j = (0); j < (maxn); j++) {
      if (can[0][i] && can[1][j] && i + j + g[2] <= dif &&
          exbt[dif - i - j - g[2]]) {
        tot++;
        i = maxn - 1;
        break;
      }
    }
  for (int i = (0); i < (maxn); i++)
    for (int j = (0); j < (maxn); j++) {
      if (can[2][i] && can[3][j] && i + j + g[0] <= dif &&
          exbt[dif - i - j - g[0]]) {
        tot++;
        i = maxn - 1;
        break;
      }
    }
  if (tot == 2) {
    cout << "YES\n";
  } else {
    cout << "NO\n";
  }
}
int main() {
  int JUDGE_ONLINE = 1;
  if (fopen("in.txt", "r")) {
    JUDGE_ONLINE = 0;
    assert(freopen("in.txt", "r", stdin));
  } else {
    ios_base::sync_with_stdio(0), cin.tie(0);
  }
  solve();
  if (!JUDGE_ONLINE) {
  }
  return 0;
}