#include <iostream>
#include <cmath>
#include <cstdlib>

using namespace std;

struct Point {
    long long x, y;

    Point(long long x, long long y) : x(x), y(y) {}
};

int main() {
    long long A, B, C, D;
    cin >> A >> B >> C >> D;

    Point pA(0, static_cast<long long>(sqrt(A)));
    long long temp;

    if (2 * pA.y == 0) {
        temp = 0;
    } else {
        temp = abs((D - B) / (2 * pA.y));
    }

    Point pC(temp, static_cast<long long>(sqrt(abs(C - temp * temp))));
    long long sideSqr = ((pA.x - pC.x) * (pA.x - pC.x) + (pA.y - pC.y) * (pA.y - pC.y)) / 2;
    long long x = (-temp + static_cast<long long>(sqrt(temp * temp - 2 * (temp * temp - sideSqr)))) / 2;

    Point p1(temp + x, pA.y + x);
    Point p2(-x, pA.y + temp + x);

    if (pC.y < pA.y) {
        p1 = Point(temp + x, pA.y - x);
        p2 = Point(-x, pA.y - temp - x);
    }

    Point pB = p2; 
    Point pD = p1;

    if (p1.x * p1.x + p1.y * p1.y == B) {
        pB = p1;
        pD = p2;
    }

    cout << pA.y << " " << pB.x << " " << pB.y << " " << pC.x << " " << pC.y << " " << pD.x << " " << pD.y << endl;

    return 0;
}