#include<bits/stdc++.h>
#define ll long long
#define eb emplace_back
#define ep emplace
#define pii pair<int,int>
#define fi first
#define se second
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define mems(arr,x) memset(arr,x,sizeof(arr))
#define memc(arr1,arr2) memcpy(arr1,arr2,sizeof(arr2))
using namespace std;
const int maxn=1e5+10,V=2e9;
int n,m,k,w,ans;
int a[maxn],b[maxn];
vector<int> v[maxn],out;
bool solve(int x){
    // printf("v[%d] : ",x);for(int i:v[x]) printf("%d ",i);puts("");
    bool flag=true;
    vector<pii> vec;vec.clear();
    for(int i=0;i<v[x].size();i++){
        int l=v[x][i];
        // printf("i = %d [%d,%d]\n",i,l,l+k);
        ans++;
        vec.eb(l,l+k-1);
        if(l+k-1>=b[x]){flag=false;break;}
        else    i=lower_bound(v[x].begin(),v[x].end(),l+k)-v[x].begin()-1;
    }
    if(flag){
        for(auto [l,r]:vec) out.eb(l);
        return true;
    }
    // puts("fuck");
    int L=b[x]-1;bool ret=false;
    for(int i=vec.size()-1;~i;i--){
        int l=vec[i].fi,r=vec[i].se;
        if(r<=L){ret=true;break;}
        r=L;
        l=r-k+1;
        vec[i]=pii(l,r);
        L=l-1;
    }
    if(L>=b[x-1]) ret=true;
    for(auto [l,r]:vec) out.eb(l);
    return ret;
}
void matt(){
    out.clear();
    scanf("%d%d%d%d",&n,&m,&k,&w);
    for(int i=1;i<=n;i++)   scanf("%d",&a[i]);
    for(int i=1;i<=m;i++)   scanf("%d",&b[i]);
    sort(a+1,a+n+1);sort(b+1,b+m+1);
    b[0]=0;
    b[m+1]=w+1;
    for(int i=1;i<=n;i++){
        int k=lower_bound(b+1,b+m+2,a[i])-b;
        v[k].eb(a[i]);
    }
    for(int i=1;i<=m+1;i++)
        if(!solve(i))   return puts("-1"),void();
    printf("%d\n",out.size());
    for(int i:out)  printf("%d ",i);puts("");
}
int main(){
    int T;scanf("%d",&T);while(T--)matt();
}