#include <iostream>
#include <iomanip>
#include <random>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <cmath>
#include <numeric>
#include <cstdlib>
#include <cstring>
#include <deque>
#include <sstream>
#include <bitset>
#include <cassert>
#include <fstream>
#include <queue>

#define len(X) ((int)(X).size())

#ifdef __LOCAL
	#define DEBUG_OUTPUT_ENABLED 1
	#define DBG(X) dout << #X << "=" << (X) << '\n';
	#define SAY(X) dout << (X) << '\n';
#else
	#define DEBUG_OUTPUT_ENABLED 0
	#define DBG(X) 42;
	#define SAY(X) 42;
#endif

#define dout   debug::instance

using namespace std;

template<typename T, typename S> inline ostream& operator<<(ostream& os, const pair<T, S> p) { cout << "[" << p.first << ";" << p.second << "]"; return os; }
template<typename T> inline ostream& operator<<(ostream& os, const vector<T>& v) { for(auto el: v) cout << el << " "; return os; }
namespace debug {
	struct DebugStream {
		private:
			bool is_first;
		public:
			DebugStream(bool _is_first): is_first(_is_first) {}
			template<typename T> DebugStream operator<<(const T& value) const {
				assert(DEBUG_OUTPUT_ENABLED);
				if(is_first) cout << "[DBG] ";
				cout << value;
				return DebugStream(false);
			};

			template<typename T> DebugStream printArray(T* l, T* r) {
				assert(DEBUG_OUTPUT_ENABLED);
				if(is_first) cout << "[DBG] ";
				while(l != r) {
					cout << (*l);
					++l;
					if(l == r) {
						cout << '\n';
					} else {
						cout << ' ';
					}
				}
				return DebugStream(false);
			}
	};
	DebugStream instance(true);
};

using ll = long long int;
using ull = unsigned long long int;
using ld  = long double;
using pii = pair<int, int>;
using pll = pair<ll, ll>;

const int INT_INF = (int)(2e9);
const ll  LL_INF  = (ll)(2e18);

const int NIL = -1;
static mt19937 _g(time(nullptr));

inline ll randint(ll a, ll b) { ll w = (_g() << 31LL) ^ _g(); return a + w % (b - a + 1); }
inline void fast_io() { ios_base::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr); };
template<typename T> inline T sign(T x) { return T(x > 0) - T(x < 0); }
template<typename T> inline T fetch() { T ret; cin >> ret; return ret; }
template<typename T> inline vector<T> fetch_vec(int sz) { vector<T> ret(sz); for(auto& elem: ret) cin >> elem; return ret; }
const int MAXN = 507;
int a[MAXN * MAXN], P[MAXN];
vector<pii> Evnt[MAXN * MAXN];
int n;
inline void fail(const string& why) {
	if(DEBUG_OUTPUT_ENABLED) {
		cout << "No: " << why << '\n';
	} else {
		cout << "No\n";
	}
	exit(0);
}
inline void registerEvent(int c, int l, int r) {
	if(l > r) fail("l > r");
	Evnt[l].emplace_back(r, c);
}
int cnt[MAXN];
void solve() {
	cin >> n;
	for(int c = 1; c <= n; ++c) P[c] = fetch<int>() - 1;
	for(int c = 1; c <= n; ++c) {
		int pos = P[c];
		for(int i = 0; i < c - 1; ++i) registerEvent(c, 0, pos - 1);
		registerEvent(c, pos, pos);
		for(int i = 0; i < n - c; ++i) registerEvent(c, pos + 1, n * n - 1);
	}
	multiset<pii> setik;
	for(int i = 0; i < n * n; ++i) {
		DBG(i);
		for(auto& PR: Evnt[i]) setik.insert(PR);
		if(setik.empty()) fail("setik.empty()");
		auto iter = setik.begin();
		if(iter->first < i) fail("iter->first > i");
		a[i] = iter->second;
		setik.erase(iter);
		DBG(a[i]);
	}
	if(!setik.empty()) fail("!setik.empty()");
	for(int i = 0; i < n * n; ++i) if(!a[i]) fail("!a[i]");
	cout << "Yes\n";
	for(int i = 0; i < n * n; ++i) cout << a[i] << ' ';
	cout << '\n';
}

int main() {
	fast_io();
	solve();

	return 0;
}


Measured memory_bytes: 19316