#define LOCAL
#include <bits/stdc++.h>
#define fir first
#define sec second
#define el '\n'

#ifdef LOCAL
#define FINISH cerr << "FINISH" << endl;
#else
#define FINISH ;
#endif

#ifdef LOCAL
#define debug(x) cerr << setw(4) << #x << " == " << x << endl
#else
#define debug(x)
#endif

#ifdef LOCAL
#define debugv(x)                   \
    cerr << setw(4) << #x << ":: "; \
    for (auto i : x)                \
        cerr << i << " ";           \
    cerr << endl
#else
#define debugv(x)
#endif

using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
ostream& operator<<(ostream& out, PII& x)
{
    out << x.fir << " " << x.sec << endl;
    return out;
}
#define int long long
typedef long double ld;
const int mod = 998244353;
const int inf = 0x3f3f3f3f;
const int N = 200020;
struct point {
	int x, y;
};
void solve()
{
    int n;
    cin >> n;
    vector<point> p(n + 1);
    for (int i = 1; i <= n; i++) {
    	cin >> p[i].x >> p[i].y;
    }
    sort(next(p.begin()), p.end(), [&](point x, point y) {
    	return x.x < y.x;
    });	

    int flag = 0;
    ld las = 0, area = 0;
    for (int i = 2; i < n; i++) {
    	if (p[i].x == p[i - 1].x) {
    		flag = 1;
    		continue;
    	}
    	
    	if (flag == 1) {
    		flag = 0;
    		int x1, y1, x2, y2;
    		if ((p[i].x + 2 * p[i - 1].x) % 3 == 0) {
	    		cout << "? " << (p[i].x + 2 * p[i - 1].x) / 3 << " " << 1 << endl;
	    	}
	    	else {
	    		cout << "? " << (p[i].x * 2 +  p[i - 1].x) << " " << 3 << endl;
	    	}
	    	cin >> x1 >> y1;
	    	if ((2 * p[i].x + p[i - 1].x) % 3 == 0) {
	    		cout << "? " << (2 * p[i].x + p[i - 1].x) / 3 << " " << 1 << endl;
	    	}
	    	else {
	    		cout << "? " << (2 * p[i].x + p[i - 1].x) << " " << 3 << endl;
	    	}
	    	cin >> x2 >> y2;

            ld x = (x1 * 1.0 / y1 + x2 * 1.0 / y2) * 0.5;
            ld nxt = x2 * 2.0 / y2 - x1 * 1.0 / y1;
            ll len = p[i].x - p[i - 1].x;
            x = x * len;

            area += x;
            las = nxt;
	    	continue;
    	}

    	int res = p[i].x;
    	if ((p[i].x + p[i - 1].x) % 2 == 0) {
    		cout << "? " << (p[i].x + p[i - 1].x) / 2 << " " << 1 << endl;
    	}
    	else {
    		cout << "? " << (p[i].x + p[i - 1].x) << " " << 2 << endl;
    	}
    	ll x, y;
    	cin >> x >> y;

        ld nxt = x * 2.0 / y - las;
        ll len = p[i].x - p[i - 1].x;
        ld xx = 1.0 * x / y * len;

        area += xx;
        las = nxt;
    }
    if (p[n].x == p[n - 1].x) {
        if (area - ll(area) <= 0.25) {
            cout << "! " << ll(area) << " 1" << endl;
        }
        else if (area - (ll)area >= 0.75){
            cout << "! " << ll(area + 0.5) << " 1" << endl;
        }
        else {
            area = area * 2;
            cout << "! " << ll(area + 0.5) << " 2" << endl;
        }
    	return;
    }
    if (flag == 1) {
    	if ((p[n].x + p[n - 1].x) % 2 == 0) {
    		cout << "? " << (p[n].x + p[n - 1].x) / 2 << " " << 1 << endl;
    	}
    	else {
    		cout << "? " << (p[n].x + p[n - 1].x) << " " << 2 << endl;
    	}
    	ll x, y;
    	cin >> x >> y;

        ll len = p[n].x - p[n - 1].x;
        ld xx = 1.0 * x / y * len;
        area += xx;


    	if (area - ll(area) <= 0.25) {
            cout << "! " << ll(area) << " 1" << endl;
        }
        else if (area - (ll)area >= 0.75){
            cout << "! " << ll(area + 0.5) << " 1" << endl;
        }
        else {
            area = area * 2;
            cout << "! " << ll(area + 0.5) << " 2" << endl;
        }
    	return;
    }

    ll len = p[n].x - p[n - 1].x;

    area += 0.5 * las * len;
    if (area - ll(area) <= 0.25) {
            cout << "! " << ll(area) << " 1" << endl;
        }
        else if (area - (ll)area >= 0.75){
            cout << "! " << ll(area + 0.5) << " 1" << endl;
        }
        else {
            area = area * 2;
            cout << "! " << ll(area + 0.5) << " 2" << endl;
        }
}
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    int T = 1;
    cin >> T;
    while (T--) {
        solve();
    }
    return 0;
}