/* -*- coding: utf-8 -*-
 *
 * 3256.cc:  No.3256 Permutation Equation - yukicoder
 */

#include<cstdio>
#include<vector>
#include<algorithm>

using namespace std;

/* constant */

const int MAX_N = 200000;

/* typedef */

using vi = vector<int>;
using vvi = vector<vi>;

/* global variables */

int ps[MAX_N], ss[MAX_N];
bool used[MAX_N];
int rs[MAX_N], qs[MAX_N];

/* subroutines */

bool cmpv(const vi &a, const vi &b) {
  return a.size() < b.size();
}

bool check(int n) {
  for (int i = 0; i < n; i++)
    if (ps[i] != qs[(qs[i] + 1) % n]) { puts("NG"); return false; }
  puts("GOOD");
  return true;
}

/* main */

int main() {
  int n;
  scanf("%d", &n);
  for (int i = 0; i < n; i++) scanf("%d", ps + i), ps[i]--;
  for (int i = 0; i < n; i++) ss[i] = (ps[i] + 1) % n;

  vvi cvs;
  for (int st = 0; st < n; st++)
    if (! used[st]) {
      vi v;
      for (int u = st; ! used[u]; u = ss[u])
	used[u] = true, v.push_back(u);
      cvs.push_back(v);
    }
  sort(cvs.begin(), cvs.end(), cmpv);
  int m = cvs.size();

  for (int i = 0; i < m;) {
    int j = i, lj = cvs[j].size();
    while (i < m && lj == cvs[i].size()) i++;

    if (lj & 1) {
      for (; j < i; j++) {
	auto &cv = cvs[j];
	for (int k = 0, u = 0; k < lj; k++) {
	  int v = (u + 2) % lj;
	  rs[cv[u]] = cv[v];
	  u = v;
	}
      }
    }
    else {
      if ((i - j) & 1) { puts("No"); return 0; }
      for (; j < i; j += 2) {
	auto &cv0 = cvs[j], &cv1 = cvs[j + 1];
	for (int k = 0, u0 = 0, u1 = 0; k < lj * 2; k++) {
	  if (! (k & 1)) { // u0->u1
	    rs[cv0[u0]] = cv1[u1];
	    u0 = (u0 + 1) % lj;
	  }
	  else { // u1->u0
	    rs[cv1[u1]] = cv0[u0];
	    u1 = (u1 + 1) % lj;
	  }
	}
      }
    }
  }

  for (int i = 0; i < n; i++) qs[i] = (rs[i] + n - 1) % n;
  
  puts("Yes");
  for (int i = 0; i < n; i++)
    printf("%d%c", qs[i] + 1, (i + 1 < n) ? ' ' : '\n');

  //check(n);
  return 0;
}