#include <bits/stdc++.h>
using namespace std;
const int maxn = 100010;
int a[maxn];
int n;
void init() {
  scanf("%d", &n);
  for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
}
int vis[maxn];
int tot;
vector<int> C[maxn];
void dfs(int x) {
  if (vis[x]) {
    return;
  }
  vis[x] = 1;
  C[tot].push_back(x);
  dfs(a[x]);
}
void task() {
  init();
  tot = 0;
  for (int i = 1; i <= n; i++)
    if (!vis[i]) {
      tot++;
      dfs(i);
    }
  for (int i = 1; i <= tot; i++)
    if (C[i].size() == 1) {
      printf("YES\n");
      int u = C[i][0];
      for (int j = 1; j <= n; j++)
        if (j != u) printf("%d %d\n", u, j);
      return;
    }
  bool flag = 0;
  int pos = 0;
  for (int i = 1; i <= tot; i++) {
    if (C[i].size() & 1) {
      flag = 0;
      break;
    }
    if (C[i].size() == 2) {
      flag = 1;
      pos = i;
    }
  }
  if (!flag)
    printf("NO\n");
  else {
    printf("YES\n");
    for (int i = 1; i <= tot; i++) {
      if (i == pos) continue;
      for (int j = 0; j < C[i].size(); j++)
        printf("%d %d\n", C[pos][j & 1], C[i][j]);
    }
    printf("%d %d\n", C[pos][0], C[pos][1]);
  }
}
int main() {
  task();
  return 0;
}
