#include <bits/stdc++.h>
using namespace std;
const int RLEN = 1 << 18 | 1;
inline char nc() {
  static char ibuf[RLEN], *ib, *ob;
  (ib == ob) && (ob = (ib = ibuf) + fread(ibuf, 1, RLEN, stdin));
  return (ib == ob) ? -1 : *ib++;
}
inline int rd() {
  char ch = nc();
  int i = 0, f = 1;
  while (!isdigit(ch)) {
    if (ch == '-') f = -1;
    ch = nc();
  }
  while (isdigit(ch)) {
    i = (i << 1) + (i << 3) + ch - '0';
    ch = nc();
  }
  return i * f;
}
const int N = 1e6 + 50, L = 21;
int n, a[N], c[N], li[N], ri[N], pre[N], lst[N], mn[N][21];
vector<int> unvis;
inline int ga(int *anc, int x) {
  return (anc[x] == x) ? x : (anc[x] = ga(anc, anc[x]));
}
inline void solve(int l, int r) {
  vector<int> vec;
  int p = -1;
  for (int i = ga(ri, l); i <= r; i = ga(ri, i + 1)) {
    vec.push_back(i);
    ++p;
    while (vec.size() >= 3) {
      if (((!a[vec[p - 2]] && a[i]) || (a[vec[p - 2]] && !a[i]) ||
           (a[vec[p - 2]] && (a[vec[p - 2]] == a[i]))) &&
          a[vec[p - 1]]) {
        if (a[vec[p - 2]] != a[vec[p]])
          a[vec[p - 2]] = a[i] = a[vec[p - 2]] + a[i];
        vec.pop_back();
        vec.pop_back();
        p -= 2;
      } else
        break;
    }
  }
  for (int j = 1; j < vec.size(); ++j)
    if (a[vec[j]] && a[vec[j - 1]]) {
      puts("no");
      exit(0);
    }
  deque<int> q;
  int nowl = vec[0];
  for (int j = 1; j < vec.size() - 1; ++j) q.push_back(vec[j]);
  while (!q.empty()) {
    int u = q.front(), v = q.back();
    if (q.size() == 1) {
      if (!a[u]) {
        if (unvis.empty()) {
          puts("no");
          exit(0);
        }
        a[u] = unvis.back();
        unvis.pop_back();
      }
      break;
    } else {
      if (!a[u] && !a[v]) {
        if (unvis.empty()) {
          puts("no");
          exit(0);
        }
        a[u] = a[v] = unvis.back();
        unvis.pop_back();
        nowl = a[u];
        q.pop_front();
        q.pop_back();
      } else if (a[u]) {
        a[q[1]] = nowl;
        q.pop_front();
        q.pop_front();
      } else {
        a[q[q.size() - 2]] = nowl;
        q.pop_back();
        q.pop_back();
      }
    }
  }
  for (int i = ga(ri, l); i <= r; i = ga(ri, i + 1))
    if (i != l) ri[i] = r + 1;
  li[r] = l;
}
auto no = []() {
  puts("no");
  exit(0);
};
vector<set<int>> g;
vector<bool> used;
void dfs(int u, int prev) {
  used[u] = true;
  for (auto v : g[u])
    if (v != prev) {
      if (used[v]) {
        no();
      } else {
        dfs(v, u);
      }
    }
}
void check(vector<int> &tour) {
  int n = ((int)tour.size() + 1) / 2;
  g.resize(n);
  used.resize(n, false);
  for (int i = 0; i < (int)tour.size() - 1; i++) {
    int a = tour[i];
    int b = tour[i + 1];
    if (g[a].count(b) == 0) {
      g[a].insert(b);
      g[b].insert(a);
    }
  }
  dfs(0, -1);
  for (int i = 0; i < n; i++)
    if (!used[i]) no();
}
int main() {
  cin >> n;
  for (int i = 1; i < 2 * n; i++) {
    cin >> a[i];
    ++c[a[i]];
  }
  for (int i = 1; i <= n; i++)
    if (!c[i]) unvis.push_back(i);
  if (a[1] && a[2 * n - 1] && a[1] != a[2 * n - 1]) {
    puts("no");
    return 0;
  } else if (!a[1] && !a[2 * n - 1]) {
    if (unvis.empty()) {
      puts("no");
      return 0;
    }
    a[1] = a[2 * n - 1] = unvis.back();
    unvis.pop_back();
  } else if (a[1] != a[2 * n - 1]) {
    a[1] = a[2 * n - 1] = a[1] + a[2 * n - 1];
  }
  for (int i = 1; i < 2 * n; i++) {
    if (!a[i]) continue;
    pre[i] = (lst[a[i]] ? lst[a[i]] : i);
    lst[a[i]] = i;
  }
  for (int i = 1; i < 2 * n; i++) {
    mn[i][0] = (pre[i] ? pre[i] : i);
    if (a[i] && ((i - pre[i]) % 2)) {
      puts("no");
      return 0;
    }
  }
  for (int i = 1; i <= 2 * n; i++) ri[i] = i, li[i] = i;
  for (int i = 2; i < 2 * n; i++)
    if (a[i] && ga(li, pre[i]) < i) solve(ga(li, pre[i]), i);
  vector<int> tour;
  for (int i = 1; i < 2 * n; i++) tour.push_back(a[i] - 1);
  check(tour);
  puts("yes");
  for (int i = 1; i < 2 * n; i++) printf("%d ", a[i]);
}

