#include<bits/stdc++.h>
#define Me(x,y) memset(x,y,sizeof x)
#define Mc(x,y) memcpy(x,y,sizeof x)
using ll=long long;
#define fi first
#define se second
using namespace std;
const int N=1e4+5,M=N*N*2+5,mod=1e9+9;const int INF=1e9+7;mt19937 rnd(time(0));
int n,m,k,x,y,z,Qh,YAI;
pair<int,int> P[N],Q[N];
ll mpow(ll x,int y=mod-2){ll Ans=1;while(y) y&1&&(Ans=Ans*x%mod),y>>=1,x=x*x%mod;return Ans;}
#define FAIL {puts("-1");return;}
ll f[N],g[N],A0,p[N];
int Find(int d);
int calc(int d){
    Qh=n;int i,j;for(i=1;i<=n;i++) Q[i]=P[i];
    for(i=1;i<=n;i++) if(P[i].fi) Q[++Qh]={mpow(P[i].fi),P[i].se*mpow(mpow(P[i].fi),d)%mod};
    sort(Q+1,Q+Qh+1);for(i=2;i<=Qh;i++) if(Q[i].fi==Q[i-1].fi&&Q[i].se^Q[i-1].se) return 0;
    if(YAI) return 0;YAI=1;
    Qh=unique(Q+1,Q+Qh+1)-Q-1;
    Me(f,0);Me(g,0);f[0]=1;for(i=1;i<=Qh;i++) for(j=Qh;~j;j--) f[j]=(f[j]*(mod-Q[i].fi)+(j?f[j-1]:0))%mod;
    for(i=1;i<=Qh;i++){
        ll Iv=mpow(mod-Q[i].fi);
        for(j=0;j<=Qh;j++) f[j]=(f[j]-(j?f[j-1]:0)+mod)%mod*Iv%mod;
        ll Ts=1;for(j=1;j<=Qh;j++) if(i^j)Ts=Ts*(Q[i].fi-Q[j].fi+mod)%mod;
        Ts=mpow(Ts)*Q[i].se%mod;
        for(j=0;j<=Qh;j++) g[j]=(g[j]+f[j]*Ts)%mod;
		for(j=Qh;~j;j--) f[j]=(f[j]*(mod-Q[i].fi)+(j?f[j-1]:0))%mod;
    }
    for(i=d+1;i<=Qh;i++) if(g[i]) return 0;
	Mc(p,g);for(i=0;i<=d;i++) p[d-i]=(p[d-i]+g[i])%mod*(mod+1)/2%mod;
    // for(i=0;i<=d;i++) cerr<<p[i]<<' ';cerr<<'\n';
    if(A0==-1&&!p[0]) {
        for(i=d+1;i<=Qh;i++) if(f[i]) return 0;
        for(i=0;i<=d;i++) p[i]=(p[i]+f[i])%mod,p[d-i]=(p[d-i]+f[i])%mod;
    }else if(~A0&&p[0]^A0){
        for(i=d+1;i<=Qh;i++) if(f[i]) return 0;
    //  for(i=0;i<=d;i++) cerr<<f[i]<<' ';cerr<<'\n';
        if((f[0]+f[d])%mod==0) return 0;
        ll s=(A0-p[0]+mod)%mod*mpow(f[0]+f[d])%mod;
        for(i=0;i<=d;i++) p[i]=(p[i]+f[i]*s)%mod,p[d-i]=(p[d-i]+f[i]*s)%mod;
    }
    if(!p[d]) return 0;
    printf("%d\n",d);
	for(i=0;i<=d;i++) printf("%lld ",p[i]);
	printf("\n");return 1;
}
void Solve(){
    int i,j;scanf("%d",&n);cerr<<n<<'\n';for(i=1;i<=n;i++) scanf("%d",&P[i].fi),cerr<<P[i].fi<<' ';cerr<<'\n';for(i=1;i<=n;i++) scanf("%d",&P[i].se),cerr<<P[i].se<<' ';cerr<<'\n';
    A0=-1;Qh=0;YAI=0;for(i=1;i<=n;i++) if(P[i].fi) Q[++Qh]=P[i];else if(A0==-1) A0=P[i].se;else if(A0!=P[i].se) FAIL;Mc(P,Q);n=Qh;
    for(i=1e4;i;i--){
        if(calc(i)) return ;
    }
    puts("-1");
}
int main(){
    int t;scanf("%d",&t);while(t--) Solve();
}