#include<bits/stdc++.h>
using namespace std;

void sol() {
    int n,m,k,w; cin>>n>>m>>k>>w;
    vector<int> a(n);
    vector<int> b(m);
    for(int i=0;i<n;i++) 
        cin>>a[i];
    for(int i=0;i<m;i++) 
        cin>>b[i];

    ranges::sort(a);
    ranges::sort(b);

    b.push_back(w+1);
    vector<int> ans;

    vector<int> tmp;
    auto solve=[&](int l,int r) {
        int n=tmp.size();
        if(!n) return true;

        vector<int> pos;
        for(int i=0,j=0;i<n;i=j) {
            pos.push_back(tmp[i]);
            while(j<n&&tmp[j]<tmp[i]+k) j++;
        }

        // cout<<"range: "<<l<<" "<<r<<endl;
        // for(int i:pos) 
        //     cout<<i<<' ';
        // cout<<endl;
        auto bias=[&](auto &&self,int x,int dis) {
            pos[x]-=dis;
            if(!x) return pos[x]>l;
            if(pos[x]<pos[x-1]+k) return self(self,x-1,pos[x-1]+k-pos[x]);
            return true;
        }; 
        if(pos.back()+k>r) {
            if(!bias(bias,n-1,pos.back()+k-r)) {
                // for(int i:pos) 
                //     cout<<i<<' ';
                // cout<<endl;
                return false;
            }
        }

        for(int x:pos) ans.push_back(x);
        return true;
    };

    for(int L=0,i=0;int R:b) {
        while(i<n&&a[i]<R) {
            tmp.push_back(a[i++]);
        }
        if(!solve(L,R)) {
            // cout<<"range: "<<L<<' '<<R<<endl;
            // for(int x:tmp) cout<<x<<' ';
            // cout<<endl;
            return cout<<-1<<'\n',void();
        }
        L=R;
        tmp.clear();
    }

    cout<<ans.size()<<'\n';
    for(int x:ans) 
        cout<<x<<' ';
    cout<<'\n';
}

int main() {
    cin.tie(nullptr)->sync_with_stdio(false);

    int T; cin>>T;
    while(T--)
        sol();

    return 0;
}