### Prompt

Please provide a CPP coded solution to the problem described below:

Euler is a little, cute squirrel. When the autumn comes, he collects some reserves for winter. The interesting fact is that Euler likes to collect acorns in a specific way. A tree can be described as n acorns connected by n - 1 branches, such that there is exactly one way between each pair of acorns. Let's enumerate the acorns from 1 to n.

The squirrel chooses one acorn (not necessary with number 1) as a start, and visits them in a way called "Euler tour" (see notes), collecting each acorn when he visits it for the last time.

Today morning Kate was observing Euler. She took a sheet of paper and wrote down consecutive indices of acorns on his path. Unfortunately, during her way to home it started raining and some of numbers became illegible. Now the girl is very sad, because she has to present the observations to her teacher.

"Maybe if I guess the lacking numbers, I'll be able to do it!" she thought. Help her and restore any valid Euler tour of some tree or tell that she must have made a mistake.

Input

The first line contains a single integer n (1 ≤ n ≤ 5 ⋅ 10^5), denoting the number of acorns in the tree.

The second line contains 2n - 1 integers a_1, a_2, …, a_{2n-1} (0 ≤ a_i ≤ n) — the Euler tour of the tree that Kate wrote down. 0 means an illegible number, which has to be guessed.

Output

If there is no Euler tour satisfying the given description, output "no" in the first line.

Otherwise, on the first line output "yes", and then in the second line print the Euler tour which satisfy the given description.

Any valid Euler tour will be accepted, since the teacher doesn't know how exactly the initial tree looks.

Examples

Input

2
1 0 0


Output

yes
1 2 1


Input

4
1 0 3 2 0 0 0


Output

yes
1 4 3 2 3 4 1


Input

5
0 1 2 3 4 1 0 0 0


Output

no

Note

An Euler tour of a tree with n acorns is a sequence of 2n - 1 indices of acorns. such that each acorn occurs at least once, the first and the last acorns are same and each two consecutive acorns are directly connected with a branch.

### Response

