#include <bits/stdc++.h>

using i64 = long long;

#define A2 std::array<i64, 2>
#define lowbit(x) ((x) & -(x))
#define NO return (void)(std::cout << "NO" << endl)
#define YES return (void)(std::cout << "YES" << endl)
#define Fast_IOS std::ios::sync_with_stdio(false), std::cin.tie(0), std::cout.tie(0)

const i64 mod = 998244353;

template <class T> T lg(T x) {return (T)log10(x);}
template <class T> T log(T x) {return (T)log2(x);}
template <class T> T abs(T x) {return x < 0 ? -x : x;}
template <class T> T mysqrt(T x) {return std::floor(sqrtl(x));}
template <class T> T MOD(T &x, i64 p = mod) {return x = (x % p + p) % p;}
template <class T> T MOD(T &&x, i64 p = mod) {return x = (x % p + p) % p;}
template <class T1, class T2> T1 getbit(T1 x, T2 t) {return x >> t & 1ll;}

class WORK {
public:
	int N;

	WORK() {}

	void solve() {
		int n;
		std::cin >> n;
		std::vector<i64> a(n + 1);
		for (int i = 1; i <= n; i++) {
			std::cin >> a[i];
		}
		if (n == 2) {
			std::cout << 1 << '\n';
			std::cout << 1 << ' ' << a[1] << '\n';
			return;
		}
		std::vector<std::vector<i64>> bit(60);
		for (int i = 1; i <= n; i++) {
			for (int j = 0; j < 60; j++) {
				if ((a[i] >> j & 1ll) == 1ll) {
					bit[j].push_back(i);
				}
			}
		}
		std::vector<A2> ans;
		for (int i = 59; i >= 0; i--) {
			if (bit[i].empty()) continue;
			if (bit[i].size() > 2) {
				std::cout << -1 << '\n';
				return;
			}
			if (a[bit[i][0]] < a[bit[i][1]]) {
				std::swap(bit[i][0], bit[i][1]);
			}
			ans.push_back({bit[i][0], a[bit[i][1]]});
		}
		std::cout << ans.size() << '\n';
		for (auto [x, y] : ans) {
			std::cout << x << ' ' << y << '\n';
		}
	}
};

int main() {
	Fast_IOS;
	WORK work;
	int T = 1;
	std::cin >> T;
	while (T--) {
		work.solve();
	}
	return 0;
}