#pragma GCC optimize(3,"Ofast","inline")
#include<iostream>
#include<string>
#include<algorithm>
#include<queue>
using namespace std;
#define ll long long
const int N=300005;
const ll mod=998244353;
ll qpow(ll a,ll b){
	ll ans=1;
	if(b==0)
	return 1;
	if(b%2)
	ans=a;
	ll t=qpow(a,b/2);
	return t*t%mod*ans%mod;
}
ll inv(ll a){
	return qpow(a,mod-2);
}
struct BN{
	vector<int>a;
	void resolve(){
		for(int i=0;i<a.size();i++){
			if(a[i]>=10){
				if(i==a.size()-1){
					a.push_back(a[i]/10);
				}
				else{
					a[i+1]+=a[i]/10;
				}
				a[i]%=10;
			}
		}
	}
};
struct pt{
	int L;
	int R;
	int S1;
	int S2;
	int D;
};
struct sv{
	BN num;
	int l;
	int r;
	char c;
};
ll n;
sv maxx(sv &a,sv &b){
	if(a.num.a.size()>b.num.a.size()){
		return a;
	}
	else if(a.num.a.size()<b.num.a.size()){
		return b;
	}
	else{
		for(int i=a.num.a.size()-1;i>=0;i--){
			if(a.num.a[i]<b.num.a[i])return b;
			else if(a.num.a[i]>b.num.a[i])return a;
		}
		return a;
	}
	return a;
}
vector<sv>ans;
string str1,str2;
void rset(int id,char c){
	if(c=='l'){
		if(str1[2*id-2]=='1'||str2[2*id-2]=='2'){
			str1[2*id-2]='1';
			str2[2*id-2]='2';	
		}
		else{
			str1[2*id-2]='2';
			str2[2*id-2]='1';
		}
		str1[2*id-1]='0';
		str2[2*id-1]='0';
	}
	else if(c=='r'){
		if(str1[2*id-1]=='1'||str2[2*id-1]=='2'){
			str1[2*id-1]='1';
			str2[2*id-1]='2';	
		}
		else{
			str1[2*id-1]='2';
			str2[2*id-1]='1';
		}
		str1[2*id-2]='0';
		str2[2*id-2]='0';
	}
	else if(c=='1'){
		if(str2[2*id-2]!='?'){
			str1[2*id-2]='0';
			str2[2*id-2]='1';
		}
		else{
			str1[2*id-2]='1';
			str2[2*id-2]='0';
		}
		if(str2[2*id-1]!='?'){
			str1[2*id-1]='0';
			str2[2*id-1]='2';
		}
		else{
			str1[2*id-1]='2';
			str2[2*id-1]='0';
		}
	}
	else if(c=='2'){
		if(str2[2*id-2]!='?'){
			str1[2*id-2]='0';
			str2[2*id-2]='2';
		}
		else{
			str1[2*id-2]='2';
			str2[2*id-2]='0';
		}
		if(str2[2*id-1]!='?'){
			str1[2*id-1]='0';
			str2[2*id-1]='1';
		}
		else{
			str1[2*id-1]='1';
			str2[2*id-1]='0';
		}
	}
}
void check(char c,int l,int r,int ext){
	BN num;
	vector<int>g;
	vector<int>p(2*(r-l+1)+2);
	for(int i=l;i<=r;i++){
		if(str2[2*i-2]!='?'){
			g.push_back(3);
		}
		else{
			if(c=='1')
			g.push_back(1);
			else
			g.push_back(2);
		}
		if(str2[2*i-1]!='?'){
			g.push_back(3);
		}
		else{
			if(c=='1')
			g.push_back(2);
			else
			g.push_back(1);
		}
	}
	int cur=0;
	for(int i=g.size()-1;i>=0;i--){
		if(g[i]!=3){
			p[cur]+=g[i];
			cur++;
		}
		else{
			cur=0;
		}
	}
	p[0]+=ext;
	for(int i=p.size()-1;i>=0;i--){
		if(p[i]!=0){
			p.resize(i+1);
			break;
		}
		if(i==0){
			p.resize(1);
		}
	}
	num.a=p;
	ans.push_back({num,l,r,c});
}
void rop(sv x){
	//cout<<x.l<<' '<<x.r<<' '<<x.c<<'\n';
	
	for(int i=x.num.a.size()-1;i>=0;i--){
		cout<<x.num.a[i];
	}
	cout<<'\n';
	if(x.c=='?'){
		for(int i=1;i<=n;i++){
			if(i<=x.r)
			rset(i,'l');
			else
			rset(i,'r');
		}
	}
	else{
		for(int i=1;i<x.l;i++){
			rset(i,'l');
		}
		for(int i=x.l;i<=x.r;i++){
			rset(i,x.c);
		}
		for(int i=x.r+1;i<=n;i++){
			rset(i,'r');
		}
	}
	cout<<str1<<'\n';
	cout<<str2<<'\n';
}
void op(int x){
	if(ans.size()==0){
		cout<<"impossible\n";
		return;
	}
	for(int i=0;i<ans.size();i++){
		ans[i].num.resolve();
	}
	sv mx=ans[0];
	for(int i=1;i<ans.size();i++){
		mx=maxx(ans[i],mx);
	}
	rop(mx);
	//sort(ans.begin(),ans.end(),cmp);
}
void solve(){
	
	cin>>n;
	
	cin>>str1>>str2;
	vector<pt>a(n+1);
	ans.clear();
	vector<int>sl(n+1),sr(n+1),sd(n+1);
	for(int i=1;i<=n;i++){
		a[i]={0,0,0,0,0};
		vector<char>cur;
		cur.push_back(str1[2*i-2]);
		cur.push_back(str1[2*i-1]);
		cur.push_back(str2[2*i-2]);
		cur.push_back(str2[2*i-1]);
		if(cur[0]!='?'&&cur[2]!='?'){
			if(cur[1]!='?'||cur[3]!='?'){
				op(0);
				return ;
			}
			if(cur[0]==cur[2]){
				op(0);
				return ;
			}
			a[i].L=1;
			a[i].D=cur[0]-'0';
		}
		else if(cur[1]!='?'&&cur[3]!='?'){
			if(cur[0]!='?'||cur[2]!='?'){
				op(0);
				return ;
			}
			if(cur[1]==cur[3]){
				op(0);
				return ;
			}
			a[i].R=1;
			a[i].D=cur[1]-'0';
		}
		else{
			char l=min(cur[0],cur[2]);
			char r=min(cur[1],cur[3]);
			if(l=='?'&&r=='?'){
				a[i]={1,1,1,1,2};
			}
			else if(l==r){
				op(0);
				return;
			}
			else if(l=='?'){
				a[i].R=1;
				if(r=='1'){
					a[i].S2=1;
				}
				else if(r=='2'){
					a[i].S1=1;
				}
				
				if(cur[1]=='1'||cur[3]=='2'){
					a[i].D=1;
					
				}
				else{
					a[i].D=2;
				}
					
			}
			else if(r=='?'){
				a[i].L=1;
				if(l=='1'){
					a[i].S1=1;
				}
				else if(l=='2'){
					a[i].S2=1;
				}
				if(cur[0]=='1'||cur[2]=='2'){
					a[i].D=1;	
				}
				else{
					a[i].D=2;
				}
			}
			else{
				if(l=='1')a[i].S1=1;
				else a[i].S2=1;	
			}
		}
		
	}
	sl[0]=sr[0]=sd[0]=0;
	for(int i=1;i<=n;i++){
		sl[i]=sl[i-1]+a[i].L;
		sr[i]=sr[i-1]+a[i].R;
		sd[i]=sd[i-1]+a[i].D;
	}
	int cur=0;
	for(int i=1;i<=n;i++){
		if(a[i].S1==1){
			if(cur==0){
				cur=i;
			}
		}
		if(i==n||a[i+1].S1==0){
			if(sl[cur-1]==cur-1&&sr[n]-sr[i]==n-i){
				check('1',cur,i,sd[n]-(sd[i]-sd[cur-1]));
			}
			cur=0;
		}
	}
	cur=0;
	for(int i=1;i<=n;i++){
		if(a[i].S2==1){
			if(cur==0){
				cur=i;
			}
		}
		if(i==n||a[i+1].S2==0){
			if(sl[cur-1]==cur-1&&sr[n]-sr[i]==n-i){
				check('2',cur,i,sd[n]-(sd[i]-sd[cur-1]));
			}
			cur=0;
		}
	}
//	for(int i=1;i<=n;i++){
//		cout<<a[i].D<<' ';
//		cout<<'\n';
//	}
	int flag=-1;
	for(int i=0;i<=n;i++){
		if(sl[i]==i&&sr[n]-sr[i]==n-i){
			flag=i;
		}
	}
	if(flag!=-1){
		BN num;
		num.a.push_back(sd[n]);
		ans.push_back((sv){num,0,flag,'?'});
	}
	op(1);
	//if('?'<'0')cout<<"QAQ";
	//else cout<<"DAD";
	
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	
	int t0;
	cin>>t0;
	for(int t=0;t<t0;t++){
		solve();
	}
}