#ifdef LOCAL
#include <header.hpp>
#else
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/priority_queue.hpp>
#endif
#define lp(i, l, r) for (int i = int(l); i <= int(r); ++i)
#define rp(i, r, l) for (int i = int(r); i >= int(l); --i)
#define ls(i, s) for (int i = 0; s[i]; ++i)
#define rs(i, s) for (int i = strlen(s) - 1; i >= 0; --i)
#define lv(i, a) for (int i = 0; i < int(a.size()); ++i)
#define rv(i, a) for (int i = int(a.size()) - 1; i >= 0; --i)
#define lm(i, m) for (auto i = m; i; i = (m) & (i - 1))
#define fast ios::sync_with_stdio(0), cin.tie(0)
#define prec(a) cout << setprecision(a) << fixed
#define freo freopen
#define many for (int C = 1, T = -1; T < 0 && cin >> T, C <= T; ++C)
#define goog "Case #" + to_string(C) + ":"
#define appl(x, f, a)                                                          \
    for (auto &x : a)                                                          \
    f
#define rang(a) begin(a), end(a)
#define gene template <class T>
#define elif else if
#define retu return
#define pufr push_front
#define pofr pop_front
#define puba push_back
#define poba pop_back
#define intm modular_arithmetic::integer
#ifdef LOCAL
#define debu(...)                                                              \
    {                                                                          \
        vector<string> S(1);                                                   \
        int A = 0;                                                             \
        for (auto V : #__VA_ARGS__) {                                          \
            if (V == ',' && !A)                                                \
                S.push_back("");                                               \
            elif (S.back().size() || V != ' ') S.back().push_back(V),          \
                A += map<char, int>{{'(', 1}, {')', -1}}[V];                   \
        }                                                                      \
        cout << "\x1b[33mDebugging @ Line " << __LINE__ << ":\x1b[0m\n";       \
        debu_(S.begin(), __VA_ARGS__);                                         \
    }
