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
debug = False

def dbg(vals):
    if debug:
        print(f"B {vals[1]:.0f} {vals[2]:.0f} ({vals[1]**2+vals[2]**2 - ob == 0})")
        print(f"C {vals[3]:.0f} {vals[4]:.0f} ({vals[3]**2+vals[4]**2 - oc == 0})")
        print(f"D {vals[5]:.0f} {vals[6]:.0f} ({vals[5]**2+vals[6]**2 - od == 0})")


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)
    '''
    _x, _y = round(_x), round(_y)
    vals = [a, _x, _y, a-_y+_x, _x+_y, a-_y, a+_x]
    dbg(vals)

    if vals[1] ** 2 + vals[2] ** 2 - ob == 0 and vals[3] ** 2 + vals[4] ** 2 - oc == 0 and vals[5] ** 2 + vals[6] ** 2 - od == 0:
        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)
        '''
        _x, _y = round(_x), round(_y)
        vals = [a, _x, _y, _x+_y-a, _y-_x, _y-a, a-_x]
        dbg(vals)

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