from fractions import Fraction
import sys

T = int(input())

while T > 0:
    T -= 1
    n = int(input())
    point = []
    for i in range(n):
        x, y = map(int, input().split())
        point.append([x, y])
    point.sort(key=lambda item: item[0])
    pre = Fraction(0)
    pre_x = point[0][0]
    res = Fraction(0)


    i = 1
    while i < n:
        x, y = point[i]
        if i == n - 1:
            p, q = 0, 1
            res += (Fraction(p, q) + pre) * (x - pre_x)
            break
        j = i
        while j < n and point[j][0] == x:
            j += 1
        cnt = j - i
        # print(f'index = {i}, {j}')
        if cnt == 1:
            print(f"? {x} {1}")
            sys.stdout.flush()
            p, q = map(int, input().split())
            res += (Fraction(p, q) + pre) * (x - pre_x)
            pre = Fraction(p, q)
            pre_x = x
        else:
            cur = Fraction(x + pre_x, 2)
            print(f"? {cur.numerator} {cur.denominator}")
            sys.stdout.flush()
            p, q = map(int, input().split())
            res += (Fraction(p, q) + pre) * (Fraction(x + pre_x, 2) - pre_x)
            res += (3 * Fraction(p, q) - pre) * (x - Fraction(x + pre_x, 2))
            print(f"? {x} {1}")
            sys.stdout.flush()
            p, q = map(int, input().split())
            pre = Fraction(p, q)
            pre_x = x

        i = j
    res /= 2
    print(f"! {res.numerator} {res.denominator}")
    sys.stdout.flush()