#int(input())
#map(int, input().split())
#list(map(int, input().split()))

import sys

x = 0
y = 0
z = 0

# a = list(map(int, input().split()))
# cnt = 0

def ques(x, y, z):
    # global cnt
    print(f"? {x} {y} {z}")
    sys.stdout.flush()
    D = int(input())
    # a = [24, -81, 33, 62, 20]
    # dx = x - a[0]
    # dy = y - a[1]
    # dz = z - a[2]
    # D = (dx**2 + dy**2 + dz**2) * a[3] + a[4]
    # cnt += 1
    return D

min1 = -150
max1 = 150

while max1 - min1 >= 3:
    h1 = min1 + (max1 - min1) // 3
    h2 = max1 - (max1 - min1) // 3
    d1 = ques(h1, y, z)
    d2 = ques(h2, y, z)

    if d1 <= d2:
        max1 = h2
    else:
        min1 = h1
    # print(min1, max1)

ansx = 0
dx = 1 << 32
for i in range(min1, max1+1):
    d = ques(i, y, z)
    if d < dx:
        dx = d
        ansx = i

x = ansx


min1 = -150
max1 = 150

while max1 - min1 >= 3:
    h1 = min1 + (max1 - min1) // 3
    h2 = max1 - (max1 - min1) // 3
    d1 = ques(x, h1, z)
    d2 = ques(x, h2, z)

    if d1 <= d2:
        max1 = h2
    else:
        min1 = h1

ansx = 0
dx = 1 << 32
for i in range(min1, max1+1):
    d = ques(x, i, z)
    if d < dx:
        dx = d
        ansx = i

y = ansx

min1 = -150
max1 = 150

while max1 - min1 >= 3:
    h1 = min1 + (max1 - min1) // 3
    h2 = max1 - (max1 - min1) // 3
    d1 = ques(x, y, h1)
    d2 = ques(x, y, h2)

    if d1 <= d2:
        max1 = h2
    else:
        min1 = h1

ansx = 0
dx = 1 << 32
for i in range(min1, max1+1):
    d = ques(x, y, i)
    if d < dx:
        dx = d
        ansx = i

z = ansx

print(f"! {x} {y} {z}")
# print(cnt)