#include <algorithm>
#include <array>
#include <iostream>
#include <limits>
#include <vector>

using namespace std;

constexpr int N = 4e3;

template<typename T = int>
struct /*Maksim1744*/ Dinic {
	struct Edge {
		int v, u;
		T f, c;
	};
	vector<Edge> ed;
	vector<vector<int> > e;
	int n, st, fin;
	T flow;

	vector<int> ptr, ds, que;

	Dinic(int n_, int st_, int fin_) : n(n_), st(st_), fin(fin_), flow(0) {
		e.resize(n);
		ptr.resize(n);
		ds.resize(n);
		que.resize(n);
	}

	void ae(int v, int u, T c, T bc = 0, T f = 0) {
		int w = int(ed.size());
		ed.push_back({v, u, f, c});
		ed.push_back({u, v, -f, bc});
		e[v].push_back(w);
		e[u].push_back(w + 1);
	}

	bool bfs(T scale) {
		fill(ds.begin(), ds.end(), -1);
		ds[st] = 0;
		int sz = 1;
		que[0] = st;
		for (int it = 0; it < sz; it++) {
			int v = que[it];
			for (int i : e[v])
				if (ed[i].f + scale <= ed[i].c && ds[ed[i].u] == -1) {
					ds[ed[i].u] = ds[v] + 1;
					if (ed[i].u == fin)
						return 1;
					que[sz++] = ed[i].u;
				}
		}
		return 0;
	}

	T dfs(int v, T w, T cs) {
		if (v == fin)
			return w;
		int &i = ptr[v];
		while (i >= 0) {
			auto &t = ed[e[v][i]];
			if (ds[t.u] == ds[v] + 1 && t.f + cs <= t.c) {
				T o = dfs(t.u, min(w, t.c - t.f), cs);
				if (o) {
					t.f += o;
					ed[e[v][i] ^ 1].f -= o;
					return o;
				}
			}
			i--;
		}
		return 0;
	}

	T get(int scale = 30) {
		for (T i = ((T(1)) << scale); i > 0; i >>= 1)
			while (bfs(i)) {
				for (int j = 0; j < n; j++)
					ptr[j] = int(e[j].size()) - 1;
				while (T w = dfs(st, numeric_limits<T>::max(), i))
					flow += w;
			}
		return flow;
	}
};

bool a[N][N], u[N][N], dr[N], dc[N];
int rt[N], ct[N], br[N][2], bc[N][2];

int main() {
  cin.tie(0)->sync_with_stdio(0);
  int n; cin >> n;
  Dinic<> mf(2 * n + 2, 2 * n, 2 * n + 1);
  for (int i = 0; i < n; ++i) if (string s; cin >> s)
  for (int j = 0; j < n; ++j) a[i][j] = s[j] == '+', br[i][1] += s[j] == '+', br[i][0] -= s[j] == '-', bc[j][1] += s[j] == '+', bc[j][0] -= s[j] == '-';
  for (int i = 0; i < n; ++i) cin >> rt[i];
  for (int i = 0; i < n; ++i) cin >> ct[i];
  vector<string> ans(n, string(n, '0'));
  auto mark = [&](int i, int j, bool w) {
    if (u[i][j]) return;
    u[i][j] = 1, ans[i][j] = '0' + (a[i][j] == w);
    br[i][1] -= a[i][j], br[i][0] += !a[i][j];
    bc[j][1] -= a[i][j], bc[j][0] += !a[i][j];
    if (a[i][j] == w) rt[i] -= 2 * a[i][j] - 1, ct[j] -= 2 * a[i][j] - 1;
  };
  for (int ch = 1; ch--; ) {
    for (int i = 0; i < n; ++i) if (!dr[i] && (rt[i] >= br[i][1] || rt[i] <= br[i][0])) if (dr[i] = ch = 1)
    for (int j = 0; j < n; ++j) mark(i, j, rt[i] >= br[i][1]);
    for (int i = 0; i < n; ++i) if (!dc[i] && (ct[i] >= bc[i][1] || ct[i] <= bc[i][0])) if (dc[i] = ch = 1)
    for (int j = 0; j < n; ++j) mark(j, i, ct[i] >= bc[i][1]);
  }
  int s[2]{};
  for (int i = 0; i < n; ++i) if (rt[i] > br[i][1] || rt[i] < br[i][0] || ct[i] > bc[i][1] || ct[i] < bc[i][0]) goto fail;
  for (int i = 0; i < n; ++i)
  for (int j = 0; j < n; ++j) {
    if (!dr[i] && !dc[j] && rt[i] > 0 && ct[j] > 0 && a[i][j]) mark(i, j, 1);
    else if (!dr[i] && !dc[j] && rt[i] < 0 && ct[j] < 0 && !a[i][j]) mark(i, j, 0);
    else if (u[i][j]) ;
    else if (a[i][j]) mf.ae(i, n + j, 1);
    else              mf.ae(n + j, i, 1);
  }
  for (int i = 0; i < n; ++i) if (rt[i] > 0) mf.ae(2 * n, i, rt[i]), s[0] += rt[i]; else mf.ae(i, 2 * n + 1, -rt[i]), s[1] -= rt[i];
  for (int i = 0; i < n; ++i) if (ct[i] < 0) mf.ae(2 * n, n + i, -ct[i]), s[0] -= ct[i]; else mf.ae(n + i, 2 * n + 1, ct[i]), s[1] += ct[i];
  if (s[0] != s[1] || s[0] != mf.get(1)) fail: return cout << "NO\n", 0;
  cout << "Yes\n";
  for (int i = 0; i < mf.ed.size(); i += 2) if (auto& e = mf.ed[i]; max(e.u, e.v) < 2 * n && e.f) {
    ans[min(e.u, e.v) % n][max(e.u, e.v) % n] = '1';
  }
  for (auto& s: ans) cout << s << '\n';
}