#include <bits/stdc++.h>
#define ll long long
#define Ma 1000005
#define mod 1e9+7
#define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#define pb push_back

using namespace std;

ll n,m,k,w;
ll cnt=0;

struct node{
	ll num,op;
	bool operator <(const node &A)const{
		return num<A.num;
	}
}t[Ma];

vector <ll> ans;

bool go(ll l,ll r)
{
	if (l>r)
		return 1;
	vector <ll> add;
	for (ll i=l;i<=r;i++)
	{
		if (!add.empty()&&add.back()+k-1>=t[i].num)
			continue;
		add.pb(t[i].num);
	}
	if (add.back()+k-1<t[r+1].num)
	{
		for (auto z:add)
			ans.pb(z);
		return 1;
	}
	vector <ll> add2;
	add2.pb(t[r+1].num-k),add.pop_back();
	for (ll i=r;i>=l;i--)
	{
		if (add2.back()>t[i].num)
			add2.pb(min(t[i].num,add.back())),add.pop_back();
	}
	if (add2.back()>t[l-1].num)
	{
		for (auto z:add2)
			ans.pb(z);
		return 1;
	}
	return 0;
}

void sol()
{
	cnt=0;
	ans.clear();
	cin>>n>>m>>k>>w;
	for (ll i=1;i<=n;i++)
	{
		ll x;
		cin>>x;
		t[++cnt]={x,0};
	}
	for (ll i=1;i<=m;i++)
	{
		ll x;
		cin>>x;
		t[++cnt]={x,1};
	}
	sort(t+1,t+cnt+1);
	t[cnt+1]={w+1,0};
	ll l=1;
	for (ll i=1;i<=cnt;i++)
	{
		if (t[i].op)
		{
			if (!go(l,i-1))
			{
				printf("-1\n");
				return;
			}
			l=i+1;
		}
	}
	go(l,cnt);
	printf("%lld\n",(ll)ans.size());
	for (auto z:ans)
		printf("%lld ",z);
	printf("\n");
}


int main()
{
	IOS
	ll tt=1;
	cin>>tt;
	while (tt--)
		sol();
	return 0;
}

/*
4
5 2 3 16
7 11 2 9 14
13 5
3 2 4 11
6 10 2
1 11
2 1 2 6
1 5
3
2 1 2 6
1 5
2

*/