#include <bits/stdc++.h>
#define int long long
#define io ios::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define pb push_back
using namespace std;
const int N=1e6+5,mod1=1e9+7,mod2=998244353,inf=1e18;
typedef pair<int,int> pi;
int a[N],b[N];
void solve()
{
    int n,m,k,w;
    cin>>n>>m>>k>>w;
    for(int i=1;i<=n;i++)cin>>a[i];
    sort(a+1,a+1+n);
    for(int i=1;i<=m;i++)cin>>b[i];
    b[m+1]=w+1;
    m++;
    sort(b+1,b+1+m);
    vector<int>ans;
    for(int i=0;i<m;i++)
    {
        int x=b[i],y=b[i+1];
        int w=lower_bound(a+1,a+1+n,x)-a;
        if(a[w]>y)continue;
        int cnt=0;
        vector<pi>q;
        int last=x;
        int sum=0;
        for(int j=w;a[j]<=y&&j<=n;j++)
        {
            if(a[j]>last)
            {
                q.pb({a[j],a[j]+k-1});
                sum+=a[j]-last-1;
                last=a[j]+k-1;
                cnt++;
            }
        }
        int s=max(0ll,last-y+1);
        if(sum<s)
        {
            cout<<-1<<'\n';
            return ;
        }
        for(int j=cnt-1;j>=0;j--)
        {
            if(s)
            {
                ans.pb(q[j].first-s);
                if(j==0)s-=q[j].first-q[j-1].second-1,s=0;
            }
            else
            {
                ans.pb(q[j].first);
            }
        }
    }
    cout<<ans.size()<<'\n';
    sort(ans.begin(),ans.end());
    for(auto i:ans)cout<<i<<' ';
    cout<<'\n';
}
signed main()
{
    io;
    int t=1;
    cin>>t;
    while(t--)
    {
        solve();
    }
    return 0;
}