```cpp
#include <bits/stdc++.h>
using namespace std;
const int RLEN = 1 << 18 | 1;
inline char nc() {
  static char ibuf[RLEN], *ib, *ob;
  (ib == ob) && (ob = (ib = ibuf) + fread(ibuf, 1, RLEN, stdin));
  return (ib == ob) ? -1 : *ib++;
}
inline int rd() {
  char ch = nc();
  int i = 0, f = 1;
  while (!isdigit(ch)) {
    if (ch == '-') f = -1;
    ch = nc();
  }
  while (isdigit(ch)) {
    i = (i << 1) + (i << 3) + ch - '0';
    ch = nc();
  }
  return i * f;
}
const int N = 1e6 + 50, L = 21;
int n, a[N], c[N], li[N], ri[N], pre[N], lst[N], mn[N][21];
vector<int> unvis;
inline int ga(int *anc, int x) {
  return (anc[x] == x) ? x : (anc[x] = ga(anc, anc[x]));
}
inline void solve(int l, int r) {
  vector<int> vec;
  int p = -1;
  for (int i = ga(ri, l); i <= r; i = ga(ri, i + 1)) {
    vec.push_back(i);
    ++p;
    while (vec.size() >= 3) {
      if (((!a[vec[p - 2]] && a[i]) || (a[vec[p - 2]] && !a[i]) ||
           (a[vec[p - 2]] && (a[vec[p - 2]] == a[i]))) &&
          a[vec[p - 1]]) {
        if (a[vec[p - 2]] != a[vec[p]])
          a[vec[p - 2]] = a[i] = a[vec[p - 2]] + a[i];
        vec.pop_back();
        vec.pop_back();
        p -= 2;
      } else
        break;
    }
  }
  for (int j = 1; j < vec.size(); ++j)
    if (a[vec[j]] && a[vec[j - 1]]) {
      puts("no");
      exit(0);
    }
  deque<int> q;
  int nowl = vec[0];
  for (int j = 1; j < vec.size() - 1; ++j) q.push_back(vec[j]);
  while (!q.empty()) {
    int u = q.front(), v = q.back();
    if (q.size() == 1) {
      if (!a[u]) {
        if (unvis.empty()) {
          puts("no");
          exit(0);
        }
        a[u] = unvis.back();
        unvis.pop_back();
      }
      break;
    } else {
      if (!a[u] && !a[v]) {
        if (unvis.empty()) {
          puts("no");
          exit(0);
        }
        a[u] = a[v] = unvis.back();
        unvis.pop_back();
        nowl = a[u];
        q.pop_front();
        q.pop_back();
      } else if (a[u]) {
        a[q[1]] = nowl;
        q.pop_front();
        q.pop_front();
      } else {
        a[q[q.size() - 2]] = nowl;
        q.pop_back();
        q.pop_back();
      }
    }
  }
  for (int i = ga(ri, l); i <= r; i = ga(ri, i + 1))
    if (i != l) ri[i] = r + 1;
  li[r] = l;
}
auto no = []() {
  puts("no");
  exit(0);
};
vector<set<int>> g;
vector<bool> used;
void dfs(int u, int prev) {
  used[u] = true;
  for (auto v : g[u])
    if (v != prev) {
      if (used[v]) {
        no();
      } else {
        dfs(v, u);
      }
    }
}
void check(vector<int> &tour) {
  int n = ((int)tour.size() + 1) / 2;
  g.resize(n);
  used.resize(n, false);
  for (int i = 0; i < (int)tour.size() - 1; i++) {
    int a = tour[i];
    int b = tour[i + 1];
    if (g[a].count(b) == 0) {
      g[a].insert(b);
      g[b].insert(a);
    }
  }
  dfs(0, -1);
  for (int i = 0; i < n; i++)
    if (!used[i]) no();
}
int main() {
  cin >> n;
  for (int i = 1; i < 2 * n; i++) {
    cin >> a[i];
    ++c[a[i]];
  }
  for (int i = 1; i <= n; i++)
    if (!c[i]) unvis.push_back(i);
  if (a[1] && a[2 * n - 1] && a[1] != a[2 * n - 1]) {
    puts("no");
    return 0;
  } else if (!a[1] && !a[2 * n - 1]) {
    if (unvis.empty()) {
      puts("no");
      return 0;
    }
    a[1] = a[2 * n - 1] = unvis.back();
    unvis.pop_back();
  } else if (a[1] != a[2 * n - 1]) {
    a[1] = a[2 * n - 1] = a[1] + a[2 * n - 1];
  }
  for (int i = 1; i < 2 * n; i++) {
    if (!a[i]) continue;
    pre[i] = (lst[a[i]] ? lst[a[i]] : i);
    lst[a[i]] = i;
  }
  for (int i = 1; i < 2 * n; i++) {
    mn[i][0] = (pre[i] ? pre[i] : i);
    if (a[i] && ((i - pre[i]) % 2)) {
      puts("no");
      return 0;
    }
  }
  for (int i = 1; i <= 2 * n; i++) ri[i] = i, li[i] = i;
  for (int i = 2; i < 2 * n; i++)
    if (a[i] && ga(li, pre[i]) < i) solve(ga(li, pre[i]), i);
  vector<int> tour;
  for (int i = 1; i < 2 * n; i++) tour.push_back(a[i] - 1);
  check(tour);
  puts("yes");
  for (int i = 1; i < 2 * n; i++) printf("%d ", a[i]);
}

```