#define DEBUG 0
#include <bits/stdc++.h>
#define loop(n) for (lint ngtkana_is_a_genius = 0; ngtkana_is_a_genius < lint(n); ngtkana_is_a_genius++)
#define rep(i, begin, end) for (lint i = lint(begin); (i) < lint(end); i++)
#define all(v) v.begin(), v.end()
#define rand(l, r) std::uniform_int_distribution<>(l, r)(mt)
using lint = long long;
auto mt = std::mt19937_64(std::random_device{}());
auto cmn = [](auto&& a, auto b){ if (a > b) {a = b; return true;} return false; };
auto cmx = [](auto&& a, auto b){ if (a < b) {a = b; return true;} return false; };
void debug_impl() { std::cerr << std::endl; }
template <typename Head, typename... Tail>
void debug_impl(Head head, Tail... tail) { std::cerr << " " << head; debug_impl(tail...); }
#if DEBUG
#define debug(...)\
  do {\
    std::cerr << std::boolalpha << "[" << #__VA_ARGS__ << "]:";\
    debug_impl(__VA_ARGS__);\
    std::cerr << std::noboolalpha;\
  } while (false)
#else
#define debug(...) {}
#endif

template < typename Container, typename Value = typename Container::value_type, std::enable_if_t<!std::is_same< Container, std::string >::value, std::nullptr_t> = nullptr>
std::istream& operator>> (std::istream& is, Container& v)
  { for (auto & x : v) { is >> x; } return is; }

template < typename Container, typename Value = typename Container::value_type, std::enable_if_t<!std::is_same< Container, std::string >::value, std::nullptr_t> = nullptr >
std::ostream& operator<< (std::ostream& os, Container const& v) {
 os << "{";
  for (auto it = v.begin(); it != v.end(); it++)
    {os << (it != v.begin() ? "," : "") << *it;}
  return os << "}";
}

template < template < typename ... > class Tuple,  typename... Args, std::size_t ... Inds, std::size_t = std::tuple_size< Tuple < Args ... > >::value >
std::istream& tuple_input_impl(std::istream& os, Tuple<Args...>& tuple, std::integer_sequence<std::size_t, Inds...>)
  { (void)std::initializer_list<int>{((void)(os >> std::get< Inds >(tuple)), 0)...}; return os; }
template < template < typename ... > class Tuple,  typename... Args, std::size_t = std::tuple_size< Tuple < Args ... > >::value >
std::istream& operator>> (std::istream& os, Tuple<Args...>& tuple)
  { return tuple_input_impl(os, tuple, std::index_sequence_for<Args...>()); }

template < template < typename ... > class Tuple,  typename... Args, std::size_t ... Inds, std::size_t = std::tuple_size< Tuple < Args ... > >::value >
std::ostream& tuple_output_impl(std::ostream& os, const Tuple<Args...>& tuple, std::integer_sequence<std::size_t, Inds...>)
  { os << "("; (void)std::initializer_list<int>{((void)(os << (Inds > 0 ? "," : "") << std::get< Inds >(tuple)), 0)...}; return os << ")"; }
template < template < typename ... > class Tuple,  typename... Args, std::size_t = std::tuple_size< Tuple < Args ... > >::value >
std::ostream& operator<< (std::ostream& os, const Tuple<Args...>& tuple)
 { return tuple_output_impl(os, tuple, std::index_sequence_for<Args...>()); }

void err() { std::cout << -1 << std::endl; exit(0); }
#define require(exp) do { if (!(exp)) { err(); } } while(false)

int main() {
  std::cin.tie(0); std::cin.sync_with_stdio(false);
  int n; std::cin >> n;
  std::vector<int> a(n); std::cin>>a;
  for (int & x : a) x--;
  debug(a);

  if (n==2 && a == std::vector<int>{1,0}) {
    err();
  }
  std::vector<int> ref(n,0);
  for (int x : a) {
    ref.at(x)++;
  }
  debug(ref);

  int prv = -1;
  std::set<int> set;
  rep(i,0,n) set.insert(i);
  std::set<std::pair<int,int>, std::greater<>> refset;
  rep(i,0,n) {
    refset.emplace(ref.at(i), i);
  }
  std::vector<int> ans;

  while (!set.empty()) {
    debug(prv, set, ref, refset);
    int sz = set.size();
    int max, i; std::tie(max, i) = *refset.begin();
    assert(max < sz);

    int nxt = *set.begin();

    if (sz == 3) {
      for (int x : set) {
        int y = a.at(x);
        if (x < y && set.find(y) != set.end() && a.at(y) == x) {
          nxt = (prv == x ? y : x);
        }
      }
    }

    if (max + 1 == sz) {
      nxt = i;
    } else if (prv == nxt) {
      nxt = *std::next(set.begin());
    }

    prv = a.at(nxt);
    ans.emplace_back(nxt);
    set.erase(nxt);

    refset.erase(std::pair<int,int>{ref.at(nxt), nxt});
    ref.at(nxt) = 0;

    auto it = refset.find(std::pair<int,int>{ref.at(a.at(nxt)), a.at(nxt)});
    ref.at(a.at(nxt))--;
    if (it!=refset.end()) {
      refset.erase(it);
      refset.emplace(ref.at(a.at(nxt)), a.at(nxt));
    }
  }

  for (int & x : ans) x++;
  rep(i,0,n) {
    std::cout << (0 < i ? " " : "") << ans.at(i);
  }
  std::cout << std::endl;
  return 0;
}

/*
最後の 2 つが、循環参照をしていたらアウトです。
また、全員に参照されているものがあったら、それはすぐに使わないと最後まで使えなくて詰みます。

各人の参照カウントを管理して、頂点数 - 1 のものがあったら使うをします。
これは、初期状態が 2 点循環参照だったとき以外は絶対に選べます。
また、そうでないときには最小のものを、それが無理なら次を選びます。

参照カウントは、添字をキーにした vector と、値だけの multiset で管理します。
頂点を消したら、それが参照していたものの参照カウントを 1 減らし、自分を 0 にします。
*/


Measured memory_bytes: 11008