//#define _GLIBCXX_DEBUG

#include<bits/stdc++.h>
using namespace std;

//#define endl '\n'
#define lfs cout<<fixed<<setprecision(10)
#define ALL(a)  (a).begin(),(a).end()
#define ALLR(a)  (a).rbegin(),(a).rend()
#define spa << " " <<
#define fi first
#define se second
#define MP make_pair
#define MT make_tuple
#define PB push_back
#define EB emplace_back
#define rep(i,n,m) for(ll i = (n); i < (ll)(m); i++)
#define rrep(i,n,m) for(ll i = (ll)(m) - 1; i >= (ll)(n); i--)
using ll = long long;
using ld = long double;
const ll MOD1 = 1e9+7;
const ll MOD9 = 998244353;
const ll INF = 1e18;
using P = pair<ll, ll>;
template<typename T1, typename T2>bool chmin(T1 &a,T2 b){if(a>b){a=b;return true;}else return false;}
template<typename T1, typename T2>bool chmax(T1 &a,T2 b){if(a<b){a=b;return true;}else return false;}
ll median(ll a,ll b, ll c){return a+b+c-max({a,b,c})-min({a,b,c});}
void ans1(bool x){if(x) cout<<"Yes"<<endl;else cout<<"No"<<endl;}
void ans2(bool x){if(x) cout<<"YES"<<endl;else cout<<"NO"<<endl;}
void ans3(bool x){if(x) cout<<"Yay!"<<endl;else cout<<":("<<endl;}
template<typename T1,typename T2>void ans(bool x,T1 y,T2 z){if(x)cout<<y<<endl;else cout<<z<<endl;}  
template<typename T>void debug(vector<vector<T>>&v,ll h,ll w){for(ll i=0;i<h;i++){cout<<v[i][0];for(ll j=1;j<w;j++)cout spa v[i][j];cout<<endl;}};
void debug(vector<string>&v,ll h,ll w){for(ll i=0;i<h;i++){for(ll j=0;j<w;j++)cout<<v[i][j];cout<<endl;}};
template<typename T>void debug(vector<T>&v,ll n){if(n!=0)cout<<v[0];for(ll i=1;i<n;i++)cout spa v[i];cout<<endl;};
template<typename T>vector<vector<T>>vec(ll x, ll y, T w){vector<vector<T>>v(x,vector<T>(y,w));return v;}
ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;}
vector<ll>dx={1,-1,0,0,1,1,-1,-1};vector<ll>dy={0,0,1,-1,1,-1,1,-1};
template<typename T>vector<T> make_v(size_t a,T b){return vector<T>(a,b);}
template<typename... Ts>auto make_v(size_t a,Ts... ts){return vector<decltype(make_v(ts...))>(a,make_v(ts...));}
template<typename T1, typename T2>ostream &operator<<(ostream &os, const pair<T1, T2>&p){return os << p.first << " " << p.second;}
template<typename T>ostream &operator<<(ostream &os, const vector<T> &v){for(auto &z:v)os << z << " ";cout<<"|"; return os;}
template<typename T>void rearrange(vector<ll>&ord, vector<T>&v){
  auto tmp = v;
  for(int i=0;i<tmp.size();i++)v[i] = tmp[ord[i]];
}
template<typename Head, typename... Tail>void rearrange(vector<ll>&ord,Head&& head, Tail&&... tail){
  rearrange(ord, head);
  rearrange(ord, tail...);
}
//mt19937 mt(chrono::steady_clock::now().time_since_epoch().count());
int popcount(ll x){return __builtin_popcountll(x);};
int poplow(ll x){return __builtin_ctzll(x);};
int pophigh(ll x){return 63 - __builtin_clzll(x);};
struct RandomNumberGenerator {
  mt19937_64 mt;

  RandomNumberGenerator() : mt(chrono::steady_clock::now().time_since_epoch().count()) {}

  ll operator()(ll a, ll b) { // [a, b)
    uniform_int_distribution< ll > dist(a, b - 1);
    return dist(mt);
  }

  ll operator()(ll b) { // [0, b)
    return (*this)(0, b);
  }
};
//#define _DEBUG

int main(){
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  ll res=0,buf=0;
  bool judge = true;
  vector<ll>a;
  ll n;cin>>n;
  auto ask=[&](ll i,ll j){
    ll p,q;
    #ifdef _DEBUG
    p=a[i]/n-a[j]/n;
    q=a[i]%n-a[j]%n;
    if(p>q)swap(p,q);
    cout<<"?" spa i+1 spa j+1<<endl;
    cout<<p spa q<<endl;
    #else
    cout<<"?" spa i+1 spa j+1<<endl;
    cin>>p>>q;
    #endif
    return MP(p,q);
  };
  map<ll,ll>mp;
  ll sz=n*n-n;
  #ifdef _DEBUG
  a.assign(sz,0);
  rep(i,0,sz)cin>>a[i];
  #endif
  ll k=0;
  //RandomNumberGenerator rng;
  //k=rng(0,sz);
  vector<P>p(sz);
  rep(i,0,sz){
    if(i==k)continue;
    p[i]=ask(k,i);
    mp[p[i].fi]++;
    mp[p[i].se]++;
  }
  ll xr=-1,yr=-1;
  vector<P>tmp(p.begin(),p.end());
  sort(ALL(tmp));
  rep(x,1,n)rep(y,0,n){
    vector<P>v;
    rep(i,1,n)rep(j,0,n){
      v.EB(x-i,y-j);
      if(v.back().fi>v.back().se)swap(v.back().fi,v.back().se);
    }
    sort(ALL(v));
    //debug(tmp,tmp.size());
    //debug(v,v.size());
    if(v==tmp){
      assert(xr==-1&&yr==-1);
      xr=x,yr=y;
    }
  }
  //return 0;
  assert(xr!=-1);
  //cout<<xr spa yr<<endl;
  vector<bool>dame(sz);
  rep(i,0,sz)rep(j,i+1,sz){
    if(p[i]==p[j]){
      dame[i]=true;
      dame[j]=true;
    }
  }
  rep(i,0,sz){
    if(!dame[i]){
      if(xr-p[i].fi<1||xr-p[i].fi>=n||yr-p[i].se<0||yr-p[i].se>=n){
        swap(p[i].fi,p[i].se);
      }
    }
  }
  ll s=-1,sx=-1,sy=-1;
  rep(i,0,sz){
    if(i!=k&&!dame[i]&&(p[i].fi==0||p[i].se==0)){
      s=i;
      sx=xr-p[i].fi;
      sy=yr-p[i].se;
    }
  }
  //cout<<p<<endl;
  rep(i,0,sz)rep(j,i+1,sz){
    if(p[i]==p[j]){
      assert(s!=-1);
      //cout<<i spa j spa s spa sx spa sy<<endl;
      auto tp=ask(s,i);
      auto xi=MP(xr-p[i].fi,yr-p[i].se);
      auto xj=MP(xr-p[i].se,yr-p[i].fi);
      auto si=MP(sx-xi.fi,sy-xi.se);
      if(si.fi>si.se)swap(si.fi,si.se);
      if(tp!=si)swap(p[i].fi,p[i].se);
      else swap(p[j].fi,p[j].se);
    }
  }
  vector<ll>ret(sz);
  rep(i,0,sz){
    ll x=xr-p[i].fi,y=yr-p[i].se;
    //cout<<x spa y spa x*n+y<<endl;
    ret[i]=x*n+y;
  }
  cout<<"! ";
  debug(ret,sz);
  return 0;
}