#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();
        for(int &x:tmp) x-=l;

        vector<int> dp(n),from(n,-1);
        dp[0]=0;

        for(int i=1,j=0;i<n;i++) {
            while(tmp[i]-tmp[j+1]+1>k) j++;
            dp[i]=max(dp[j]+k,tmp[i]);
            from[i]=j;
        }

        if(dp.back()>=r) return false;
        int x=n-1;
        vector<int> path;
        while(x) {
            path.push_back(dp[x]+l-k+1);
            x=from[x];
        }
        while(!path.empty()) {
            ans.push_back(path.back());
            path.pop_back();
        }
        return true;
    };

    for(int L=0,i=0;int R:b) {
        tmp.push_back(L);
        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;
}