#include <bits/stdc++.h>

using namespace std;

using int64 = long long;
const int mod = 1e9 + 7;
// const int mod = 998244353;

const int64 infll = (1LL << 62) - 1;
const int inf = (1 << 30) - 1;

struct IoSetup {
  IoSetup() {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(10);
    cerr << fixed << setprecision(10);
  }
} iosetup;


template< typename T1, typename T2 >
ostream &operator<<(ostream &os, const pair< T1, T2 > &p) {
  os << p.first << " " << p.second;
  return os;
}

template< typename T1, typename T2 >
istream &operator>>(istream &is, pair< T1, T2 > &p) {
  is >> p.first >> p.second;
  return is;
}

template< typename T >
ostream &operator<<(ostream &os, const vector< T > &v) {
  for(int i = 0; i < (int) v.size(); i++) {
    os << v[i] << (i + 1 != v.size() ? " " : "");
  }
  return os;
}

template< typename T >
istream &operator>>(istream &is, vector< T > &v) {
  for(T &in : v) is >> in;
  return is;
}

template< typename T1, typename T2 >
inline bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); }

template< typename T1, typename T2 >
inline bool chmin(T1 &a, T2 b) { return a > b && (a = b, true); }

template< typename T = int64 >
vector< T > make_v(size_t a) {
  return vector< T >(a);
}

template< typename T, typename... Ts >
auto make_v(size_t a, Ts... ts) {
  return vector< decltype(make_v< T >(ts...)) >(a, make_v< T >(ts...));
}

template< typename T, typename V >
typename enable_if< is_class< T >::value == 0 >::type fill_v(T &t, const V &v) {
  t = v;
}

template< typename T, typename V >
typename enable_if< is_class< T >::value != 0 >::type fill_v(T &t, const V &v) {
  for(auto &e : t) fill_v(e, v);
}

template< typename F >
struct FixPoint : F {
  FixPoint(F &&f) : F(forward< F >(f)) {}

  template< typename... Args >
  decltype(auto) operator()(Args &&... args) const {
    return F::operator()(*this, forward< Args >(args)...);
  }
};

template< typename F >
inline decltype(auto) MFP(F &&f) {
  return FixPoint< F >{forward< F >(f)};
}

int main() {
  // なもりグラフみたいになっているんですが
  // 基本的には普通にとってもよい
  // 残しちゃダメなやつ: 長さ2のループ1つ

  // 最後の数頂点はうくなので全探索をします

  // これは嘘で, ある1頂点に向かうやつだけのこってもダメです
  // そのような形になったやつは除去していく必要がある
  // ある1頂点に向かうやつ+長さ2のループがめんどい

  // 入次数が残り頂点数のやつがある か？　だと信じています

  int N;
  cin >> N;
  vector< int > A(N);

  vector< int > in(N);

  for(int i = 0; i < N; i++) {
    cin >> A[i];
    --A[i];
    in[A[i]]++;
  }
  multiset< int > XX;
  for(int i = 0; i < N; i++) {
    XX.emplace(in[i]);
  }


  auto solve_uku = [&](vector< int > rest, int last) {

    sort(begin(rest), end(rest));
    do {
      int las = last;
      bool ng = false;
      for(int i = 0; i < rest.size(); i++) {
        if(rest[i] == las) ng = true;
        las = A[rest[i]];
      }
      if(!ng) return rest;
    } while(next_permutation(begin(rest), end(rest)));

    if(N == 2) {
      cout << -1 << endl;
      exit(0);
    }

    throw 0;
    exit(0);
  };

  set< int > st;
  for(int i = 0; i < N; i++) st.emplace(i);
  vector< int > ans;
  int ng = -1;
  while(st.size() >= 10 && *XX.rbegin() + 1 < st.size()) {
    auto it = st.begin();
    if(*it == ng) ++it;
    ans.emplace_back(*it);
    ng = A[*it];
    XX.erase(XX.find(in[A[*it]]));
    in[A[*it]]--;
    XX.emplace(in[A[*it]]);
    st.erase(it);
  }

  if(st.size() > 10) {
    int p = -1;
    for(int i : st) {
      if(in[i] + 1 == st.size()) p = i;
    }
    assert(A[p] != ng);
    bool first = true;
    while(st.size() >= 10) {
      auto it = first ? st.find(p) : st.begin();
      if(*it == ng) ++it;
      ans.emplace_back(*it);
      ng = A[*it];
      in[A[*it]]--;
      st.erase(it);
      first = false;
    }
  }


  vector< int > B;
  for(auto &rest : st) {
    B.emplace_back(rest);
  }
  auto ans2 = solve_uku(B, ng);

  for(auto &p : ans) cout << p + 1 << " ";
  for(auto &p : ans2) cout << p + 1 << " ";

  cout << endl;
}



Measured memory_bytes: 11008