// Slower C++ solution:
#define _CRT_SECURE_NO_WARNINGS

#include <iostream>

#include <sstream>

#include <iomanip>

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cctype>

#include <cstring>

#include <vector>

#include <list>

#include <queue>

#include <deque>

#include <stack>

#include <map>

#include <set>

#include <algorithm>

#include <iterator>

#include <bitset>

#include <ctime>

#include<complex>

using namespace std;



#define FOR(i,a,b) for (int i = (a); i < (b); i++)

#define RFOR(i,b,a) for (int i = (b)-1; i >= (a); i--)

#define ITER(it,a) for (__typeof(a.begin()) it = a.begin(); it != a.end(); it++)

#define FILL(a,value) memset(a, value, sizeof(a))



#define SZ(a) (int)a.size()

#define ALL(a) a.begin(), a.end()

#define PB push_back

#define MP make_pair



typedef long long LL;

typedef vector<int> VI;

typedef pair<int, int> PII;



const double PI = acos(-1.0);

const int INF = 1000 * 1000 * 1000 + 7;

const LL LINF = INF * (LL)INF;



const int MAX = 1000 * 1000 + 47;

struct Item

{

	int i, pos, left;

	Item(int i1, int pos1, int left1)

	{

		i = i1;

		pos = pos1;

		left = left1;

	}



	Item() {}

};



int A[MAX];

Item X[MAX];



bool cmp(Item a, Item b)

{

	return a.pos < b.pos;

}



void no()

{

	cout << "No" << endl;

	exit(0);

}



int main()

{

	//freopen("in.txt", "r", stdin);

	//freopen("out.txt", "w", stdout);

	ios::sync_with_stdio(false); cin.tie(0);

	int n;

	cin >> n;

	FOR(i, 0, n)

	{

		int x;

		cin >> x;

		x--;

		X[i] = Item(i, x, i);

	}



	sort(X, X + n, cmp);

	FILL(A, -1);

	int p = 0;

	FOR(i, 0, n)

	{

		//cout << X[i].i << " " << X[i].left << endl;

		while (X[i].left != 0)

		{

			if (A[p] == -1)

			{

				X[i].left--;

				A[p] = X[i].i;

			}



			p++;

		}



		if (A[X[i].pos] != -1) no();

		A[X[i].pos] = X[i].i;

	}



	p = 0;

	FOR(i, 0, n * n)

	{

		while (p < n && X[p].left == n - X[p].i - 1) p++;

		if (p == n) break;



		if (A[i] == -1)

		{

			if (X[p].pos > i) no();

			A[i] = X[p].i;

			X[p].left++;

		}

	}



	//FOR(i, 0, n)

	//{

	//	FOR(k, 0, n - i - 1)

	//	{

	//		while (A[p] != -1) p++;

	//		A[p] = X[k].i;

	//	}

	//}



	cout << "Yes" << endl;

	FOR(i, 0, n * n)

	{

		cout << A[i] + 1 << " ";

	}



	cout << endl;

}


// Faster C++ solution for the same problem (measured speedup: 5.031068358797098x):
#include <bits/stdc++.h>



using namespace std;



typedef long long ll;

typedef long double ld;



#define mp make_pair

#define ff first

#define ss second

#define pb push_back

#define pf push_front

#define popb pop_back

#define popf pop_front

#define all(v) v.begin(),v.end()



const int maxn = 5e5 + 100;

const ll inf = 2e18, mod = 1e9 + 7;



int n, x[maxn], a[maxn];

bool alpha[maxn], valid = true;



int main()

{

	ios_base::sync_with_stdio(false);

	cin.tie(0);



	cin >> n;

	for (int i = 0; i < n; i++)

		cin >> x[i], x[i]--;



	for (int i = 0; i < n; i++)

	{

		if (a[x[i]])

			valid = false;

		a[x[i]] = i + 1;

		alpha[x[i]] = true;

	}

	int pos = 0;

	for (int i = 0; i < n * n; i++)

	{

		if (a[i])

		{

			int cnt = a[i] - 1;

			while (cnt && pos <= i)

			{

				if (a[pos])

					pos++;

				else

					a[pos] = a[i], cnt--, pos++;

			}

			if (cnt)

				valid = false;

		}

	}



	pos = n * n - 1;

	for (int i = n * n - 1; i >= 0; i--)

	{

		if (alpha[i])

		{

			int cnt = n - a[i];

			while (cnt && pos >= i)

			{

				if (a[pos])

					pos--;

				else

					a[pos] = a[i], cnt--, pos--;

			}

			if (cnt)

				valid = false;

		}

	}



	if (valid)

	{

		cout << "Yes" << endl;

		for (int i = 0; i < n * n; i++)

			cout << a[i] << ' ';

		cout << endl;

	}

	else

		cout << "No" << endl;

}