#include <bits/stdc++.h>
using namespace std;
int vis[123456];
int A[123456];
int dfs(int u, int ori) {
  vis[u] = 1;
  if (ori == A[u]) return 1;
  return 1 + dfs(A[u], ori);
}
int main() {
  int n;
  scanf("%d", &n);
  for (int i = 1; i <= n; ++i) {
    scanf("%d", &A[i]);
  }
  int t;
  bool ok = 1;
  bool hasodd = 0;
  int a = -1;
  int pa = -1, pb = -1;
  for (int i = 1; i <= n; ++i) {
    if (!vis[i]) {
      t = dfs(i, i);
      if (t > 2 && (t & 1))
        hasodd = 1;
      else if (t == 2) {
        pa = i;
        pb = A[i];
      } else if (t == 1) {
        a = i;
      }
    }
  }
  if (pa == -1 && a == -1) {
    puts("NO");
  } else if (a == -1 && hasodd) {
    puts("NO");
  } else {
    puts("YES");
    memset(vis, 0, sizeof(vis));
    if (a != -1) {
      for (int i = 1; i <= n; ++i) {
        if (i != a) {
          printf("%d %d\n", a, i);
        }
      }
    } else {
      printf("%d %d\n", pa, pb);
      for (int i = 1; i <= n; ++i) {
        if (pa != i && pb != i) {
          int k1 = i, k2 = A[i];
          while (!vis[k1]) {
            printf("%d %d\n", pa, k1);
            printf("%d %d\n", pb, k2);
            vis[k1] = vis[k2] = 1;
            k1 = A[k2];
            k2 = A[k1];
          }
        }
      }
    }
  }
  return 0;
}
