#include<bits/stdc++.h> 
#define L(i, j, k) for(int i = (j); i <= (k); ++i)
#define R(i, j, k) for(int i = (j); i >= (k); --i)
#define ll long long 
#define vi vector < int > 
#define sz(a) (int) (a).size()
#define me(a, x) memset(a, x, sizeof(a))
using namespace std;
template < int N, int Ne >  struct flows { 
	using F = int; // flow type
	F inf = 1e9;
	int n, s, t; // Remember to assign n, s and t !
	int ehd[N], cur[N], ev[Ne << 1], enx[Ne << 1], eid = 1;
	void clear() {
		eid = 1;
		L(i, 1, n) ehd[i] = 0;
	}
	F ew[Ne << 1], dis[N];
	void Eadd(int u, int v, F w) {
		++eid, enx[eid] = ehd[u], ew[eid] = w, ev[eid] = v, ehd[u] = eid;
	}
	void add(int u, int v, F w) {
		Eadd(u, v, w), Eadd(v, u, 0);
	}
	int que[N], ql, qr;
	bool bfs() {
		ql = 1, qr = 0;
		L(i, 1, n) dis[i] = inf, cur[i] = ehd[i]; 
		que[++qr] = s, dis[s] = 0;
		while(ql <= qr) {
			int u = que[ql++];
			for(int i = ehd[u]; i; i = enx[i]) if(ew[i] && dis[ev[i]] == inf) {
				dis[ev[i]] = dis[u] + 1, que[++qr] = ev[i];
			}
		}
		return dis[t] < inf;
	}
	F dfs(int x, F now) {
		if(!now || x == t) return now;
		F res = 0, f;
		for(int i = cur[x]; i; i = enx[i]) {
			cur[x] = i;
			if(ew[i] && dis[ev[i]] == dis[x] + 1) {
				f = dfs(ev[i], min(now, ew[i])), ew[i] -= f, now -= f, ew[i ^ 1] += f, res += f;
				if(!now) break;
			}
		}
		return res;
	}
	F max_flow() {
		F res = 0;
		while(bfs()) res += dfs(s, inf);
		return res;
	}
} ; 

const int N = 5007;

flows < N * N, N * N * 3 / 2 > G;

int n;
char s[N][N];
int Pl[N][N], Pr[N][N];
bool win[N][N], gol[N][N];
short ans[N][N];

void Main() {
	G.clear();
	cin >> n;
	L(i, 1, n) {
		cin >> (s[i] + 1); 
	}
	L(i, 0, n + 1) {
		L(j, 0, i + 1) {
			Pl[i][j] = 0;
			Pr[i][j] = 0;
			win[i][j] = 0;
			gol[i][j] = 0;
			ans[i][j] = 0;
		}
	}
	G.n = 0;
	L(i, 1, n) {
		L(j, 1, i) {
			Pl[i][j] = ++G.n; 
			Pr[i][j] = ++G.n;
		}
	}
	G.s = ++G.n, G.t = ++G.n;
	L(i, 1, n) 
		Pr[i][0] = G.t;
	L(i, 1, n) 
		L(j, 1, i) 
			if(s[i][j] == '0') 
				G.add(G.s, Pr[i][j], 1);
	int cur = G.eid;
	L(i, 1, n) {
		L(j, 1, i) {
			G.add(Pr[i][j], Pl[i][j], 1);
			G.add(Pl[i][j], Pr[i][j - 1], 1);
			if(i > 1 && j > 1) G.add(Pl[i][j], Pr[i - 1][j - 1], 1);
		}
	} 
	if(G.max_flow() < n) {
		cout << "Impossible!\n";
		return ;
	}
	L(i, 1, n) {
		L(j, 1, i) {
			win[i][j] = 0;
			if(G.ew[cur + 2]) {
				win[i][j] = 1;
			} 
			cur += 2;
			if(win[i][j]) {
				if(G.ew[cur + 2]) {
					gol[i][j] = true;
				} else {
					gol[i][j] = false;
				}
			}
			cur += 2;
			if(i > 1 && j > 1) cur += 2;
		}
	}
	L(i, 1, n) {
		L(j, 1, i) {
			if(s[i][j] != '0') {
				ans[i][j] = 1;
			}
		}
	}
	L(i, 1, n) {
		L(j, 1, i) {
			if(s[i][j] == '0') {
				int x = i, y = j;
				while(y > 0) {
					if(gol[x][y]) {
						if(y > 1) ans[x][y - 1] = 3;
						--y;
					} else {
						int cy = y - 1;
						while(!win[x][cy]) {
//							cout << x << ' ' << y << " : " << x << " and " << cy << endl;
							if(!cy) {
								assert(false);
//								cout << "jinhya\n";
								continue;
							}
							ans[x][cy] = 1, --cy;
						}
						ans[x - 1][y - 1] = 2;
						--x, --y;
					}
				}
			}
		}
	}
	L(i, 1, n) {
		L(j, 1, i) {
			if(s[i][j] == '0') cout << "-";
			else cout << ans[i][j];
		}
		cout << '\n';
	}
}
int main() {
	ios :: sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	int t; cin >> t; while(t--) Main();
	return 0;
}