# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef vector<ll> vl;
int ctz(ll x){return __builtin_ctzll(x);}
int get(ll x){return (x>>1)&1;}
void Do(ll &x,ll y){printf("%lld %lld\n",x,y);x=(x+y)/2;}
void simple(const vl &A)
{
	vl vec[2];
	for(ll x:A) vec[get(x)].push_back(x);
	for(ll x,y;vec[0].size()>1 || vec[1].size()>1;)
	{
		int k=(vec[1].size()>1);
		x=vec[k].back();vec[k].pop_back();
		y=vec[k].back();vec[k].pop_back();
		Do(x,y);vec[get(x)].push_back(x);
	}
	if(!vec[0].empty() && !vec[1].empty()) Do(vec[0][0],vec[1][0]);
}
bool check(vl &A,vl &spc)
{
	int n=A.size(),p=0,k=0;
	for(int i=1;i<n;i++)
		if(A[0]!=A[i] && (!p || ctz(A[0]^A[i])<k)) p=i,k=ctz(A[0]^A[i]);
	if(!p || n<k+1) return false;
	ll x=A[0],y=A[p];
	A.erase(A.begin()+p);A.erase(A.begin());
	for(int i=1;i<k;i++) Do(x,y),y=A.back(),A.pop_back();
	assert(get(x)!=get(y));spc={x,y};
	return true;
}
int main()
{
	int T;
	cin>>T;
	while(T--)
	{
		int n;ll x;
		scanf("%d",&n);
		vl vec[2];
		for(int i=1;i<=n;i++) scanf("%lld",&x),vec[x&1].push_back(x);
		if(vec[0].empty()) simple(vec[1]);
		else if(vec[1].empty()) simple(vec[0]);
		else
		{
			vl spc;
			if(!check(vec[0],spc) && !check(vec[1],spc)) puts("-1");
			else
			{
				for(ll x,y;vec[0].size()>1 || vec[1].size()>1;)
				{
					int k=(vec[1].size()>1);
					x=vec[k].back();vec[k].pop_back();
					y=vec[k].back();vec[k].pop_back();
					Do(x,y);vec[x&1].push_back(x);
				}
				if(spc[0]&1) swap(vec[0],vec[1]);
				if(!vec[0].empty() && !vec[1].empty())
				{
					if(get(spc[0])!=get(vec[0][0])) swap(spc[0],spc[1]);
					if(spc[0]%8==vec[0][0]%8) Do(spc[0],vec[0][0]),Do(spc[0],spc[1]),Do(spc[0],vec[1][0]);
					else if((((spc[0]+spc[1])/2+vec[1][0])/2+vec[0][0])%2==0) Do(spc[1],spc[0]),Do(spc[1],vec[1][0]),Do(spc[1],vec[0][0]);
					else Do(spc[1],vec[0][0]),Do(spc[1],vec[1][0]),Do(spc[1],spc[0]);
				}
				else if(!vec[0].empty())
				{
					if(get(spc[0])!=get(vec[0][0])) swap(spc[0],spc[1]);
					Do(spc[0],vec[0][0]);Do(spc[0],spc[1]);
				}
				else if(!vec[1].empty()) Do(spc[0],spc[1]),Do(spc[0],vec[1][0]);
				else Do(spc[0],spc[1]);
			}
		}
	}
	return 0;
}