#include<bits/stdc++.h>
#define pd push_back
#define all(A) A.begin(),A.end()
#define lb lower_bound
#define ve std::vector
typedef long long ll;
typedef long long ll;
typedef __int128 Int;
typedef unsigned long long ul;
typedef long double LD;
bool FileIfstream(std::string name){
	std::ifstream f(name.c_str());
	return f.good();
}
namespace Math{
	ll QP(ll x,ll y,ll mod){ll ans=1;for(;y;y>>=1,x=x*x%mod)if(y&1)ans=ans*x%mod;return ans;}
	ll inv(ll x,ll mod){return QP(x,mod-2,mod);}
}
const int N=2e5+10;
const int mod=998244353;
std::mt19937_64 rnd(time(NULL));
void solve(){
	//don't forget to open long long
	int n;std::cin>>n;
	ve<ll>v(n+1);
	for(int i=1;i<=n;i++)
		std::cin>>v[i];
	int T;
	if(n<=10)T=400;
	else if(n<=100)T=100;
	else if(n<=1e4)T=50;
	else if(n<=1e5)T=20;
	else T=3;
	while(T--){
		std::multiset<std::pair<ul,ll>>st[2];
		for(int i=1;i<=n;i++)
			st[v[i]&1].insert({rnd(),v[i]});
		ve<std::array<ll,2>>ans;
		bool flg=true;
		while(st[0].size()+st[1].size()>1){
			if(st[0].size()==1&&st[1].size()==1){
				flg=false;
				break;
			}
			int p;
			if(st[0].size()<=1)p=1;
			else if(st[1].size()<=1)p=0;
			else p=rnd()&1;
			auto x=*st[p].begin();st[p].erase(st[p].find(x));
			auto y=*st[p].begin();st[p].erase(st[p].find(y));
			st[((x.second+y.second)/2)&1].insert({rnd(),(x.second+y.second)/2});
			ans.pd({x.second,y.second});
		}
		if(flg){
			for(auto &p:ans)
				std::cout<<p[0]<<' '<<p[1]<<'\n';
			return;
		} 
	}
	std::cout<<-1<<'\n';
}
int main(){
	std::ios::sync_with_stdio(false);
	std::cin.tie(0);
	std::cout.tie(0);
	int T=1;
	std::cin>>T;
	while(T--)solve();

#ifndef ONLINE_JUDGE
	std::cerr<<std::fixed<<std::setprecision(10)<<1.0*clock()/CLOCKS_PER_SEC<<'\n';
#endif

	return 0;
}