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;
void setmax(int &x, int y) { x = max(x, y); }
void setmin(int &x, int y) { x = max(x, y); }
void setmax(long long &x, long long y) { x = max(x, y); }
void setmin(long long &x, long long y) { x = max(x, y); }
long long AR = 19, BR = 13, CR = 23, XR = 228, YR = 322, MOD = 1e9 + 993;
const int mod = 1e9 + 7;
long long myrand() {
  long long ZR = (XR * AR + YR * BR + CR) % MOD;
  XR = YR;
  YR = ZR;
  return ZR;
}
string prm(string s) {
  string t;
  for (int i = 1; i < s.length(); i++) t.push_back(s[i]);
  t.push_back(s[0]);
  return t;
}
long long bpow(long long x, long long y) {
  if (y == 0) return 1;
  if (y == 1) return x;
  long long ret = bpow(x, y >> 1);
  ret = (ret * ret) % mod;
  if (y & 1) ret = (ret * x) % mod;
  return ret;
}
long long bdiv(long long x, long long y) {
  return (x * bpow(y, mod - 2)) % mod;
}
int md(int x) {
  if (x >= mod) x -= mod;
  if (x < 0) x += mod;
  return x;
}
long long MD(long long x) {
  x %= MOD;
  if (x < 0) x += MOD;
  return x;
}
long long sqr(long long x) { return x * x; }
const long long llinf = 2e18;
const long double EPS = 1e-6;
const int maxn = 5e3 + 100, inf = 1e9 + 100, sq = 300;
int n;
vector<int> e[maxn];
vector<int> q[maxn];
int a, b, c, d;
int sz[maxn];
int rt[maxn];
bitset<maxn> g[2];
int dfs(int v) {
  if (e[v].empty()) sz[v] = 1;
  if (v == a) {
    rt[v] = 1;
    return 1;
  }
  if (v == b) {
    rt[v] = 2;
    return 2;
  }
  if (v == c) {
    rt[v] = 3;
    return 3;
  }
  if (v == d) {
    rt[v] = 4;
    return 4;
  }
  int ret = 0;
  for (int i = 0; i < e[v].size(); i++)
    q[v][i] = dfs(e[v][i]), sz[v] += sz[e[v][i]], ret = max(ret, q[v][i]);
  rt[v] = ret;
  if (v == 0) return ret;
  if (ret == 1) {
    for (int i = 0; i < q[v].size(); i++)
      if (q[v][i] == 0) g[0] |= g[0] << sz[e[v][i]];
  }
  if (ret == 3) {
    for (int i = 0; i < q[v].size(); i++)
      if (q[v][i] == 0) g[0] |= g[0] << sz[e[v][i]];
  }
  if (ret == 2) {
    for (int i = 0; i < q[v].size(); i++)
      if (q[v][i] == 0) g[1] |= g[1] << sz[e[v][i]];
  }
  if (ret == 4) {
    for (int i = 0; i < q[v].size(); i++)
      if (q[v][i] == 0) g[1] |= g[1] << sz[e[v][i]];
  }
  return ret;
}
int main() {
  ios::sync_with_stdio(0);
  g[0][0] = 1;
  g[1][0] = 1;
  cin >> n >> a >> b >> c >> d;
  a--;
  b--;
  c--;
  d--;
  for (int i = 1; i < n; i++) {
    int x;
    cin >> x;
    e[x - 1].push_back(i);
  }
  for (int i = 0; i < n; i++) q[i].resize(e[i].size());
  dfs(0);
  int szb, szc;
  for (int i = 0; i < e[0].size(); i++)
    if (rt[e[0][i]] == 0)
      g[0] |= g[0] << sz[e[0][i]], g[1] |= g[1] << sz[e[0][i]];
    else if (rt[e[0][i]] == 2)
      szb = sz[e[0][i]];
    else if (rt[e[0][i]] == 3)
      szc = sz[e[0][i]];
  if (sz[0] % 2 == 1) {
    cout << "NO";
    return 0;
  }
  if ((sz[0] - 2) / 2 >= max(szb, szc) && g[0][(sz[0] - 2) / 2 - szb] &&
      g[1][(sz[0] - 2) / 2 - szc])
    cout << "YES";
  else
    cout << "NO";
}