#include <bits/stdc++.h>
using namespace std;
template <class T>
using VV = vector<vector<T>>;
template <class T>
inline T SMIN(T& a, const T b) {
  return a = (a > b) ? b : a;
}
template <class T>
inline T SMAX(T& a, const T b) {
  return a = (a < b) ? b : a;
}
template <class Iter>
void dumpc(Iter begin, Iter end) {
  if (true) {
    for (; begin != end; ++begin) {
      cerr << *begin << ',';
    }
    cerr << endl;
  }
}
int N;
void yup(int d, int e) {
  fprintf(stdout, "YES %d", d), putchar('\n');
  ;
  fprintf(stdout, "%d", e), putchar('\n');
  ;
}
void nope() {
  fprintf(stdout, "NO"), putchar('\n');
  ;
  exit(0);
}
vector<pair<int, int>> es;
queue<int> Q1;
priority_queue<pair<int, int>> PQ;
int f(int u, int d) {
  if (int((PQ).size()) == 0) {
    bool added = false;
    while (int((Q1).size()) && d--) {
      int v = Q1.front();
      Q1.pop();
      es.emplace_back(u, v);
      added = true;
    }
    if (int((Q1).size())) return -1e8;
    return added;
  } else {
    int u2, d2;
    tie(d2, u2) = PQ.top();
    PQ.pop();
    es.emplace_back(u, u2);
    --d, --d2;
    while (int((Q1).size()) > 1 && d--) {
      int v = Q1.front();
      Q1.pop();
      es.emplace_back(u, v);
    }
    return 1 + f(u2, d2);
  }
}
void solve() {
  cin >> N;
  for (int i = 0; i < (N); ++i) {
    int d;
    cin >> d;
    if (d == 1)
      Q1.emplace(i);
    else
      PQ.emplace(d, i);
  }
  if (int((PQ).size()) == 0) {
    if (int((Q1).size()) == 2) {
      int u = Q1.front();
      Q1.pop();
      int v = Q1.front();
      Q1.pop();
      es.emplace_back(u, v);
      yup(1, int((es).size()));
    } else {
      nope();
    }
  } else {
    int u, d;
    if (int((Q1).size()) == 0) {
      tie(d, u) = PQ.top();
      PQ.pop();
    } else {
      u = Q1.front();
      Q1.pop();
      d = 1;
    }
    int D = f(u, d);
    if (D < 0) nope();
    yup(D, int((es).size()));
  }
  for (auto& e : es)
    fprintf(stdout, "%d %d", e.first + 1, e.second + 1), putchar('\n');
  ;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout << setprecision(12) << fixed;
  solve();
  return 0;
}
