#include <bits/stdc++.h>
using namespace std;
const int maxn = 25;
long long gcd(long long a, long long b) { return b ? gcd(b, a % b) : a; }
long long qPow(long long a, long long b, long long mod) {
  long long ans = 1, bas = a;
  while (b) {
    if (b & 1) ans = ans * bas % mod;
    bas = bas * bas % mod;
    b >>= 1;
  }
  return ans;
}
map<int, int> mp;
priority_queue<pair<int, int>> q;
int ans[maxn][maxn];
int main() {
  int n, x;
  cin >> n;
  for (int i = 0; i < n * n; i++) {
    cin >> x;
    mp[x]++;
  }
  for (auto p : mp) {
    q.push(make_pair(p.second, p.first));
  }
  for (int i = 0; i < n / 2; i++) {
    for (int j = 0; j < n / 2; j++) {
      auto p = q.top();
      q.pop();
      if (p.first < 4) {
        cout << "NO" << endl;
        return 0;
      }
      ans[i][j] = p.second;
      q.push(make_pair(p.first - 4, p.second));
    }
  }
  if (n % 2) {
    for (int i = 0; i < n / 2; i++) {
      auto p = q.top();
      q.pop();
      if (p.first < 2) {
        cout << "NO" << endl;
        return 0;
      }
      ans[i][n / 2] = p.second;
      q.push(make_pair(p.first - 2, p.second));
    }
    for (int i = 0; i < n / 2; i++) {
      auto p = q.top();
      q.pop();
      if (p.first < 2) {
        cout << "NO" << endl;
        return 0;
      }
      ans[n / 2][i] = p.second;
      q.push(make_pair(p.first - 2, p.second));
    }
    if (!q.top().first) {
      cout << "NO" << endl;
      return 0;
    } else
      ans[n / 2][n / 2] = q.top().second;
  }
  for (int i = 0; i < n / 2; i++)
    for (int j = 0; j < n / 2; j++)
      ans[n - i - 1][j] = ans[i][n - j - 1] = ans[n - i - 1][n - j - 1] =
          ans[i][j];
  if (n % 2) {
    for (int i = 0; i < n / 2; i++) {
      ans[n - i - 1][n / 2] = ans[i][n / 2];
      ans[n / 2][n - i - 1] = ans[n / 2][i];
    }
  }
  cout << "YES" << endl;
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < n; j++) cout << ans[i][j] << " ";
    cout << endl;
  }
  return 0;
}