#else
#define endl '\n'
#undef assert
#define assert
#define debu(...)
#endif
using namespace std;
using namespace __gnu_pbds;
template <class I, class T> void debu_(I i, const T &a) {
    cout << "\x1b[33m    " << *i << " = " << a << "\x1b[0m" << endl;
}
template <class I, class A, class... B>
void debu_(I i, const A &a, const B &... b) {
    cout << "\x1b[33m    " << *i << " = " << a << "\x1b[0m" << endl;
    debu_(++i, b...);
}
gene int barr(const T &a, const int &i) { return a >> i & 1; }
gene int bcnt(const T &a) { return __builtin_popcount(a); }
template <> int bcnt(const int64_t &a) { return __builtin_popcountll(a); }
gene int blen(const T &a) { return a ? 1 + blen(a / 2) : 0; }
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
gene T rint(T l, T r) { return uniform_int_distribution<T>(l, r)(rng); }
template <class A, class B> B conv(const A &a, B b = B()) {
    stringstream s;
    s << a, s >> b;
    return b;
}
vector<string> spli(const string &s) {
    vector<string> r;
    for (int i = 0, j; i < int(s.size()); ++i)
        if (s[i] != ' ') {
            r.push_back("");
            for (j = i; j < int(s.size()) && s[j] != ' '; ++j)
                r.back().push_back(s[j]);
            i = j - 1;
        }
    return r;
}
template <class T = int> T inf() { return numeric_limits<T>::max() / 2; }
gene T sign(const T &a) { return a == 0 ? 0 : (a < 0 ? -1 : 1); }
gene T floor(const T &a, const T &b) {
    assert(b != 0);
    return sign(a) * sign(b) > 0 ? abs(a) / abs(b)
                                 : -(abs(a) + abs(b) - 1) / abs(b);
}
gene T ceil(const T &a, const T &b) {
    assert(b != 0);
    return sign(a) * sign(b) > 0 ? (abs(a) + abs(b) - 1) / abs(b)
                                 : -abs(a) / abs(b);
}
gene bool tmin(T &a, T b) { return b < a ? a = b, true : false; }
gene bool tmax(T &a, T b) { return b > a ? a = b, true : false; }
#if __cplusplus > 201103L
gene auto min(const T &a) { return *min_element(rang(a)); }
gene auto max(const T &a) { return *max_element(rang(a)); }
gene auto sum(const T &a) {
    return accumulate(rang(a), (typename T::value_type)0);
}
template <> auto sum(const vector<int> &a) { return accumulate(rang(a), 0ll); }
template <> auto sum(const vector<string> &a) {
    return accumulate(rang(a), string());
}
template <class T, int N, class I>
auto gath(const T (&a)[N], const vector<I> &b) {
    vector<T> r;
    for (auto i : b)
        r.push_back(a[i]);
    return r;
}
template <class T, class I> auto gath(const T &a, const vector<I> &b) {
    vector<typename T::value_type> r;
    for (auto i : b)
        r.push_back(a[i]);
    return r;
}
#endif
template <class F, class S> pair<F, S> operator-(pair<F, S> a) {
    a.first = -a.first, a.second = -a.second;
    return a;
}
template <class F, class S>
pair<F, S> operator+(pair<F, S> a, const pair<F, S> &b) {
    a.first += b.first, a.second += b.second;
    return a;
}
template <class F, class S>
pair<F, S> operator-(pair<F, S> a, const pair<F, S> &b) {
    a.first -= b.first, a.second -= b.second;
    return a;
}
template <class F, class S>
pair<F, S> operator*(pair<F, S> a, const pair<F, S> &b) {
    a.first *= b.first, a.second *= b.second;
    return a;
}
template <class F, class S>
pair<F, S> operator/(pair<F, S> a, const pair<F, S> &b) {
    a.first /= b.first, a.second /= b.second;
    return a;
}
template <class F, class S, class T>
pair<F, S> operator*=(pair<F, S> a, const T &b) {
    a.first *= b, a.second *= b;
    return a;
}
template <class F, class S, class T>
pair<F, S> operator/=(pair<F, S> a, const T &b) {
    a.first /= b, a.second /= b;
    return a;
}
gene vector<T> &operator--(vector<T> &a) {
    for (auto &v : a)
        --v;
    return a;
}
template <class T, class C = less<typename T::value_type>>
void sort(T &a, C c = C()) {
    sort(rang(a), c);
}
gene void reve(T &a) { reverse(rang(a)); }
gene void uniq(T &a) { sort(a), a.erase(unique(rang(a)), end(a)); }
gene void shuf(T &a) { shuffle(rang(a), rng); }
gene void merg(T l, T m, T r) { inplace_merge(l, m, r); }
#if __cplusplus > 201103L
gene auto vect(const vector<T> &v, int n) {
    assert(v.size() <= 1);
    return v.size() ? vector<T>(n, v[0]) : vector<T>(n);
}
template <class T, class... D> auto vect(const vector<T> &v, int n, D... m) {
    return vector<decltype(vect(v, m...))>(n, vect(v, m...));
}
#endif
#if __cplusplus > 201103L
template <class F, class S> struct zip_ {
    zip_(F &f_, S &s_) : f(f_), s(s_) {}
    template <class I, class J> struct iterator {
        iterator(I i_, J j_) : i(i_), j(j_) {}
        auto operator++() {
            ++i, ++j;
            return *this;
        }
        bool operator!=(const iterator &a) const {
            assert((i != a.i && j != a.j) || (i == a.i && j == a.j));
            return i != a.i;
        }
        auto operator*() const { return make_pair(*i, *j); }
        I i;
        J j;
    };
    template <class I, class J> auto iterator_(I i, J j) const {
        return iterator<I, J>(i, j);
    }
    auto begin() const { return iterator_(f.begin(), s.begin()); }
    auto end() const { return iterator_(f.end(), s.end()); }
    F &f;
    S &s;
};
template <class F, class S> auto zip(F &f, S &s) { return zip_<F, S>(f, s); }
#endif
template <class F, class S> istream &operator>>(istream &s, pair<F, S> &a) {
    return s >> a.first >> a.second;
}
template <class F, class S>
ostream &operator<<(ostream &s, const pair<F, S> &a) {
    return s << a.first << " " << a.second;
}
gene istream &operator>>(istream &s, vector<T> &a) {
    for (int i = -1; i < 0 && (!a.size() && (cin >> i, a.resize(i), 0), i = 0),
             i < int(a.size());)
        s >> a[i++];
    return s;
}
gene ostream &operator<<(ostream &s, const vector<T> &a) {
    lv(i, a) cout << a[i] << (i + 1 == int(a.size()) ? "" : " ");
    return s;
}
gene T in() {
    T a;
    cin >> a;
    return a;
}
gene void in(T &a) { cin >> a; }
template <class A, class... B> void in(A &a, B &... b) { cin >> a, in(b...); }
void ou() { cout << endl; }
gene void ou(const T &a) { cout << a << endl; }
template <class A, class... B> void ou(const A &a, const B &... b) {
    cout << a << ' ', ou(b...);
}
gene struct hash_safe {};
template <> struct hash_safe<int> {
    size_t operator()(uint64_t a) const {
        static uint64_t d =
            chrono::steady_clock::now().time_since_epoch().count();
        a += d + 0x9e3779b97f4a7c15, a = (a ^ (a >> 30)) * 0xbf58476d1ce4e5b9;
        a = (a ^ (a >> 27)) * 0x94d049bb133111eb;
        return a ^ (a >> 31);
    }
};
template <> struct hash_safe<int64_t> {
    size_t operator()(uint64_t a) const {
        static uint64_t d =
            chrono::steady_clock::now().time_since_epoch().count();
        a += d + 0x9e3779b97f4a7c15, a = (a ^ (a >> 30)) * 0xbf58476d1ce4e5b9;
        a = (a ^ (a >> 27)) * 0x94d049bb133111eb;
        return a ^ (a >> 31);
    }
};
gene using unordered_set_safe = unordered_set<T, hash_safe<T>>;
gene using unordered_multiset_safe = unordered_multiset<T, hash_safe<T>>;
template <class K, class V>
using unordered_map_safe = unordered_map<K, V, hash_safe<K>>;
template <class K, class V>
using unordered_multimap_safe = unordered_multimap<K, V, hash_safe<K>>;
#define unordered_set unordered_set_safe
#define unordered_multiset unordered_multiset_safe
#define unordered_map unordered_map_safe
#define unordered_multimap unordered_multimap_safe
gene using bbst =
    tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
