#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+4;
int n;vector<ll>a[4];
vector<pair<ll,ll> >ans;
void upd(int s,int t){
	ll x=a[s].back();a[s].pop_back();
	ll y=a[t].back();a[t].pop_back();
	ans.emplace_back(x,y);
	x=x+y>>1,a[x&3].emplace_back(x);
}
void sub(int s){
	int t=s+2;
	while(a[s].size()>1||a[t].size()>1){
		while(a[s].size()>1)upd(s,s);
		while(a[t].size()>1)upd(t,t);
	}
	if(!a[s].empty()&&!a[t].empty())upd(s,t);
	for(auto&[x,y]:ans)cout<<x<<' '<<y<<'\n';
}
bool check(int s){
	if(!a[s].empty()&&!a[s+2].empty())return 1;
	if(a[s].empty())s+=2;
	for(int i=1;i<a[s].size();i++){
		ll x=a[s][i]^a[s][i-1];
		if(x&&__builtin_ctz(x)<a[s].size())return 1;
	}
	return 0;
}
void calc(int s){
	if(!a[s].empty()&&!a[s+2].empty())return;
	if(a[s].empty())s+=2;
	int ret=100,pos;
	for(int i=1;i<a[s].size();i++){
		ll x=a[s][i]^a[s][i-1];
		if(x&&__builtin_ctz(x)<ret)ret=__builtin_ctz(x),pos=i;
	}
	ans.emplace_back(a[s][pos-1],a[s][pos]);
	ll A=a[s][pos-1]+a[s][pos]>>1;
	a[s].erase(a[s].begin()+pos-1),a[s].erase(a[s].begin()+pos-1),a[A&3].emplace_back(A);
	for(int i=ret-1;i>1;i--)upd(s,s);
}
void sol(int s){
	int t=s+2;
	ll x=a[s].back(),y=a[t].back();a[s].pop_back(),a[t].pop_back();
	while(a[0].size()+a[2].size()>1||a[1].size()+a[3].size()>1){
		while(a[0].size()>1)upd(0,0);
		while(a[1].size()>1)upd(1,1);
		while(a[2].size()>1)upd(2,2);
		while(a[3].size()>1)upd(3,3);
		while(!a[0].empty()&&!a[2].empty())upd(0,2);
		while(!a[1].empty()&&!a[3].empty())upd(1,3);
	}
	if(a[0].size()+a[2].size()>0&&a[1].size()+a[3].size()>0){
		ll S=s,T=s^1;if(a[S].empty())S^=2;if(a[T].empty())T^=2;
		S=a[S][0],T=a[T][0];vector<ll>vec;
		if(~((x>>2)^(S>>2))&1)ans.emplace_back(x,S),vec={x+S>>1,y,T};
		else if(~((y>>2)^(S>>2))&1)ans.emplace_back(y,S),vec={y+S>>1,x,T};
		else ans.emplace_back(x,y),vec={x+y>>1,S,T};
		if((vec[0]&1)==(vec[1]&1))ans.emplace_back(vec[0],vec[1]),ans.emplace_back(vec[0]+vec[1]>>1,vec[2]);
		else if((vec[0]&1)==(vec[2]&1))ans.emplace_back(vec[0],vec[2]),ans.emplace_back(vec[0]+vec[2]>>1,vec[1]);
		else ans.emplace_back(vec[1],vec[2]),ans.emplace_back(vec[1]+vec[2]>>1,vec[0]);
	}else{
		if(a[s].size()+a[t].size()>0){
			if(!a[s].empty()){
				ans.emplace_back(x,a[s][0]),x=x+a[s][0]>>1;
				ans.emplace_back(x,y);
			}else{
				ans.emplace_back(y,a[t][0]),y=y+a[t][0]>>1;
				ans.emplace_back(x,y);
			}
		}else if(a[s^1].size()+a[t^1].size()>0){
			ans.emplace_back(x,y),x=x+y>>1;
			if(!a[s^1].empty())ans.emplace_back(x,a[s^1][0]);
			if(!a[t^1].empty())ans.emplace_back(x,a[t^1][0]);
		}else ans.emplace_back(x,y);
	}
	for(auto&[x,y]:ans)cout<<x<<' '<<y<<'\n';
}
void sol(){
	cin>>n,a[0].clear(),a[1].clear(),a[2].clear(),a[3].clear(),ans.clear();
	for(int i=0;i<n;i++){ll x;cin>>x;a[x&3].emplace_back(x);}
	if(a[1].empty()&&a[3].empty()){sub(0);return;}
	if(a[0].empty()&&a[2].empty()){sub(1);return;}
	if(check(0)){calc(0),sol(0);return;}
	if(check(1)){calc(1),sol(1);return;}
	cout<<-1<<'\n';
}
int main(){
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	int tc;cin>>tc;while(tc--)sol();
}