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

#include<cstdio>
#include<cassert>
#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 ts[MAX_N], 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; i++) {
  //  printf("%d:", (int)cvs[i].size());
  //  for (auto u: cvs[i]) printf(" %d", u);
  //  putchar('\n');
  //}

  fill(rs, rs + n, -1);
  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++, u = (u + 2) % lj)
	  ts[u] = cv[k];
	for (int k = 0; k < lj; k++)
	  rs[ts[k]] = ts[(k + 1) % lj];
      }
    }
    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; k < lj; k++) {
	  rs[cv0[k]] = cv1[k];
	  rs[cv1[k]] = cv0[(k + 1) % lj];
	}
      }
    }
  }
  for (int i = 0; i < n; i++) assert(rs[i] >= 0);

  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;
}