#include <bits/stdc++.h>
using i64 = long long;
void solve() {
	int n, m, k, w;
	std::cin >> n >> m >> k >> w;
	// std::map<int, int> map;
	std::vector<int> a(n), b(m + 2);
	for (int i = 0; i < n; i++) {
		// int a; std::cin >> a;
		// map[a] = 1;
		std::cin >> a[i];
	}
	for (int i = 0; i < m; i++) {
		// int b; std::cin >> b;
		// map[b] = -1;
		std::cin >> b[i];
	}
	b[m] = 0, b[m + 1] = w + 1;
	std::ranges::sort(a), std::ranges::sort(b);
	std::vector<int> end;

	for (int i = 0; i < n; i++) {
		if (end.empty() || end.back() <= a[i]) {
			end.push_back(a[i] + k);
		}
	}

	// std::cout << end.size() << '\n';
	// for (int i = 0; i < end.size(); i++) {
	// 	std::cout << end[i] - k << ' ';
	// }
	// std::cout << '\n';

	int i = end.size() - 1, j = m + 1, curr = 0;
	while (i >= 0) {
		while (j > 0 && b[j] >= end[i]) {
			j--;
			// std::cout << j << '\n';
		}
		while (i >= 0 && j > 0 && end[i] > b[j] && end[i] - k <= b[j]) {
			// curr = end[i] - b[j];
			// end[i] -= curr, j--;
			end[i] = b[j];
			j--;
		}

		i--;
		while (i >= 0 && i < end.size() - 2 && end[i] > end[i + 1] - k) {
			//end[i] -= curr;
			end[i] = end[i + 1] - k;
			if (end[i] - k <= b[j] && end[i] > b[j]) {
				std::cout << -1 << '\n';
				return;
			}
			i--;
		}
		curr = 0;
		while (i >= 0 && end[i] - k >= b[j]) {
			i--;
		}
	}
	for (int i = 0; i < end.size() - 1; i++) {
		if (end[i + 1] - end[i] < k) {
			std::cout << -1 << '\n';
			return;
		}
	}
	if (end[0] - k <= 0) {
		std::cout << -1 << '\n';
		return;
	}
	std::cout << end.size() << '\n';
	for (int i = 0; i < end.size(); i++) {
		std::cout << end[i] - k << ' ';
	}
	std::cout << '\n';
}
int main() {
	std::ios::sync_with_stdio(false);
	std::cin.tie(0), std::cout.tie(0);
	int t;
	std::cin >> t;
	while (t--) {
		solve();
	}
	return 0;
}