#include<bits/stdc++.h>
using namespace std;
int n,m,k,w;
int a[100010];
int b[100010];
struct node {
	int f,x;
} nd[200020];
int ans[100010];
int res=0;
bool cmp(node a,node b) {
	return a.x<b.x;
}
bool solve(int l,int r,int stf,int edf) {
	int ed=-1;
	int q=res;
	for(int i=l; i<r; i++) {
		if(ed<a[i]) {
			ans[res++]=a[i];
			ed=a[i]+k-1;
		}
	}
	if(ed>=edf) {
		int mv=ed-edf+1;
		for(int i=res-1; i>=q; i--) {
			ans[i]-=mv;
			if(i==q) {
				if(ans[i]<=stf) return false;
			} else {
				mv=ans[i-1]+k-1-ans[i];
				if(mv<=0) break;
			}
		}
	}
	return true;
}
int main() {
	ios::sync_with_stdio(false);
	int T;
	cin>>T;
	while(T--) {
		cin>>n>>m>>k>>w;
		res=0;
		int ff=0;
		for(int i=0; i<n; i++) {
			cin>>a[i];
			nd[i].f=0;
			nd[i].x=a[i];
		}
		for(int i=0; i<m; i++) {
			cin>>b[i];
			nd[i+n].f=1;
			nd[i+n].x=b[i];
		}
		sort(a,a+n);
		sort(nd,nd+n+m,cmp);
		int cnt=0;
		int las=0;
		int cntf=0;
		for(int i=0; i<n+m; i++) {
			if(nd[i].f==0) {
				cnt++;
			} 
			else {
				if(cntf!=cnt) {
					if(!solve(cntf,cnt,las,nd[i].x)) {
						ff=-1;
						break;
					}
				}
				cntf=cnt;
				las=nd[i].x;
			}
		}
		if(cntf!=cnt) {
			if(!solve(cntf,cnt,las,w+1)) {
				ff=-1;
			}
		}
		if(ff==-1) cout<<-1<<endl;
		else {
			cout<<res<<endl;
			for(int i=0; i<res-1; i++) {
				cout<<ans[i]<<" ";
			}
			cout<<ans[res-1]<<endl;
		}
//		int st=0;
//		while(st<n){
//			int c=1;
//			while(st+c<n&&ed[])
//		}
	}
}