#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;
    });	
    ll lasx = 0, lasy = 1;
    ll areax = 0, areay = 1;
    int lft = 0;
    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;

	    	// ll nxtx = 2 * x2 * y1 - x1 * y2;
	    	// ll nxty = y1 * y2;
	    	// ll g = __gcd (nxtx, nxty);
	    	// nxtx /= g, nxty /= g;
	    	// lasx = nxtx, lasy = nxty;

	    	// ll x = x2 * y1 + x1 * y2;
	    	// ll y = y1 * y2 * 2;
	    	// g = __gcd(x, y);
	    	// x = x / g, y = y / g;

	    	// ll len = p[i].x - p[i - 1].x;
	    	// x = x * len;

	    	// ll nxtareax = x * areay + y * areax;
	    	// ll nxtareay = y * areay;
	    	// g = __gcd (nxtareax, nxtareay);
	    	// areax = nxtareax / g;
	    	// areay = nxtareay / g;

            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;

    	// ll nxtx = 2 * x * lasy - y * lasx;
    	// ll nxty = y * lasy;
    	// ll g = __gcd (nxtx, nxty);
    	// nxtx /= g, nxty /= g;
    	// lasx = nxtx, lasy = nxty;

    	// ll len = p[i].x - p[i - 1].x;
    	// x = x * len;

    	// ll nxtareax = x * areay + y * areax;
    	// ll nxtareay = y * areay;
    	// g = __gcd (nxtareax, nxtareay);
    	// areax = nxtareax / g;
    	// areay = nxtareay / g;

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

        area += x;
        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;
    	// x = x * len;

    	// ll nxtareax = x * areay + y * areax;
    	// ll nxtareay = y * areay;
    	// ll g = __gcd (nxtareax, nxtareay);
    	// areax = nxtareax / g;
    	// areay = nxtareay / g;

        ll len = p[n].x - p[n - 1].x;
        x = 1.0 * x / y * len;
        area += 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;
    }
    // ll len = p[n].x - p[n - 1].x;
	// ll nxtx = len * lasx;
    // ll nxty = 2 * lasy;
    // ll g = __gcd (nxtx, nxty);
    // nxtx /= g, nxty /= g;
    // ll nxtareax = nxty * areax + nxtx * areay;
    // ll nxtareay = nxty * areay;
    // g = __gcd (nxtareax, nxtareay);
    // areax = nxtareax / g;
    // areay = nxtareay / g;

    ll len = p[n].x - p[n - 1].x;
    area += 0.5 * las * len / 2;

    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;
}