#include<bits/stdc++.h>
#define rp(i,a,b) for(int i=a;i<=b;++i)
#define pr(i,a,b) for(int i=a;i>=b;--i)
#define pb push_back
#define ls (x<<1)
#define rs (ls|1)
#define md ((l+r)>>1)
using namespace std;
const int N=5e3+5,M=1<<14|5;
char s[N];
int T,n,a[N][N],b[N][N],mx[M],tag[M],posn,pos[N];
vector<int> lft[N];
inline void pushup(int x){
	mx[x]=max(mx[ls],mx[rs]);
}
inline void addtag(int x,int y){
	mx[x]+=y;
	tag[x]+=y;
}
inline void pushdown(int x){
	addtag(ls,tag[x]),addtag(rs,tag[x]),tag[x]=0;
}
void build(int x,int l,int r){
	mx[x]=r-1;
	tag[x]=0;
	if(l!=r){
		build(ls,l,md),build(rs,md+1,r);
	}
}
void modify(int x,int l,int r,int lx,int rx){
	if(l>=lx&&r<=rx){
		addtag(x,1);
		return;
	}
	pushdown(x);
	if(lx<=md)modify(ls,l,md,lx,rx);
	if(rx>md)modify(rs,md+1,r,lx,rx);
	pushup(x);
}
int query(int x,int l,int r,int lx,int rx){
	if(l>=lx&&r<=rx){
		return mx[x];
	}
	pushdown(x);
	int ret=0;
	if(lx<=md)ret=max(ret,query(ls,l,md,lx,rx));
	if(rx>md)ret=max(ret,query(rs,md+1,r,lx,rx));
	return ret;
}
void mian(){
	scanf("%d",&n);
	int zero=0;
	rp(i,1,n){
		scanf("%s",s+1);
		rp(j,1,i){
			a[i][j]=s[j]&1;
			b[i][j]=0;
			if(!a[i][j]){
				++zero;
			}
		}
	}
	if(zero!=n){
		puts("Impossible!");
		return;
	}
	rp(i,1,n){
		lft[i].clear();
	}
	rp(i,1,n){
		rp(j,i,n){
			if(!a[j][j+1-i]){
				lft[i].pb(j+1-i);
			}
		}
	}
	build(1,1,n);
	rp(i,1,n){
		for(int l:lft[n+1-i]){
			// printf("[%d %d]\n",l,i);
			modify(1,1,n,1,l);
		}
		if(query(1,1,n,1,i)>i){
			puts("Impossible!");
			return;
		}
	}
	// puts("ok");
	pr(i,n,1){
		if(i!=n){
			posn=0;
			rp(j,1,i+1)if(!a[i+1][j]||b[i+1][j]==2){
				pos[++posn]=j;
			}
			if(!posn){
				puts("Impossible!");
				return;
			}
			// printf("row=%d\n",i);
			// rp(j,1,posn){
			// 	printf("%d:%d\n",j,pos[j]);
			// }
			rp(j,1,pos[1]-1){
				b[i+1][j]=3;
			}
			rp(j,pos[posn]+1,i+1){
				b[i+1][j]=1;
			}
			build(1,1,i);
			int posp=0,low=0;
			// printf("solve(%d)\n",i);
			// _sleep(100);
			rp(j,1,i){
				if(j==pos[1]){
					posp=2;
					low=pos[1];
				}
				for(int l:lft[i+1-j]){
					// printf("[%d %d]\n",l,j);
					modify(1,1,i,1,l);
				}
				if(query(1,1,i,1,j)>j){
					puts("Impossible!");
					return;
				}
				if(low&&query(1,1,i,1,low)==j){
					low=j+1;
				}
				// printf("low=%d\n",low);
				if(posp<=posn&&j+1==pos[posp]){
					if(low==j+1){
						puts("Impossible!");
						return;
					}
					// puts("before...");
					// rp(ii,1,n)rp(jj,1,ii){
					// 	printf("%d",b[ii][jj]);
					// 	if(jj==ii)puts("");
					// }
					rp(k,pos[posp-1]+1,low){
						b[i+1][k]=1;
					}
					b[i][low]=2;
					rp(k,low+1,pos[posp]-1){
						b[i+1][k]=3;
					}
					// printf("omg low=%d\n",low);
					// rp(ii,1,n)rp(jj,1,ii){
					// 	printf("%d",b[ii][jj]);
					// 	if(jj==ii)puts("");
					// }
					modify(1,1,i,1,low);
					low=pos[posp++];
				}
			}
		}
		rp(j,1,i)if(!a[i][i+1-j]){
			lft[j].pop_back();
		}
	}
	rp(i,1,n){
		rp(j,1,i)s[j]="-123"[b[i][j]];
		s[i+1]=0;
		puts(s+1);
	}
}
int main(){

	//freopen("1.in","r",stdin);

	scanf("%d",&T);
	while(T--){
		mian();
	}
	return 0;
}