#include <bits/stdc++.h>
using namespace std;
template <class T>
void _R(T &x) {
  cin >> x;
}
void _R(int &x) { scanf("%d", &x); }
void _R(long long &x) { scanf("%lld", &x); }
void _R(double &x) { scanf("%lf", &x); }
void _R(char &x) { scanf(" %c", &x); }
void _R(char *x) { scanf("%s", x); }
void R() {}
template <class T, class... U>
void R(T &head, U &...tail) {
  _R(head);
  R(tail...);
}
template <class T>
void _W(const T &x) {
  cout << x;
}
void _W(const int &x) { printf("%d", x); }
void _W(const long long &x) { printf("%lld", x); }
void _W(const double &x) { printf("%.16f", x); }
void _W(const char &x) { putchar(x); }
void _W(const char *x) { printf("%s", x); }
template <class T, class U>
void _W(const pair<T, U> &x) {
  _W(x.first);
  putchar(' ');
  _W(x.second);
}
template <class T>
void _W(const vector<T> &x) {
  for (auto i = x.begin(); i != x.end(); _W(*i++))
    if (i != x.cbegin()) putchar(' ');
}
void W() {}
template <class T, class... U>
void W(const T &head, const U &...tail) {
  _W(head);
  putchar(sizeof...(tail) ? ' ' : '\n');
  W(tail...);
}
int MOD = 1e9 + 7;
void ADD(long long &x, long long v) {
  x = (x + v) % MOD;
  if (x < 0) x += MOD;
}
const int SIZE = 1e6 + 10;
set<int> e[SIZE];
int p[SIZE];
int h[SIZE];
vector<int> AA;
vector<int> an;
void f(int x) {
  AA.push_back(x);
  if (!x) return;
  while (((int)(e[p[x]]).size()) > 1) {
    an.push_back(x);
    int new_p = -1;
    for (auto y : e[p[x]]) {
      if (y != x) {
        new_p = y;
        break;
      }
    }
    e[p[x]].erase(x);
    p[x] = new_p;
    e[new_p].insert(x);
  }
  f(p[x]);
}
int main() {
  int n;
  R(n);
  int ma = 0, ma_id = 0;
  for (int i = (1); i < (n); ++i) {
    R(p[i]);
    e[p[i]].insert(i);
    h[i] = h[p[i]] + 1;
    if (h[i] > ma) {
      ma = h[i];
      ma_id = i;
    }
  }
  f(ma_id);
  reverse((AA).begin(), (AA).end());
  W(AA);
  W(((int)(an).size()));
  reverse((an).begin(), (an).end());
  W(an);
  return 0;
}