gene using heap = __gnu_pbds::priority_queue<T, greater<T>, pairing_heap_tag>;
typedef int64_t ll;
typedef long double ld;
typedef string str;
typedef pair<int, int> pii;
typedef pair<int, ll> pil;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<pii> vpii;
typedef vector<pil> vpil;
typedef vector<pli> vpli;
typedef vector<pll> vpll;
typedef vector<vi> vvi;
typedef vector<vl> vvl;
namespace maximum_flow {
template <class T> struct maximum_flow {
    struct edge {
        int v;
        T c, l;
        edge(int _v, T _c) : v(_v), c(_c), l(_c) {}
    };
    maximum_flow(int _n, int _src, int _snk)
        : n(_n), src(_src), snk(_snk), bge(n), hei(n, _n), gap(n + 1), cur(n),
          frm(n) {}
    int push(int u) {
        for (int i = 0; i < int(bge[u].size()); ++i) {
            edge &e = egs[bge[u][i]];
            if (e.l && e.c != e.l) {
                ++e.c;
                if (e.v == snk)
                    return u;
                return push(e.v);
            }
        }
    }
    void label() {
        hei[snk] = 0;
        queue<int> qu;
        qu.push(snk);
        for (int u; qu.empty() ? 0 : (u = qu.front(), qu.pop(), 1);)
            for (int i = 0; i < int(bge[u].size()); ++i) {
                edge &e = egs[bge[u][i]], &ev = egs[bge[u][i] ^ 1];
                if (ev.c > 0 && hei[e.v] == n)
                    hei[e.v] = hei[u] + 1, qu.push(e.v);
            }
        for (int i = 0; i < n; ++i)
            ++gap[hei[i]];
    }
    T augment() {
        T f = 0;
        for (int u = snk; u != src; u = egs[frm[u] ^ 1].v)
            if (f <= 0 || f > egs[frm[u]].c)
                f = egs[frm[u]].c;
        for (int u = snk; u != src; u = egs[frm[u] ^ 1].v)
            egs[frm[u]].c -= f, egs[frm[u] ^ 1].c += f;
        return f;
    }
    void add_edge(int u, int v, T c) {
        bge[u].push_back(int(egs.size()));
        egs.push_back(edge(v, c));
        bge[v].push_back(int(egs.size()));
        egs.push_back(edge(u, 0));
    }
    T run() {
        label();
        T r = 0;
        for (int u = src; hei[src] != n;) {
            if (u == snk)
                r += augment(), u = src;
            int f = 0;
            for (int i = cur[u]; i < int(bge[u].size()); ++i) {
                edge &e = egs[bge[u][i]];
                if (e.c > 0 && hei[u] == hei[e.v] + 1) {
                    f = 1;
                    frm[e.v] = bge[u][i];
                    u = e.v;
                    break;
                }
            }
            if (!f) {
                int mh = n - 1;
                for (int i = 0; i < int(bge[u].size()); ++i) {
                    edge &e = egs[bge[u][i]];
                    if (e.c > 0 && mh > hei[e.v])
                        mh = hei[e.v];
                }
                if (!--gap[hei[u]])
                    break;
                ++gap[hei[u] = mh + 1];
                cur[u] = 0;
                if (u != src)
                    u = egs[frm[u] ^ 1].v;
            }
        }
        return r;
    }
    int n, src, snk;
    vector<edge> egs;
    vector<vector<int>> bge;
    vector<int> hei, gap, cur, frm;
};
} // namespace maximum_flow
// The templates end here.
int main() {
    fast, prec(12);
    int n;
    in(n);
    vi x(n);
    in(x);
    vi a(n * n, -1);
    lv(i, x) a[x[i] - 1] = i + 1;
    vector<tuple<int, int, int>> segments;
    lv(i, a) {
        int j = i;
        while (j + 1 < a.size() && a[j + 1] == a[j])
            ++j;
        segments.emplace_back(i, j, a[i]);
        i = j;
    }
    int m = 0;
    int source = m++, sink = m++;
    vi segment_index(segments.size());
    for (auto &v : segment_index)
        v = m++;
    vi prefix_index(segments.size());
    for (auto &v : prefix_index)
        v = m++;
    vi suffix_index(segments.size());
    for (auto &v : suffix_index)
        v = m++;
    vi number_index(n);
    for (auto &v : number_index)
        v = m++;
    maximum_flow::maximum_flow<int> flow(m, source, sink);
    lv(i, segments) {
        flow.add_edge(segment_index[i], sink,
                      get<1>(segments[i]) - get<0>(segments[i]) + 1);
        flow.add_edge(prefix_index[i], segment_index[i], inf());
        flow.add_edge(suffix_index[i], segment_index[i], inf());
        if (i)
            flow.add_edge(prefix_index[i], prefix_index[i - 1], inf());
        if (i + 1 < segments.size())
            flow.add_edge(suffix_index[i], suffix_index[i + 1], inf());
    }
    lp(i, 1, n) {
        flow.add_edge(source, number_index[i - 1], n);
        int j = 0;
        while (get<2>(segments[j]) != i)
            ++j;
        flow.add_edge(number_index[i - 1], segment_index[j], 1);
        if (j)
            flow.add_edge(number_index[i - 1], prefix_index[j - 1], i - 1);
        if (j + 1 < segments.size())
            flow.add_edge(number_index[i - 1], suffix_index[j + 1], n - i);
    }
    if (flow.run() != n * n)
        ou("No");
    else {
        ou("Yes");
        vi current(m);
        lv(i, segments) { current[segment_index[i]] = get<0>(segments[i]); }
        vi answer(n * n);
        lp(i, 1, n) {
            lp(j, 1, n) {
                answer[current[flow.push(number_index[i - 1])]++] = i;
            }
        }
        ou(answer);
    }
    return 0;
}


Measured memory_bytes: 3712