#include<bits/stdc++.h>
using namespace std;
#define ins insert
#define era erase
#define lb lower_bound
#define ub upper_bound
#define vec vector
#define pb push_back
#define eb emplace_back
#define siz(vec) (int)((vec).size())
#define aint(vec) (vec).begin(),(vec).end()
template<class T>
void operator +=(vec<T> &a,T b){a.pb(b);}
#define pii pair<int,int>
#define x first
#define y second
#define mp make_pair
#define exc(exp) if(exp)continue;
#define stop(exp) if(exp)break;
#define ret(exp) if(exp)return;
#define quit(sth) {sth;return;}
#define deb(var) cerr<<#var<<'='<<(var)<<"; "
#define debl(var) cerr<<#var<<'='<<(var)<<";\n"
#define ll long long
#define inf (long long)(1e18)
template<class T>
bool Min(T &x,T y){return x>y?x=y,1:0;}
template<class T>
bool Max(T &x,T y){return x<y?x=y,1:0;}
const int mod=1e9+7;
void Add(int &x,int y){x=x+y<mod?x+y:x+y-mod;}
void Dec(int &x,int y){x=x>=y?x-y:x-y+mod;}
int fpm(int x,int y){
    int ans=1;for(;y;y>>=1,(x*=x)%=mod)if(y&1)(ans*=x)%=mod;return ans;
}

const int B=50;
int _gc[B][B];
vec<int> fac[B];

int k;

unordered_map<ll,int> f[B+2];

ll S;

int CNT;
ll reduce(int x,ll T){
    ll S=T>>x<<x;
    while(S){
        int i=__lg(S&(-S));
        S^=S&(-S);
        
        if(T>>i&1){
            bool fl=1;
            for(auto d:fac[i]){
                stop(d>=x);
                fl&=T>>d&1;
            }
            if(fl){
                T^=1ll<<i;
            }
        }
    }			return T;
}
bool dfs(int x,ll T,int g){
   ++CNT;
   bool TWC=0;
    if(f[x][T]){
        if(f[x][T]<k){
            k-=f[x][T];
            return 0;
        }else{
        	TWC=1;
        }
    }

    if(g==x){
		if(!TWC)++f[x][T];
        if(!(--k)){
            // debl(CNT);
            cout<<__builtin_popcountll(S)<<' ';
            for(int i=1;i<B;i++)
                if(S>>i&1)cout<<i<<' ';
            cout<<'\n';
            f[x][T]=0;
            return 1;
        }
    }
    if(x==1)return 0;
    for(int v=max(1,63^__builtin_clzll(T&((1ll<<x)-1)));v<x;v++){
        ll R=T,Q=T;
        while(Q){
            int x=__lg(Q&(-Q));
            R|=1ll<<_gc[x][v];
            Q^=Q&(-Q);
        }
//		ll R=T|(1ll<<v);
//		for(int i=0;i<B;i++)
//			if(T>>i&1)R|=1ll<<_gc[v][i];
		R=reduce(v,R);

        S^=1ll<<v;
        auto res=dfs(v,R,_gc[g][v]);
        S^=1ll<<v;
		if(!TWC)f[x][T]+=f[v][R];

        if(res){
            if(!TWC)f[x][T]=0;
            return 1;
        }
    }                               return 0;
}
void work(){
    cin>>k;
    if(!k){
        cout<<"0\n";
    }else{
        S=0;
        dfs(B,1,0);
    }
}
signed main(){
	for(int i=0;i<B;i++){
		for(int j=0;j<B;j++){
			_gc[i][j]=__gcd(i,j);
		}
	}
    for(int i=1;i<B;i++){
        for(int j=i+i;j<B;j+=i){
            fac[j]+=i;
        }
    }
    ios::sync_with_stdio(0),    
    cin.tie(0),cout.tie(0);
    int T=1;cin>>T;while(T--)work();
}