#include <bits/stdc++.h>

#define _overload(_1,_2,_3,name,...) name
#define _rep(i,n) _range(i,0,n)
#define _range(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload(__VA_ARGS__,_range,_rep,)(__VA_ARGS__)

#define _rrep(i,n) _rrange(i,n,0)
#define _rrange(i,a,b) for(int i=int(a)-1;i>=int(b);--i)
#define rrep(...) _overload(__VA_ARGS__,_rrange,_rrep,)(__VA_ARGS__)

#define _all(arg) begin(arg),end(arg)
#define uniq(arg) sort(_all(arg)),(arg).erase(unique(_all(arg)),end(arg))
#define getidx(ary,key) lower_bound(_all(ary),key)-begin(ary)
#define clr(a,b) memset((a),(b),sizeof(a))
#define bit(n) (1LL<<(n))
#define popcount(n) (__builtin_popcountll(n))

using namespace std;

template<class T>bool chmax(T &a, const T &b) { return (a<b)?(a=b,1):0;}
template<class T>bool chmin(T &a, const T &b) { return (b<a)?(a=b,1):0;}

using ll=long long;
using R=long double;
const R EPS=1e-9L; // [-1000,1000]->EPS=1e-8 [-10000,10000]->EPS=1e-7
inline int sgn(const R& r){return(r > EPS)-(r < -EPS);}
inline R sq(R x){return sqrt(max(x,0.0L));}

const int dx[8]={1,0,-1,0,1,-1,-1,1};
const int dy[8]={0,1,0,-1,1,1,-1,-1};

// Problem Specific Parameter:
using ll=long long;

ll a[3][3];
ll x[50010],y[50010];

int main(void){
	int m;
	cin >> m;
	rep(i,m) cin >> x[i] >> y[i];

	cout << "? 0 0" << endl << flush;
	cin >> a[0][2] >> a[1][2];

	cout << "? 0 1" << endl << flush;
	cin >> a[0][1] >> a[1][1];
	
	a[0][1]-=a[0][2];
	a[1][1]-=a[1][2];

	cout << "? 1 1" << endl << flush;
	cin >> a[0][0] >> a[1][0];

	a[0][0]-=a[0][1]+a[0][2];
	a[1][0]-=a[1][1]+a[1][2];

	cout << "!" << endl;
	rep(i,m){
		ll ax=a[0][0]*x[i]+a[0][1]*y[i]+a[0][2];
		ll ay=a[1][0]*x[i]+a[1][1]*y[i]+a[1][2];
		cout << ax << " " << ay << endl;
	}	

	return 0;
}