#include<bits/stdc++.h>
#define int long long
using namespace std;
inline int read()
{
	int n=0,f=1,ch=getchar();
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		n=n*10+ch-'0';
		ch=getchar();
	}
	return n*f;
}
int nx[105],ny[105],tmp;
bool dfs(vector<int> v)
{
	if(v.size()==1)
	{
		for(int i=1;i<=tmp;i++)printf("%lld %lld\n",nx[i],ny[i]);
		return true;
	}
	for(int i=0;i<v.size();i++)
	{
		for(int j=i+1;j<v.size();j++)
		{
			if((v[i]+v[j])%2!=0)continue;
			vector<int>nv;
			for(int k=0;k<v.size();k++)
			{
				if(k!=i&&k!=j)nv.push_back(v[k]);
			}
			nx[++tmp]=v[i];
			ny[tmp]=v[j];
			nv.push_back((v[i]+v[j])/2);
			if(dfs(nv)==true)
			{
				return true;
			}
			tmp--;
		}
	}
	return false;
}
int a[1000005];
int sta1[1000005],top1;
int sta2[1000005],top2;
signed main()
{
	int t,n;
	t=read();
	for(int greg=1;greg<=t;greg++)
	{
		n=read();
		for(int i=1;i<=n;i++)
		{
			a[i]=read();
		}
		tmp=0;
		if(n<=7)
		{
			vector<int>v;
			for(int i=1;i<=n;i++)v.push_back(a[i]);
			if(dfs(v)==false)printf("-1\n");
			continue;
		}
		bool flag=true;
		for(int i=1;i<=n;i++)
		{
			if(a[i]%2!=0)
			{
				flag=false;
				break;
			}
		}
		if(flag==true)
		{
			top1=top2=0;
			for(int i=1;i<=n;i++)
			{
				if(a[i]%4==0)sta1[++top1]=a[i];
				else sta2[++top2]=a[i];
			}
			for(int i=1;i<=n-2;i++)
			{
				if(top1>=2)
				{
					printf("%lld %lld\n",sta1[top1],sta1[top1-1]);
					int sth=(sta1[top1]+sta1[top1-1])/2;
					top1-=2;
					sta1[++top1]=sth;
				}
				else
				{
					printf("%lld %lld\n",sta2[top2],sta2[top2-1]);
					int sth=(sta2[top2]+sta2[top2-1])/2;
					top2-=2;
					sta2[++top2]=sth;
				}
			}
			if(top1==2)printf("%lld %lld\n",sta1[1],sta1[2]);
			else if(top2==2)printf("%lld %lld\n",sta2[1],sta2[2]);
			else printf("%lld %lld\n",sta1[1],sta2[1]);
		}
	}
	return 0;
}