#include <bits/stdc++.h>

#define LL long long
#define ull unsigned long long
#define F(i, j, k) for (int i = j; i <= k; ++i)
#define DF(i, j, k) for (int i = j; i >= k; --i)

using namespace std;

template <typename T> inline void read(T &n) {
    T w = 1;
    n = 0;
    char ch = getchar();
    while (!isdigit(ch) && ch != EOF) {
        if (ch == '-') w = -1;
        ch = getchar();
    }
    while (isdigit(ch) && ch != EOF) {
        n = (n << 3) + (n << 1) + (ch & 15);
        ch = getchar();
    }
    n *= w;
}

template <typename T> inline void write(T x) {
    T l = 0;
    ull y = 0;
    if (!x) { putchar(48); return; }
    if (x < 0) { x = -x; putchar('-'); }
    while (x) { y = y * 10 + x % 10; x /= 10; ++l; }
    while (l) { putchar(y % 10 + 48); y /= 10; --l; }
}

template <typename T> inline void writes(T x) {
    write(x);
    putchar(' ');
}

template <typename T> inline void writeln(T x) {
    write(x);
    puts("");
}

template <typename T> inline void checkmax(T &a, T b) { a = a > b ? a : b; }

template <typename T> inline void checkmin(T &a, T b) { a = a < b ? a : b; }

const double Eps = 1e-8;

LL A, B, C, D;

inline double check(double now) {
    double xa = 0, ya = sqrt(A);
    double csita = (D - A - now * now) / (2.0 * sqrt(A) * now);
    double yd = ya + csita * now;
    double xd = sqrt(now * now - csita * now * csita * now);
    double xb1 = yd - ya, yb1 = ya - xd;
    double xb2 = ya - yd, yb2 = ya + xd;
    return min(xb1 * xb1 + yb1 * yb1, xb2 * xb2 + yb2 * yb2);
}

inline void work(double now) {
    double xa = 0, ya = sqrt(A);
    double csita = (D - A - now * now) / (2.0 * sqrt(A) * now);
    double yd = ya + csita * now;
    double xd = sqrt(now * now - csita * now * csita * now);
    double xb1 = yd - ya, yb1 = ya - xd;
    double xb2 = ya - yd, yb2 = ya + xd;
    double xb, yb;
    if (xb1 * xb1 + yb1 * yb1 < xb2 * xb2 + yb2 * yb2) {
        xb = xb1;
        yb = yb1;
    }
    else {
        xb = xb2;
        yb = yb2;
    }
    double xc = xb + xd - xa, yc = yb + yd - ya;
    if (fabs(xc * xc + yc * yc - C) <= Eps) {
        cout << round(ya) << ' ';
        cout << round(xb) << ' ' << round(yb) << ' ';
        cout << round(xc) << ' ' << round(yc) << ' ';
        cout << round(xd) << ' ' << round(yd) << '\n';
        exit(0);
    }
}

int main() {
    //freopen(".in", "r", stdin);
    //freopen(".out", "w", stdout);
    cin >> A >> B >> C >> D;
    if (B > D) swap(B, D);
    double l = fabs(sqrt(A) - sqrt(D)) + Eps, r = sqrt(A) + sqrt(D) - Eps;
    F(i, 1, 1000) {
        double mid1 = l + (r - l) / 3.0, mid2 = r - (r - l) / 3.0;
        if (check(mid1) < check(mid2)) r = mid2;
        else l = mid1;
    }
    double true_mid = l;
    l = fabs(sqrt(A) - sqrt(D)) + Eps, r = true_mid;
    F(i, 1, 1000) {
        double mid = (l + r) / 2.0;
        if (check(mid) > B) l = mid;
        else r = mid;
    }
    //work(l);
    l = true_mid, r = sqrt(A) + sqrt(D) - Eps;
    F(i, 1, 1000) {
        double mid = (l + r) / 2.0;
        if (check(mid) < B) l = mid;
        else r = mid;
    }
    work(l);
    return 0;
}