// superyijin AK IOI
// wangsiyuanZP AK IOI
// #pragma GCC optimize(2)
// #pragma GCC optimize(3)
#include <bits/stdc++.h>
#define int long long
#define ll long long
#define ull unsigned long long
#define i128 __int128
#define ld long double
#define pii pair<int, int>
#define vi vector<int>
#define vpii vector<pair<int, int>>
#define pb push_back
#define m_p make_pair
#define fi first
#define se second
#define pcnt __builtin_popcountll
using namespace std;
const int inf = (sizeof(int) == 4 ? 0x3f3f3f3f : 0x3f3f3f3f3f3f3f3f);
int mod = 998244353;
const ld eps = 1e-7;
void chkmin(int &x, int y) { x = min(x, y); }
void chkmax(int &x, int y) { x = max(x, y); }
void inc(int &x, int y) { x = (x + y >= mod ? x + y - mod : x + y); }
void dec(int &x, int y) { x = (x - y < 0 ? x - y + mod : x - y); }
int ksm(int a, int b)
{
	int t = 1;
	while (b)
	{
		if (b & 1) t = t * a % mod;
		a = a * a % mod, b >>= 1;
	}
	return t;
}
int ksm(int a, int b, int p)
{
	int t = 1;
	while (b)
	{
		if (b & 1) t = t * a % p;
		a = a * a % p, b >>= 1;
	}
	return t;
}
bool Mbe;
const int N = 3e5 + 10;
struct Bit
{
	int t[N];
	Bit() { memset(t, 0, sizeof(t)); }
	int lowbit(int x) { return x & -x; }
	void add(int x, int y) { while (x < N) t[x] += y, x += lowbit(x); }
	int sum(int x)
	{
		int ans = 0;
		while (x >= 1) ans += t[x], x -= lowbit(x);
		return ans;
	}
};
int fac[N], invfac[N];
void init(int n)
{
	fac[0] = 1;
	for (int i = 1; i <= n; i++) fac[i] = fac[i - 1] * i % mod;
	invfac[n] = ksm(fac[n], mod - 2);
	for (int i = n - 1; i >= 0; i--) invfac[i] = invfac[i + 1] * (i + 1) % mod;
}
int A(int n, int m)
{
	if (m < 0 || m > n) return 0;
	return fac[n] * invfac[n - m] % mod;
}
int C(int n, int m)
{
	if (m < 0 || m > n) return 0;
	return fac[n] * invfac[m] % mod * invfac[n - m] % mod;
}
int calc(int now, vi vec)
{
	int sum = 0;
	while (now)
	{
		int x = __lg(now);
		for (auto &i : vec) i &= (1ll << x + 1) - 1;
		sort(vec.begin(), vec.end());
		int pos = lower_bound(vec.begin(), vec.end(), (1ll << x)) - vec.begin() - 1;
		if (pos < 0) return inf;
		now ^= (1ll << x) ^ vec[pos], sum += (1ll << x) - vec[pos];
		vec[pos] = 0;
	}
	return sum;
}
int n, m, ans;
int id[60][100010];
vpii v[60];
vector<vpii> V;
int check(int now, vpii vec)
{
	int sum = 0;
	while (now)
	{
		int x = __lg(now);
		vector<bool> f(min((int)v[x].size(), 65ll));
		for (auto i : vec) if (id[x][i.fi] < 65) f[id[x][i.fi]] = true;
		bool flag = false; 
		for (int i = 0; i < v[x].size() && i < 65; i++)
		{
			if (f[i]) continue;
			now ^= (1ll << x) ^ v[x][i].fi, sum += (1ll << x) - v[x][i].fi;
			vec.pb({v[x][i].se, (1ll << x) - v[x][i].fi});
			flag = true;
			break;
		}
		if (!flag) now ^= (1ll << x), sum += (1ll << x);
	}
	return sum;
}
void solve(int now, int sum, vpii vec)
{
	if (V.size() == m) return;
	if (check(now, vec) + sum > ans) return;
	if (now == 0)
	{
		V.pb(vec);
		return;
	}
	int x = __lg(now);
	vi vid;
	for (auto i : vec) vid.pb(id[x][i.fi]);
	sort(vid.begin(), vid.end());
	int pp = 0;
	for (int i = 0; i < v[x].size(); i++)
	{
		if (pp < vid.size() && vid[pp] == i)
		{
			pp++;
			continue;
		}
		int w = now ^ (1ll << x) ^ v[x][i].fi, ws = sum + (1ll << x) - v[x][i].fi;
		vec.pb({v[x][i].se, (1ll << x) - v[x][i].fi});
		solve(w, ws, vec);
		vec.pop_back();
		if (V.size() == m) return;
	}
	for (int i = 0; i < vec.size(); i++)
	{
		int w = now ^ (1ll << x), ws = now + (1ll << x);
		vec[i].se += (1ll << x);
		solve(w, ws, vec);
		vec[i].se -= (1ll << x);
		if (V.size() == m) return;
	}
}
int a[100010];
void Main(int c)
{
	V.clear();
	cin >> n >> m;
	int now = 0;
	vi vec;
	for (int i = 1; i <= n; i++)
	{
		cin >> a[i];
		now ^= a[i];
		vec.pb(a[i]); 
	}
	for (int x = 59; x >= 0; x--)
	{
		v[x].clear();
		for (int i = 1; i <= n; i++) id[x][i] = inf;
		for (int i = 1; i <= n; i++) if (!((a[i] >> x) & 1)) v[x].pb({a[i] & ((1 << x) - 1), i});
		sort(v[x].begin(), v[x].end());
		reverse(v[x].begin(), v[x].end());
		for (int i = 0; i < v[x].size(); i++) id[x][v[x][i].se] = i;
	}
	ans = calc(now, vec);
	if (ans == inf)
	{
		int p = __lg(now);
		for (int x = 60; x > p; x--)
		{
			vi v = vec;
			for (auto &i : v) i &= (1ll << x + 1) - 1;
			sort(v.begin(), v.end());
			int pos = lower_bound(v.begin(), v.end(), (1ll << x)) - v.begin() - 1;
			if (pos < 1) continue;
			int w = now ^ v[pos] ^ v[pos - 1], ws = (1ll << x + 1) - v[pos] - v[pos - 1];
			v[pos] = v[pos - 1] = 0;
			ans = min(ans, ws + calc(w, v));
		}
	}
	else solve(now, 0, {});
	cout << ans << '\n';
	cout << V.size() << '\n';
	for (auto v : V)
	{
		cout << v.size() << '\n';
		for (auto i : v) cout << i.fi << " ";
		cout << '\n';
		for (auto i : v) cout << i.se << " ";
		cout << '\n';
	}
}
bool Med;
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0), cout.tie(0);
	cerr << (&Mbe - &Med) / 1048576.0 << "MB" << endl;
	int c, T;
	cin >> c >> T;
	while (T--) Main(c);
	return 0;
}
// superyijin AK IOI
// wangsiyuanZP AK IOI