a,b,c,d=map(int,input().split())
base=(0,int(round(a**0.5)))

def add(p,q):
	return p[0]+q[0],p[1]+q[1]

def dist(q):
	return q[0]**2+q[1]**2

def eval(x,y):
	delt=(x,y)
	rot=(y,-x)
	return max(abs(dist(add(base,delt))-b),abs(dist(add(base,add(delt,rot)))-c),abs(dist(add(base,rot))-d))

def evalx(x,loy,hiy):
	for _ in range(200):
		m1=loy+(hiy-loy)/3
		m2=loy+2*(hiy-loy)/3
		if eval(x,m1)<eval(x,m2):hiy=m2
		else:loy=m1
	y=(loy+hiy)/2
	return eval(x,y),y

def search(lox,hix,loy,hiy):
	for _ in range(200):
		m1=lox+(hix-lox)/3
		m2=lox+2*(hix-lox)/3
		if evalx(m1,loy,hiy)[0]<evalx(m2,loy,hiy)[0]:hix=m2
		else:lox=m1
	x=(lox+hix)/2
	return int(round(x)),int(round(evalx(x,loy,hiy)[1]))

c1=search(-10**9,0,-10**9,0)
c2=search(-10**9,0,0,10**9)
c3=search(0,10**9,-10**9,0)
c4=search(0,10**9,0,10**9)
ops=[c1,c2,c3,c4]
# print(*[(x,y,eval(x,y)) for x,y in ops])

x,y=min(ops,key=lambda t:eval(*t))

# score,y=evalx(x)
# y=int(round(y))
A=base
B=add(base,(x,y))
C=add(base,(x+y,y-x))
D=add(base,(y,-x))
print(A[1],*B,*C,*D)
# print(eval(x,y))