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=-10**9
	hiy=10**9
	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

lox=-10**9
hix=10**9
for _ in range(200):
	m1=lox+(hix-lox)/3
	m2=lox+2*(hix-lox)/3
	if evalx(m1)[0]<evalx(m2)[0]:hix=m2
	else:lox=m1
x=int(round((lox+hix)/2))
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(score)