/*
60 + 0 + 100 + 64 = 224.
*/

#include <bits/stdc++.h>
using namespace std;
//#define int long long
#define ll long long
#define L(i, j, k) for (int i = (j); i <= (k); i++)
#define R(i, j, k) for (int i = (j); i >= (k); i--)
#define pb push_back
#define pii pair<int, int>
inline int read()
{
	int sum = 0, nega = 1;
	char ch = getchar();
	while (ch > '9'||ch < '0')
	{
	    if (ch == '-') nega = -1;
		ch = getchar();
	}
	while (ch <= '9' && ch >= '0') sum = sum * 10 + ch - '0', ch = getchar();
	return sum * nega;
}
const int N = 4e7 + 9, mod = 998244353;
inline void add(int &x, int y) {x = (x + y) % mod;}
inline void del(int &x, int y) {x = (x - y + mod) % mod;}
int T, x, y, siz, a, b, ca, cb;
int G[2][N], F[109], tmp[N], id = 0;
inline void Write(int X)
{
	if(X > 9) Write(X / 10);
	putchar(X % 10 + '0');
}
inline void dfs(int fx, int fy, int fa, int fb, int nw)
{
	if(nw == siz + 1) return ;
	int fac = F[nw];
	if(x % fac != 0 && y % fac != 0) {dfs(fx, fy, fa, fb, nw + 1); return ;}
	if(x % fac != 0) swap(x, y), swap(fx, fy);
	if(a % fac != 0) id ^= 1, swap(fa, fb), swap(a, b), swap(ca, cb);
	x /= fac, a /= fac;
	int val = 1, cn = 0;
	L(i, 1, fac)
	{
		L(j, 1, ca) tmp[++cn] = G[id][j] + val - 1;
		val += fx;
	}
	L(i, 1, cn) G[id][i] = tmp[i]; ca = cn;
	fx *= fac, fa *= fac;
	dfs(fx, fy, fa, fb, nw); return ;
}
inline void work()
{
	ca = cb = 1;
	G[0][1] = 1, G[1][1] = 1; id = 0;
	dfs(1, 1, 1, 1, 1);
	Write(ca); putchar(' ');
	L(i, 1, ca)
	{
		Write(G[0][i]); putchar(' ');
	}
	puts("");
	Write(cb); putchar(' ');
	L(i, 1, cb)
	{
		Write(G[1][i]); putchar(' ');
	}
	puts(""); return ;
}
inline void solve()
{
	x = read(), y = read();
	ll t = x * y, f = sqrt(t); siz = 0;
	ll ff = t;
	L(i, 2, f)
		if(ff % i == 0)
		{
			F[++siz] = i;
			while(ff % i == 0) ff /= i;
		}
	if(ff != 1)
	{
		if(ff == t) {puts("0"); puts("0"); return ;}
		F[++siz] = ff;
	}
	R(i, f, 1)
	{
		if(t % i != 0) continue;
		a = i, b = t / i;
		if(a == x || a == y) continue;
		work(); return ;
	}
	puts("0"); puts("0");
	return ;
}
signed main()
{
	T = read();
	L(i, 1, T) solve();
	return 0;
}