#include<bits/stdc++.h>
using namespace std;
const int MAXN=5005,N=1<<13;
int tr[N<<1];
void psu(int p) {
	int t=max(tr[p<<1],tr[p<<1|1]);
	tr[p]+=t,tr[p<<1]-=t,tr[p<<1|1]-=t;
}
void build(int n) {
	memset(tr,0,sizeof(tr));
	for(int i=1;i<=n;++i) tr[i+N]=i-1;
	for(int i=N-1;i;--i) psu(i);
}
void add(int l,int r) {
	for(l+=N-1,r+=N+1;l^r^1;l>>=1,r>>=1) {
		if(~l&1) ++tr[l^1];
		if(r&1) ++tr[r^1];
		psu(l>>1),psu(r>>1);
	}
	while(l>1) psu(l>>=1);
}
int qry(int l,int r) {
	int sl=0,sr=0;
	for(l+=N,r+=N;(l^r)>1;l>>=1,r>>=1) {
		sl+=tr[l],sr+=tr[r];
		if(~l&1) sl=max(sl,tr[l^1]);
		if(r&1) sr=max(sr,tr[r^1]);
	}
	int s=max(sl+tr[l],sr+tr[r]);
	while(l>1) s+=tr[l>>=1];
	return s;
}
char s[MAXN][MAXN],t[MAXN][MAXN];
vector <int> g[MAXN],o[MAXN];
void init(int x) {
	build(x);
	for(int i=1;i<=x;++i) o[i].clear();
	for(int i=1;i<=x;++i) for(int j:g[i]) {
		o[j+x-i].push_back(j);
	}
}
int n;
void solve() {
	scanf("%d",&n);
	for(int i=1;i<=n;++i) {
		scanf("%s",s[i]+1);
		for(int j=1;j<=i;++j) if(s[i][j]=='0') g[i].push_back(j),t[i][j]='-';
	}
	init(n);
	for(int i=1;i<=n;++i) {
		for(int j:o[i]) add(1,j);
		if(qry(1,i)>i) return puts("Impossible!"),void();;
	}
	for(int i=n;i>1;--i) {
		if(g[i].empty()) return puts("Impossible!"),void();;
		init(i-1);
		int lst=0,pos=0;
		for(int j=1;j<=i;++j) {
			for(int k:o[j]) add(1,k);
			if(s[i][j]=='0') {
				s[i-1][pos]='0',t[i-1][pos]='2',g[i-1].push_back(pos);
				for(int k=lst+1;k<=pos;++k) t[i][k]='1';
				for(int k=pos+1;k<j;++k) t[i][k]='3';
				add(1,pos),lst=pos=j;
			}
			if(j==i) break;
			if(qry(1,j)>j) return puts("Impossible!"),void();;
			if(lst&&qry(1,lst)>=j) pos=j+1;
		}
		for(int j=lst+1;j<=i;++j) t[i][j]='1';
	}
	for(int i=1;i<=n;++i) printf("%s\n",t[i]+1);
}
signed main() {
	int cas; scanf("%d",&cas);
	while(cas--) {
		solve();
		for(int i=1;i<=n;++i) {
			g[i].clear(),o[i].clear();
			for(int j=1;j<=i;++j) s[i][j]=t[i][j]=0;
		}
	}
}