Problem: Bytelandian Tree Factory produces trees for all kinds of industrial applications. You have been tasked with optimizing the production of a certain type of tree for an especially large and important order.

The tree in question is a rooted tree with n vertices labelled with distinct integers from 0 to n - 1. The vertex labelled 0 is the root of the tree, and for any non-root vertex v the label of its parent p(v) is less than the label of v.

All trees at the factory are made from bamboo blanks. A bamboo is a rooted tree such that each vertex has exactly one child, except for a single leaf vertex with no children. The vertices of a bamboo blank can be labelled arbitrarily before its processing is started.

To process a bamboo into another tree a single type of operation can be made: choose an arbitrary non-root vertex v such that its parent p(v) is not a root either. The operation consists of changing the parent of v to its parent's parent p(p(v)). Note that parents of all other vertices remain unchanged, in particular, the subtree of v does not change.

Efficiency is crucial, hence you have to minimize the number of operations to make the desired tree from a bamboo blank. Construct any optimal sequence of operations to produce the desired tree.

Note that the labelling of the resulting tree has to coincide with the labelling of the desired tree. Formally, the labels of the roots have to be equal, and for non-root vertices with the same label the labels of their parents should be the same.

It is guaranteed that for any test present in this problem an answer exists, and further, an optimal sequence contains at most 10^6 operations. Note that any hack that does not meet these conditions will be invalid.

Input

The first line contains a single integer n — the number of vertices in the tree (2 ≤ n ≤ 10^5).

The second line contains n - 1 integers p(1), …, p(n - 1) — indices of parent vertices of 1, …, n - 1 respectively (0 ≤ p(i) < i).

Output

In the first line, print n distinct integers id_1, …, id_n — the initial labelling of the bamboo blank starting from the root vertex (0 ≤ id_i < n).

In the second line, print a single integer k — the number of operations in your sequence (0 ≤ k ≤ 10^6).

In the third line print k integers v_1, …, v_k describing operations in order. The i-th operation consists of changing p(v_i) to p(p(v_i)). Each operation should be valid, i.e. neither v_i nor p(v_i) can be the root of the tree at the moment.

Examples

Input


5
0 0 1 1


Output


0 2 1 4 3
2
1 3


Input


4
0 1 2


Output


0 1 2 3
0

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
int read() {
  int x = 0;
  bool flg = false;
  char ch = getchar();
  for (; !isdigit(ch); ch = getchar())
    if (ch == '-') flg = true;
  for (; isdigit(ch); ch = getchar()) x = (x << 3) + (x << 1) + (ch ^ 48);
  return flg ? -x : x;
}
int n;
int p[100010];
vector<int> e[100010];
struct C {
  int o, u, v;
  C(int O, int U, int V) { o = O, u = U, v = V; }
};
vector<C> G;
priority_queue<pair<int, int>, vector<pair<int, int> >,
               greater<pair<int, int> > >
    q[100010];
void merge(priority_queue<pair<int, int>, vector<pair<int, int> >,
                          greater<pair<int, int> > > &x,
           priority_queue<pair<int, int>, vector<pair<int, int> >,
                          greater<pair<int, int> > > &y) {
  if (x.size() < y.size()) swap(x, y);
  while (y.size()) x.push(y.top()), y.pop();
}
int gnd[17][100010], dep[100010];
int lift(int x, int d) {
  d = dep[x] - d;
  for (int i(0), _i(16); i <= _i; i++)
    if (d >> i & 1) x = gnd[i][x];
  return x;
}
void dfs0(int x, int d) {
  gnd[0][x] = p[x];
  dep[x] = d;
  for (int j = 1; j < 17 && gnd[j - 1][x]; ++j)
    gnd[j][x] = gnd[j - 1][gnd[j - 1][x]];
  for (const int &y : e[x]) dfs0(y, d + 1);
}
void dfs(int x, int d) {
  if (e[x].empty()) return q[x].push(pair<int, int>(d, x)), void();
  set<int> S;
  for (const int &y : e[x]) {
    dfs(y, d + 1);
    S.insert(y);
    merge(q[x], q[y]);
  }
  const int lim = e[x].size() - 1;
  if (!lim) return;
  vector<pair<int, int> > T;
  {
    for (int _(1), __(lim); _ <= __; _++) T.push_back(q[x].top()), q[x].pop();
  }
  S.erase(lift(T[0].second, d + 1));
  for (int _(0), __(lim - 1); _ <= __; _++) {
    const auto &w = T[_];
    const int u = lift(w.second, d + 1);
    assert(!S.count(u));
    if (_ < lim - 1) {
      const int v = lift(T[_ + 1].second, d + 1);
      if (S.count(v)) {
        S.erase(v);
        G.emplace_back(v, w.second, x);
        continue;
      }
    }
    const int v = *S.begin();
    S.erase(v);
    G.emplace_back(v, w.second, x);
  }
}
void solve() {
  n = read();
  for (int i(2), _i(n); i <= _i; i++) p[i] = read() + 1, e[p[i]].push_back(i);
  dfs0(1, 1);
  dfs(1, 1);
  reverse((G).begin(), (G).end());
  for (const auto &w : G) p[w.o] = w.u;
  reverse((G).begin(), (G).end());
  const function<void(int)> work = [&](int x) {
    static bool vis[100010];
    if (vis[x]) return;
    if (p[x]) work(p[x]);
    vis[x] = 1;
    printf("%d ", x - 1);
  };
  for (int i(1), _i(n); i <= _i; i++) work(i);
  puts("");
  vector<int> H;
  for (const auto &w : G) {
    assert(p[w.o] == w.u);
    while (p[w.o] != w.v) H.push_back(w.o), p[w.o] = p[p[w.o]];
  }
  printf("%d\n", int(H.size()));
  for (int x : H) printf("%d ", x - 1);
  puts("");
}
int main() {
  for (int T = 1; T--;) solve();
  return 0;
}