#include<bits/stdc++.h>
#define pb push_back
#define fi first 
#define se second 
using namespace std;
typedef pair<int,int> pii;

const int N=2e5+5;
int n1,n2,k,n,m;
int R[N];
pii p[N];
vector<int> ans;

void input()
{
    scanf("%d%d%d%d",&n1,&n2,&k,&m),n=n1+n2;
    for(int i=1,x;i<=n1;i++) scanf("%d",&x),p[i]={x,1};
    for(int i=1,x;i<=n2;i++) scanf("%d",&x),p[i+n1]={x,2};
    sort(p+1,p+n+1);
    ans.clear();
}

int calc()
{
    int nw=0,lst=0,res=0;
    for(int i=1;i<=n;i++)
    {
        if(p[i].se==2) {lst=p[i].fi; continue;}
        if(nw>=p[i].fi) continue;
        int r=R[i],x=p[i].fi;
        if(r-k<=lst||r-k<=nw) return -1;
        if(r>x+k-1) ans.pb(x),nw=x+k-1,res++;
        else ans.pb(r-k),nw=r-1,res++;
    }
    return res;
}

void solve()
{
    input();
    R[n+1]=m+1;
    for(int i=n;i>=1;i--) R[i]=p[i].se==2?p[i].fi:R[i+1];
    // for(int i=1;i<=n;i++)
    // {
    //     cerr<<' '<<p[i].fi<<' '<<p[i].se<<' '<<R[i]<<endl;

    // }
    int res=calc();
    printf("%d\n",res);
    if(res==-1) return;
    for(int x:ans) printf("%d ",x);
    printf("\n");
}

int main()
{
    int T=1;
    scanf("%d",&T);
    while(T--) solve();
}