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;
struct fast_ios {
  fast_ios() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(20);
  };
} fast_ios_;
template <typename T>
istream &operator>>(istream &is, vector<T> &v) {
  for (auto &x : v) is >> x;
  return is;
}
template <typename T>
ostream &operator<<(ostream &os, vector<T> &v) {
  for (long long i = 0; i < v.size(); i++) {
    cout << v[i];
    if (i != v.size() - 1) cout << endl;
  };
  return os;
}
template <typename T>
void Out(T x) {
  cout << x << endl;
}
template <typename T1, typename T2>
void Ans(bool first, T1 y, T2 n) {
  if (first)
    Out(y);
  else
    Out(n);
}
using vec = vector<long long>;
using mat = vector<vec>;
using Pii = pair<long long, long long>;
using PiP = pair<long long, Pii>;
using PPi = pair<Pii, long long>;
using bools = vector<bool>;
using pairs = vector<Pii>;
const long long mod = 1000000007;
struct edge {
  long long to, cost, id;
};
class Graph {
 public:
  long long N;
  vector<vector<edge>> G;
  vec depth;
  vec par;
  vec ans, vs;
  Graph(long long N) : N(N) {
    G = vector<vector<edge>>(N, vector<edge>(0));
    depth.resize(N);
    par.resize(N);
  }
  void add_Directed_edge(long long from, long long to, long long cost = 1,
                         long long id = 0) {
    G[from].push_back(edge({to, cost, id}));
  }
  void add_Undirected_edge(long long v1, long long v2, long long cost = 1,
                           long long id = 0) {
    add_Directed_edge(v1, v2, cost, id);
    add_Directed_edge(v2, v1, cost, id);
  }
  void dfs_for_depth(long long v = 0, long long p = -1) {
    depth[v] = 0;
    par[v] = p;
    for (auto e : G[v]) {
      dfs_for_depth(e.to, v);
      depth[v] = max(depth[v], depth[e.to] + 1);
    }
  }
  void dfs(long long v = 0, long long prev = -1) {
    pairs p;
    for (auto e : G[v]) p.push_back(Pii(depth[e.to], e.to));
    sort(p.begin(), p.end());
    reverse(p.begin(), p.end());
    for (auto e : p) {
      if (prev != -1) ans.push_back(prev);
      dfs(e.second, prev);
      prev = e.second;
    }
    vs.push_back(v);
    if (v == 0) {
      reverse(vs.begin(), vs.end());
      reverse(ans.begin(), ans.end());
    }
  }
};
signed main() {
  long long n;
  cin >> n;
  Graph G(n);
  for (long long i = (1); i < (n); i++) {
    long long p;
    cin >> p;
    G.add_Directed_edge(p, i);
  }
  G.dfs_for_depth();
  G.dfs();
  for (long long i = (0); i < (n); i++) cout << G.vs[i] << ' ';
  cout << endl;
  Out(((long long)G.ans.size()));
  for (long long i = (0); i < (((long long)G.ans.size())); i++)
    cout << G.ans[i] << ' ';
  cout << endl;
  return 0;
}