from math import sqrt

oa, ob, oc, od = map(int, input().split())
a = sqrt(oa)
t = oc - 2*ob - oa

def solve_quad(a, b, c):
    d = b*b - 4*a*c
    if d < 0:
        return None
    d = sqrt(d)
    return (-b+d)/(2*a), (-b-d)/(2*a)

'''
x^2 + y^2 = ob
2a*(x-y) = oc - 2ob - oa
=> x = t/(2a) + y
=> t^2/(4a^2) + 2t/(2a)*y + 2y^2 = ob
'''

'''
x^2 + y^2 = ob
-2a*(x+y) = oc - 2ob - oa
=> x = -t/(2a) - y
'''

y1, y2 = solve_quad(2, t/a, t*t/(4*a*a) - ob)
x1, x2 = t/(2*a) + y1, t/(2*a) + y2
x3, x4 = -t/(2*a) - y1, -t/(2*a) - y2

ok = False

for _x, _y in [(x1, y1), (x2, y2)]:
    '''
    A (0, a)
    B (x, y) 
    AB = (x, y-a)
    BC = (a-y, x) = AD
    C (a-y+x, x+y)
    D (a-y, a+x)
    '''
    vals = [a, _x, _y, a-_y+_x, _x+_y, a-_y, a+_x]

    if vals[1] ** 2 + vals[2] ** 2 == ob and vals[3] ** 2 + vals[4] ** 2 == oc and vals[5] ** 2 + vals[6] ** 2 == od:
        ok = True
        for val in vals:
            print("%.0f" % val, end=" ")
        break

if not ok:
    for _x, _y in [(x3, y1), (x4, y2)]:
        '''
        A (0, a)
        B (x, y)
        AB = (x, y-a)
        BC = (y-a, -x) = AD
        C (x+y-a, y-x)
        D (y-a, a-x)
        '''
        vals = [a, _x, _y, _x+_y-a, _y-_x, _y-a, a-_x]

        if vals[1] ** 2 + vals[2] ** 2 == ob and vals[3] ** 2 + vals[4] ** 2 == oc and vals[5] ** 2 + vals[6] ** 2 == od:
            for val in vals:
                print("%.0f" % val, end=" ")
            break