// Slower C++ solution:
#include <bits/stdc++.h>

using namespace std;



using i32 = std::int_fast32_t;

using i64 = std::int_fast64_t;

template <typename T> constexpr T inf = numeric_limits<T>::has_infinity ? numeric_limits<T>::infinity() : numeric_limits<T>::max() / 4;



#define REP(i, stop) FOR(i, 0, stop)

#define FOR(i, start, stop) for (int i = start, i##_len = stop; i < i##_len; ++i)

#define RREP(i, n) for (int i = n; i-- > 0;)

#define ALL(v) (v).begin(), (v).end()



template<typename T, typename U> istream& operator>>(istream& is, pair<T, U>& a) { return is >> a.first >> a.second; }

template<typename T> istream& operator>>(istream& is, vector<T>& a) { for (T& x : a) is >> x; return is; }



void LOG(const int& x) { cerr << x << "\n"; }

void LOG(const vector<int>& a) { for (int x : a) cerr << x << " "; cerr << "\n"; }

template<typename T, std::size_t S> void LOG(const T (&a)[S]) { REP(i, S) cerr << a[i] << " "; cerr << "\n"; }

template<typename T, std::size_t S, std::size_t R> void LOG(const T (&a)[S][R]) { REP(i, S) { REP(j, R) cerr << a[i][j] << " "; cerr << "\n"; } }



int ceil(int a, int b) { return (a - 1) / b + 1; }

i64 gcd(i64 a, i64 b) { while (b != 0) { i64 t = b; b = a % b; a = t; } return a; }

i64 lcm(i64 a, i64 b) { return a / gcd(a, b) * b; }

i64 gcd(const vector<i64>& v) { return accumulate(ALL(v), 1LL, static_cast<i64(*)(i64, i64)>(gcd)); }

i64 lcm(const vector<i64>& v) { return v.empty() ? 0LL : accumulate(v.begin() + 1, v.end(), v[0], static_cast<i64(*)(i64, i64)>(lcm)); }

i64 pow(i64 x, i64 y, i64 z) { i64 a = 1; while (y > 0) { if (y & 1) a = a * x % z; x = x * x % z; y /= 2; } return a; }

i64 inv(i64 x, i64 m) { return pow(x, m - 2, m); }



#define COMP(key) [](const auto& a, const auto& b) { return a.key < b.key; }



struct InitCpp { InitCpp() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(10); } } initCpp;



signed main() {

  i32 x, y, z, k;

  cin >> x >> y >> z >> k;

  vector<i64> A(x), B(y), C(z);

  cin >> A >> B >> C;

  vector<i64> AB(x * y);

  REP(i, x) {

    REP(j, y) {

      AB[i * y + j] = A[i] + B[j];

    }

  }

  sort(ALL(AB));

  i32 w = min(k, x * y);

  vector<i64> ABC(w * z);

  REP(i, w) {

    REP(j, z) {

      ABC[i * z + j] = AB[AB.size() - i - 1] + C[j];

    }

  }

  sort(ALL(ABC));

  auto iter = ABC.end();

  REP(_, k) {

    iter--;

    cout << *iter << endl;

  }

}

// Faster C++ solution for the same problem (measured speedup: 1.104306259499911x):
#include <cstdio>

#include <algorithm>

#include <set>

using namespace std;



using i32 = int;

using i64 = long long;



#define REP(i, stop) for (int i = 0; i < stop; ++i)



template<typename T>

T read() {

  T a = 0;

  register int c;

  for(;;) {

    c = getchar_unlocked();

    if (c & 16) {

      a = a * 10 + c - 48;

    } else {

      return a;

    }

  }

}



void writeln(i64 a) {

  char buf[20];

  register int i = 0;

  // if (a <= 0)

  while (a > 0) {

    buf[i] = a % 10;

    a = a / 10;

    ++i;

  }

  while (i > 0) {

    --i;

    putchar_unlocked(buf[i] + 48);

  }

  putchar_unlocked('\n');

}



i64 A[1000], B[1000], C[1000];

i64 AB[3000];

signed main() {

  i32 X = read<i32>(),

      Y = read<i32>(),

      Z = read<i32>(),

      K = read<i32>();

  REP(i, X) A[i] = -read<i64>();

  REP(i, Y) B[i] = -read<i64>();

  REP(i, Z) C[i] = -read<i64>();

  sort(A, A + X);

  sort(B, B + Y);

  sort(C, C + Z);

  using P = pair<i64, i32>;

  set<P> Q, R;

  Q.insert({A[0] + B[0], 0});

  i32 L = min(K, X * Y);

  REP(i, L) {

    i64 p = (*Q.begin()).first;

    i64 q = (*Q.begin()).second;

    AB[i] = p;

    Q.erase(Q.begin());

    i32 u = q >> 10;

    i32 v = q & 1023;

    if (u + 1 < X && (u + 1) * v < L) Q.insert({A[u + 1] + B[v], q + 1024});

    if (v + 1 < Y && u * (v + 1) < L) Q.insert({A[u] + B[v + 1], q + 1});

  }

  R.insert({AB[0] + C[0], 0});

  REP(_, K) {

    i64 p = (*R.begin()).first;

    i64 q = (*R.begin()).second;

    // printf("%lld\n", -p);

    writeln(-p);

    R.erase(R.begin());

    i32 u = q >> 10;

    i32 v = q & 1023;

    if (u + 1 < L && (u + 1) * v < K) R.insert({AB[u + 1] + C[v], q + 1024});

    if (v + 1 < Z && u * (v + 1) < K) R.insert({AB[u] + C[v + 1], q + 1});

  }

}