X,Z,S,Y=map(int,input().strip().split(' '))
alpha=2*Y-X-S
beta=2*Z-X-S
import math
A = int(round(math.sqrt(X)))
assert A*A==X
t_real_minus_t_imag = (2*Y-X-S) // (-2*A)
for f in [1,-1]:
    d = t_real_minus_t_imag*f
    # print(f,d)
    aa = 2
    bb = 2*d-2*A
    cc = d*d+X-Y
    delta = bb*bb-4*aa*cc
    if delta<0:
        continue
    dd = int(round(math.sqrt(delta)))
    if dd*dd!=delta:
        continue
    for s in [-dd,dd]:
        uu = -bb+s
        if uu%(2*aa)!=0:
            continue
        ty = uu//(2*aa)
        tx = ty+d
        coords = [
            (0,f*A),
            (0-ty,f*A+tx),
            (0+tx-ty,f*A+ty+tx),
            (0+tx,f*A+ty),
        ]
        if f<0:
            coords = [(x,-y) for x,y in coords]
        ll = [x*x+y*y for x,y in coords]
        if ll == [X,Z,S,Y]:
            import sys
            print(coords[0][1],end=' ')
            for u in coords[1:]:
                print(u[0],u[1],end=' ')
            sys.exit(0)