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;
template <typename T>
inline T qmin(const T a, const T b) {
  return a < b ? a : b;
}
template <typename T>
inline T qmax(const T a, const T b) {
  return a > b ? a : b;
}
template <typename T>
inline void getmin(T &a, const T b) {
  if (a > b) a = b;
}
template <typename T>
inline void getmax(T &a, const T b) {
  if (a < b) a = b;
}
inline void fileio(string s) {
  freopen((s + ".in").c_str(), "r", stdin);
  freopen((s + ".out").c_str(), "w", stdout);
}
const int inf = (int)1e9 + 7;
const long long linf = (long long)1e17 + 7;
const int N = 5005, logN = 21;
int v[N], nxt[N];
int head[N], dep[N], par[N][logN], sz[N];
bool kicked[N];
int n;
inline void addedge(int p, int x, int y) {
  v[p] = y, nxt[p] = head[x], head[x] = p;
}
inline void init(int u, int p) {
  dep[u] = dep[p] + 1;
  par[u][0] = p;
  if (!head[u]) ++sz[u];
  for (int e = head[u]; e; e = nxt[e]) {
    init(v[e], u);
    sz[u] += sz[v[e]];
  }
}
inline int lca(int x, int y) {
  if (dep[x] < dep[y]) return lca(y, x);
  for (int i = logN - 1; i >= 0; --i) {
    if (dep[x] - dep[y] >> i & 1) x = par[x][i];
  }
  if (x == y) return x;
  for (int i = logN - 1; i >= 0; --i) {
    if (par[x][i] != par[y][i]) x = par[x][i], y = par[y][i];
  }
  return par[x][0];
}
inline bool chk(const vector<int> &a, int tot) {
  static bool f[N];
  memset(f, 0, sizeof(bool) * (tot + 1));
  f[0] = 1;
  for (int i = 0; i < a.size(); ++i) {
    for (int j = tot; j >= a[i]; --j) f[j] |= f[j - a[i]];
  }
  return f[tot];
}
int a, b, c, d;
int main() {
  scanf("%d", &n);
  scanf("%d%d%d%d", &a, &b, &c, &d);
  for (register int(i) = (2); (i) <= (n); ++(i)) {
    int p;
    scanf("%d", &p);
    addedge(i - 1, p, i);
  }
  init(1, 0);
  if (sz[1] & 1) return puts("NO"), 0;
  for (register int(i) = (1); (i) <= (logN - 1); ++(i)) {
    for (register int(j) = (1); (j) <= (n); ++(j)) {
      par[j][i] = par[par[j][i - 1]][i - 1];
    }
  }
  {
    int anc = lca(a, b);
    vector<int> all;
    if (lca(anc, c) != anc) return puts("NO"), 0;
    memset(kicked, false, sizeof(kicked));
    kicked[a] = true;
    for (int t = a, p = par[t][0]; p != anc; t = p, p = par[t][0]) {
      kicked[p] = true;
      for (int e = head[p]; e; e = nxt[e]) {
        if (v[e] == t) continue;
        all.push_back(sz[v[e]]);
      }
    }
    kicked[b] = true;
    for (int t = b, p = par[t][0]; p != anc; t = p, p = par[t][0]) {
      kicked[p] = true;
      for (int e = head[p]; e; e = nxt[e]) {
        if (v[e] == t) continue;
        all.push_back(sz[v[e]]);
      }
    }
    for (int t = c; t != anc; t = par[t][0]) kicked[t] = true;
    for (int e = head[anc]; e; e = nxt[e]) {
      if (kicked[v[e]]) continue;
      all.push_back(sz[v[e]]);
    }
    if (!chk(all, sz[1] / 2 - 1)) return puts("NO"), 0;
  }
  {
    int anc = lca(c, d);
    vector<int> all;
    if (lca(anc, b) != anc) return puts("NO"), 0;
    memset(kicked, false, sizeof(kicked));
    kicked[c] = true;
    for (int t = c, p = par[t][0]; p != anc; t = p, p = par[t][0]) {
      kicked[p] = true;
      for (int e = head[p]; e; e = nxt[e]) {
        if (v[e] == t) continue;
        all.push_back(sz[v[e]]);
      }
    }
    kicked[d] = true;
    for (int t = d, p = par[t][0]; p != anc; t = p, p = par[t][0]) {
      kicked[p] = true;
      for (int e = head[p]; e; e = nxt[e]) {
        if (v[e] == t) continue;
        all.push_back(sz[v[e]]);
      }
    }
    for (int t = b; t != anc; t = par[t][0]) kicked[t] = true;
    for (int e = head[anc]; e; e = nxt[e]) {
      if (kicked[v[e]]) continue;
      all.push_back(sz[v[e]]);
    }
    if (!chk(all, sz[1] / 2 - 1)) return puts("NO"), 0;
  }
  return puts("YES"), 0;
}