#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int MAXN = 1e6 + 5;

int n;
ll a[MAXN];

namespace Subtask1
{
	vector< pair<ll,ll> > seq;
	bool del[10];
	bool valid;
	
	void dfs(int dep)
	{
		if(valid) return;
		if(dep == n-1)
		{
			for(auto t: seq)
				printf("%lld %lld\n",t.first,t.second);
			valid = 1;
			return;
		}
		for(int i=1; i<=n; ++i) if(!del[i])
			for(int j=1; j<i; ++j) if(!del[j] && (a[i] & 1) == (a[j] & 1))
			{
				ll x = a[i], y = a[j];
				del[j] = 1; a[i] = (x + y) >> 1;
				seq.emplace_back(x, y);
				
				dfs(dep + 1);
				
				if(valid) return;
				
				del[j] = 0; a[i] = x; a[j] = y;
				seq.pop_back();
			}
	}
	
	void solve(void)
	{
		seq.clear();
		memset(del, 0, sizeof(del));
		valid = 0;
		dfs(0);
		if(!valid)
			printf("-1\n");
	}
}

namespace Subtask2
{
	pii seq[MAXN];
	bool valid;
	
	void dfs(int len)
	{
		if(valid) return;
		if(len == 1)
		{
			for(int i=1; i<=n-1; ++i)
				printf("%d %d\n",seq[i].first,seq[i].second);
			valid = 1;
			return;
		}
		for(int i=1; i<=len; ++i)
		{
			swap(a[i], a[len]);
			for(int j=1; j<i; ++j) if((a[i] & 1) == (a[j] & 1))
			{
				swap(a[j], a[len-1]);
				
				int x = a[len], y = a[len - 1], z = (x + y) >> 1;
				a[len - 1] = z;
				seq[n-len+1] = {x, y};
				
				dfs(len - 1);
				
				a[len - 1] = y;
				swap(a[j], a[len-1]);
			}
			swap(a[i], a[len]);
		}
	}
	
	void solve(void)
	{
		valid = 0;
		
		dfs(n);
		
		if(!valid) printf("-1\n");
	}
}

void solve(void)
{
	scanf("%d",&n);
	for(int i=1; i<=n; ++i)
		scanf("%lld",&a[i]);
	
	vector<ll> vec[4];
	for(int i=1; i<=n; ++i)
		vec[a[i] % 4].emplace_back(a[i]);
	
	if(vec[1].size() || vec[3].size())
	{
		if(n <= 7)
			Subtask1 :: solve();
		else
			Subtask2 :: solve();
		return;
	}
	
	auto upd = [&] (int i,int j)
	{
		ll x = vec[i].back(); vec[i].pop_back();
		ll y = vec[j].back(); vec[j].pop_back();
		
		printf("%lld %lld\n",x,y);
		
		ll z = (x + y) / 2;
		vec[z % 4].emplace_back(z);
	};
	while(vec[0].size() >= 2 || vec[2].size() >= 2)
	{
		if(vec[0].size() >= 2) upd(0, 0);
		if(vec[2].size() >= 2) upd(2, 2);
	}
	while(vec[0].size() && vec[2].size()) upd(0, 2);
}

int main(void)
{
	int T;
	scanf("%d",&T);
	while(T--) solve();
	return 0;
}