/*
60 + 0 + 100 + 64 = 224.
*/
#pragma GCC optimize("Ofast")
#pragma GCC target("avx,avx2,fma")
#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 = 7e7 + 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;
//  puts("----------------------");
//  cout << ca << " " << cb << endl;
//  L(i, 1, ca) cout << G[id][i] << " "; cout << endl;
//  L(i, 1, cb) cout << G[id ^ 1][i] << " "; cout << endl;
//  puts("------------------------");
    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[id][i]); putchar(' ');
    }
    puts("");
    Write(cb); putchar(' ');
    L(i, 1, cb)
    {
        Write(G[id ^ 1][i]); putchar(' ');
    }
    puts(""); return ;
}
inline void get(int x, int y, int a, int b)
{
    if(x % a != 0 && x % b != 0) swap(x, y);
    if(x % a != 0) swap(a, b);
    Write(a); putchar(' ');
    L(i, 1, a)
    {
        Write(i); putchar(' ');
    }
    puts("");
    Write(b); putchar(' ');
    int nw = 1;
    L(i, 1, b / y)
    {
        L(j, 1, y)
        {
            Write(nw + j - 1); putchar(' ');
        }
        nw += a;
    }
    puts(""); return ;
}
inline void solve()
{
    x = read(), y = read();
    ll t = 1ll * 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, max(1ll, f), 1)
    {
        if(t % i != 0) continue;
        a = i, b = t / i;
        if(a == x || a == y) continue;
        if(x % a == 0 && y % a == 0 && x % b == 0 && y % b == 0) get(x, y, a, b);
        else work();
        return ;
    }
    puts("0"); puts("0");
    return ;
}
signed main()
{
    // freopen(".in", "r", stdin);
    T = read();
    L(i, 1, T) solve();
    return 0;
}