// #pragma GCC target("avx")  // CPU 処理並列化
// #pragma GCC optimize("O3")  // CPU 処理並列化
// #pragma GCC optimize("unroll-loops")  // 条件処理の呼び出しを減らす
#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<deque>
#include<stack>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<bitset>
#include<stdlib.h>
#include<cassert>
#include<time.h>
#include<bitset>
#include<numeric>
#include<unordered_set>
#include<complex>
using namespace std;
const long long mod=1000000007;
const long long inf=mod*mod;
const long long d2=(mod+1)/2;
const double EPS=1e-10;
const double INF=1e+10;
const double PI=acos(-1.0);
const int C_SIZE = 3100000;
const int UF_SIZE = 210000;
namespace{
	long long fact[C_SIZE];
	long long finv[C_SIZE];
	long long inv[C_SIZE];
	long long Comb(int a,int b){
	 	if(a<b||b<0)return 0;
	 	return fact[a]*finv[b]%mod*finv[a-b]%mod;
	}
	void init_C(int n){
		fact[0]=finv[0]=inv[1]=1;
		for(int i=2;i<n;i++){
			inv[i]=(mod-(mod/i)*inv[mod%i]%mod)%mod;
		}
		for(int i=1;i<n;i++){
			fact[i]=fact[i-1]*i%mod;
			finv[i]=finv[i-1]*inv[i]%mod;
		}
	}
	long long pw(long long a,long long b){
		if(a<0LL)return 0;
		if(b<0LL)return 0;
		long long ret=1;
		while(b){
			if(b%2)ret=ret*a%mod;
			a=a*a%mod;
			b/=2;
		}
		return ret;
	}
	long long pw_mod(long long a,long long b,long long M){
		if(a<0LL)return 0;
		if(b<0LL)return 0;
		long long ret=1;
		while(b){
			if(b%2)ret=ret*a%M;
			a=a*a%M;
			b/=2;
		}
		return ret;
	}
	int pw_mod_int(int a,int b,int M){
		if(a<0)return 0;
		if(b<0)return 0;
		int ret=1;
		while(b){
			if(b%2)ret=(long long)ret*a%M;
			a=(long long)a*a%M;
			b/=2;
		}
		return ret;
	}
	int ABS(int a){return max(a,-a);}
	long long ABS(long long a){return max(a,-a);}
	double ABS(double a){return max(a,-a);}
	int sig(double r) { return (r < -EPS) ? -1 : (r > +EPS) ? +1 : 0; }
	int UF[UF_SIZE];
	void init_UF(int n){
		for(int i=0;i<n;i++)UF[i]=-1;
	}
	int FIND(int a){
		if(UF[a]<0)return a;
		return UF[a]=FIND(UF[a]);
	}
	void UNION(int a,int b){
		a=FIND(a);b=FIND(b);if(a==b)return;
		if(UF[a]>UF[b])swap(a,b);
		UF[a]+=UF[b];UF[b]=a;
	}
}
// ここから編集しろ
int p[110000];
int ret[110000];
typedef unsigned long long wolf;
wolf v[2000];
int N;
int sk[110000];
int LV;
int LA;
int rem;
int dfs(int a,int b){
	if(a==N){
		return 1;
	}

	sk[a]=-1;
	if(rem==0&&(v[LA/64]&(1LL<<(LA%64)))){
		sk[a]=0;
		v[LA/64]-=(1LL<<(LA%64));
		int tmp=0;
		ret[a]=LA;
		rem-=1000000;
		if(a==0||p[b]!=LA)tmp=dfs(a+1,LA);
		rem+=1000000;
		v[LA/64]+=(1LL<<(LA%64));
		return tmp;
	}
	for(int i=(a?sk[a-1]:0);i<(N/64+1);i++){
		if(!v[i])continue;
		if(sk[a]==-1)sk[a]=i;
		for(int j=0;j<64;j++){
			if(v[i]&(1LL<<j)){
				if(b!=-1&&i*64+j==p[b])continue;
				v[i]-=1LL<<j;
				ret[a]=i*64+j;
				if(p[i*64+j]!=LA)rem--;
				int tmp=dfs(a+1,i*64+j);
				if(p[i*64+j]!=LA)rem++;
				if(tmp)return 1;
				v[i]+=1LL<<j;
			}
		}
	}
	return 0;
}
int h[110000];
int main(){
	int a;scanf("%d",&a);
	N=a;
	for(int i=0;i<a;i++){
		scanf("%d",p+i);
		p[i]--;
		h[p[i]]++;
	}

	if(a==2){
		printf("-1\n");return 0;
	}
	for(int i=0;i<a;i++){
		if(h[i]>LV){
			LV=h[i];
			LA=i;
		}
	}
	rem=N-1-LV;
	// if(rem>100)rem=1145141919;
	// for(int i=0;i<a;i++)S.insert(i);
	for(int i=0;i<a;i++){
		v[i/64]|=(1LL<<i);
	}
	// if(tar!=-1){
	// 	for(int i=0;i<a;i++){
	// 		v[i/64]|=(1LL<<i);
	// 	}
	// 	v[tar/64]-=(1LL<<(tar%64));
	// 	ret[0]=tar;
	// 	int res=dfs(1,tar);
	// 	if(res==0){
	// 		printf("-1\n");return 0;
	// 	}
	// 	for(int i=0;i<a;i++){
	// 		if(i)printf(" ");
	// 		printf("%d",ret[i]+1);
	// 	}printf("\n");
	// 	return 0;

	// }
	int res=dfs(0,-1);
	if(res==0){
		printf("-1\n");return 0;
	}
	for(int i=0;i<a;i++){
		if(i)printf(" ");
		printf("%d",ret[i]+1);
	}printf("\n");
}

Measured memory_bytes: 11648