//Author: Kevin
#include<bits/stdc++.h>
#pragma GCC optimize(3)
using namespace std;
#define ll long long
#define ull unsigned ll
#define pb emplace_back
#define mp make_pair
#define ALL(x) (x).begin(),(x).end()
#define rALL(x) (x).rbegin(),(x).rend()
#define srt(x) sort(ALL(x))
#define rev(x) reverse(ALL(x))
#define rsrt(x) sort(rALL(x))
#define sz(x) (int)(x.size())
#define inf 0x3f3f3f3f
#define pii pair<int,int>
#define lb(v,x) (int)(lower_bound(ALL(v),x)-v.begin())
#define ub(v,x) (int)(upper_bound(ALL(v),x)-v.begin())
#define uni(v) v.resize(unique(ALL(v))-v.begin())
#define longer __int128_t
void die(string S){puts(S.c_str());exit(0);}
mt19937_64 rnd(20210448);
vector<pair<ll,ll>> op;
bool dfs(vector<ll> vec)
{
	if(sz(vec)==1) return true;
	for(int i=0;i<sz(vec);i++)
		for(int j=i+1;j<sz(vec);j++) if((vec[i]^vec[j]^1)&1)
		{
			op.pb(vec[i],vec[j]);
			vector<ll> vec2;
			for(int k=0;k<sz(vec);k++)
				if(i!=k&&j!=k)
					vec2.pb(vec[k]);
			vec2.pb(vec[i]+vec[j]>>1);
			if(dfs(vec2)) return true;
			op.pop_back();
		}
	return false;
}
void brute_force(vector<ll> vec)
{
	dfs(vec);
}
int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int t;
	cin>>t;
	while(t--)
	{
		op.clear();
		int n;
		cin>>n;
		vector<ll> vec;
		unordered_multiset<ll> ms[4];
		bool Flag=false;
		if(n>7) Flag=true;
		bool Flag2=false;
		for(int i=0;i<n;i++)
		{
			ll x;
			cin>>x;
			if(x>10) Flag2=true;
			vec.pb(x);
			ms[x&3].insert(x);
		}
		Flag&=Flag2;
		const int b=Flag?1:7;
		bool Killed=false;
		while(n-sz(op)>b)
		{
			if(sz(ms[0])&&sz(ms[2])) break;
			if(sz(ms[1])&&sz(ms[3])) break;
			int d=0;
			while(d<60)
			{
				int msk=0;
				for(auto val:ms[0])
					msk|=1<<(val>>d&1);
				for(auto val:ms[2])
					msk|=1<<(val>>d&1);
				if(msk==3)
				{
					ll A,B;
					for(auto val:ms[0])
						if(val>>d&1)
							A=val;
						else
							B=val;
					for(auto val:ms[2])
						if(val>>d&1)
							A=val;
						else
							B=val;
					if(sz(ms[0]))
					{
						ms[0].erase(ms[0].find(A));
						ms[0].erase(ms[0].find(B));
					}
					else
					{
						ms[2].erase(ms[2].find(A));
						ms[2].erase(ms[2].find(B));
					}
					ll C=A+B>>1;
					ms[C&3].insert(C);
					op.pb(A,B);
					break;
				}
				d++;
			}
			if(d==60)
			{
				d=0;
				while(d<60)
				{
					int msk=0;
					for(auto val:ms[1])
						msk|=1<<(val>>d&1);
					for(auto val:ms[3])
						msk|=1<<(val>>d&1);
					if(msk==3)
					{
						ll A,B;
						for(auto val:ms[1])
							if(val>>d&1)
								A=val;
							else
								B=val;
						for(auto val:ms[3])
							if(val>>d&1)
								A=val;
							else
								B=val;
						if(sz(ms[1]))
						{
							ms[1].erase(ms[1].find(A));
							ms[1].erase(ms[1].find(B));
						}
						else
						{
							ms[3].erase(ms[3].find(A));
							ms[3].erase(ms[3].find(B));
						}
						ll C=A+B>>1;
						ms[C&3].insert(C);
						op.pb(A,B);
						break;
					}
					d++;
				}
				if(d==60)
				{
					Killed=true;
					break;
				}
			}
		}
		if(sz(ms[0])&&sz(ms[2]))
		{
			vector<ll> vec;
			vec.pb(*ms[0].begin());
			vec.pb(*ms[2].begin());
			ms[0].erase(ms[0].begin());
			ms[2].erase(ms[2].begin());
			multiset<ll> ms2[2];
			for(int i=0;i<4;i++)
				for(auto x:ms[i])
					ms2[i&1].insert(x);
			while(sz(ms2[0])>1||sz(ms2[1])>1)
			{
				if(sz(ms2[0])>1)
				{
					ll A=*ms2[0].begin();
					ll B=*ms2[0].rbegin();
					ms2[0].erase(ms2[0].lower_bound(A));
					ms2[0].erase(ms2[0].lower_bound(B));
					ll C=A+B>>1;
					op.pb(A,B);
					ms2[C&1].insert(C);
				}
				else
				{
					ll A=*ms2[1].begin();
					ll B=*ms2[1].rbegin();
					ms2[1].erase(ms2[1].lower_bound(A));
					ms2[1].erase(ms2[1].lower_bound(B));
					ll C=A+B>>1;
					op.pb(A,B);
					ms2[C&1].insert(C);
				}
			}
			for(int i=0;i<2;i++)
				for(auto x:ms2[i])
					vec.pb(x);
			brute_force(vec);
		}
		else if(sz(ms[1])&&sz(ms[3]))
		{
			vector<ll> vec;
			vec.pb(*ms[1].begin());
			vec.pb(*ms[3].begin());
			ms[1].erase(ms[1].begin());
			ms[3].erase(ms[3].begin());
			multiset<ll> ms2[2];
			for(int i=0;i<4;i++)
				for(auto x:ms[i])
					ms2[i&1].insert(x);
			while(sz(ms2[0])>1||sz(ms2[1])>1)
			{
				if(sz(ms2[0])>1)
				{
					ll A=*ms2[0].begin();
					ll B=*ms2[0].rbegin();
					ms2[0].erase(ms2[0].lower_bound(A));
					ms2[0].erase(ms2[0].lower_bound(B));
					ll C=A+B>>1;
					op.pb(A,B);
					ms2[C&1].insert(C);
				}
				else
				{
					ll A=*ms2[1].begin();
					ll B=*ms2[1].rbegin();
					ms2[1].erase(ms2[1].lower_bound(A));
					ms2[1].erase(ms2[1].lower_bound(B));
					ll C=A+B>>1;
					op.pb(A,B);
					ms2[C&1].insert(C);
				}
			}
			for(int i=0;i<2;i++)
				for(auto x:ms2[i])
					vec.pb(x);
			brute_force(vec);
		}
		else if(n-sz(op)<=b&&!Killed)
		{
			vector<ll> vec;
			for(int i=0;i<4;i++)
				for(auto x:ms[i])
					vec.pb(x);
			brute_force(vec);
		}
		if(sz(op)!=n-1)
			cout<<-1<<'\n';
		else for(auto pr:op)
			cout<<pr.first<<" "<<pr.second<<'\n';
	}
	return 0;
}