#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=1e3;
const int mod=1e9+9;
int add(int x,int y){ return (x+=y)>=mod?x-mod:x; }
int sub(int x,int y){ return (x-=y)<0?x+mod:x; }
int ksm(ll x,int tp,int s=1){
    for(;tp;x=x*x%mod,tp>>=1) if(tp&1) s=x*s%mod;
    return s;
}
int n;
int x[N+5],y[N+5];
int a[N+5][N+5];
int ans[N+5];
bool work(){
    // for(int i=1;i<=n;i++){
    //     for(int j=0;j<=n+1;j++)
    //         printf("%d ",a[i][j]);
    //     puts("");
    // }
    int now=1;
    for(int i=0;i<=n;i++){
        int j;
        for(j=now;j<=n;j++)
            if(a[j][i]) break;
        if(j==n+1) continue;
        if(now!=j) swap(a[now],a[j]);
        int iv=ksm(a[now][i],mod-2);
        for(j=now+1;j<=n;j++){
            int tmp=1ll*a[j][i]*iv%mod;
            for(int k=i;k<=n+1;k++)
                a[j][k]=(a[j][k]-1ll*a[now][k]*tmp)%mod;
        }
        now++;
    }
    for(int i=now;i<=n;i++)
        if(a[i][n+1]) return 0;
    // for(int i=1;i<=n;i++){
    //     for(int j=0;j<=n+1;j++)
    //         printf("%d ",a[i][j]);
    //     puts("");
    // }
    for(int i=now-1,j=n;i;i--){
        int k;
        for(k=0;k<=n;k++)
            if(a[i][k]) break;
        while(j>k){
            ans[j]=1;
            for(int l=i;l;l--)
                a[l][n+1]=(a[l][n+1]-1ll*a[l][j]*ans[j])%mod;
            j--;
        }
        ans[k]=ksm(a[i][k],mod-2,a[i][n+1]);
        for(int l=i;l;l--)
            a[l][n+1]=(a[l][n+1]-1ll*a[l][j]*ans[j])%mod;
        j--;
    }
    for(int i=0;i<=n;i++)
        ans[i]=(ans[i]+mod)%mod;
    return 1;
}
void solve(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&x[i]);
    for(int i=1;i<=n;i++)
        scanf("%d",&y[i]);
    for(int i=1;i<=n;i++)
        memset(a[i],0,sizeof a[i]);
    for(int i=1;i<=n;i++){
        a[i][n+1]=y[i];
        int kp=1;
        for(int j=0;j<=n;j++){
            a[i][j]=add(a[i][j],kp);
            kp=1ll*kp*x[i]%mod;
        }
        for(int j=n+1;j<=2*n;j++){
            a[i][2*n-j]=add(a[i][2*n-j],kp);
            kp=1ll*kp*x[i]%mod;
        }
    }
    if(work()){
        printf("%d\n",2*n);
        for(int i=0;i<=n;i++)
            printf("%d ",ans[i]);
        for(int i=n-1;i>=0;i--)
            printf("%d ",ans[i]);
        puts("");
        return;
    }
    
    for(int i=1;i<=n;i++)
        memset(a[i],0,sizeof a[i]);
    for(int i=1;i<=n;i++){
        a[i][n+1]=y[i];
        int kp=1;
        for(int j=0;j<n;j++){
            a[i][j]=add(a[i][j],kp);
            kp=1ll*kp*x[i]%mod;
        }
        for(int j=n;j<2*n;j++){
            a[i][2*n-j-1]=add(a[i][2*n-j-1],kp);
            kp=1ll*kp*x[i]%mod;
        }
    }
    
    if(work()){
        printf("%d\n",2*n-1);
        for(int i=0;i<n;i++)
            printf("%d ",ans[i]);
        for(int i=n-1;i>=0;i--)
            printf("%d ",ans[i]);
        puts("");
        return;
    }
    
    puts("-1");
}
int main(){
    int t; scanf("%d",&t);
    while(t--) solve();
    return 0;
}