#include<bits/stdc++.h>
using namespace std;
using LL=long long;
using ULL=unsigned long long;
#define rep(i,n) for(int i=0;i<(n);i++)

int N;

struct Node{
  int p,q;
  int val()const{return p*N+q;}
  Node operator+(Node r)const{return {p+r.p,q+r.q};}
  Node operator-(Node r)const{return {p-r.p,q-r.q};}
  bool operator==(Node r)const{return p==r.p && q==r.q;}
  Node swap()const{return {q,p};}
};

Node f(Node l,Node r){ Node x={l.p-r.p,l.q-r.q}; if(x.p>x.q) swap(x.p,x.q); return x; }
Node query(int i,int j){ cout<<"? "<<i<<" "<<j<<endl; Node q; cin>>q.p>>q.q; return q;}

const int sign_offset=100;

int main(){
  cin>>N;
  int NN=N*N-N;
  Node R[2500];
  rep(i,NN) R[i] = query(i+1,1);
  int C[sign_offset*2]={}; rep(i,NN){ C[R[i].p+sign_offset]++; C[R[i].q+sign_offset]++; }

  Node ans1={-1,-1};

  rep(i,sign_offset*2-1){
    if(C[i]+N <= C[i+1]){ ans1.p = -(i-sign_offset)-1; for(int j=i+1; j<=i+N-1; j++) C[j]-=N; }
    if(C[i]+N-1 <= C[i+1]){ ans1.q = -(i-sign_offset)-1; for(int j=i+1; j<=i+N; j++) C[j]-=(N-1); }
  }

  Node corner_val= {0,N-1};
  if(ans1.q > N-1-ans1.q) corner_val={N-2,0};
  int corner_p = -1;
  rep(i,NN) if(f(corner_val,ans1)==R[i]) corner_p=i;

  Node ans[2500]; rep(i,NN) ans[i]={-1,-1};
  Node min_box = Node{0,0}-ans1;
  Node max_box = Node{N-2,N-1}-ans1;

    rep(i,NN){
      if(ans[i].p!=-1) continue;
      Node buf = R[i];
      if(buf.p==buf.q){ ans[i]=ans1+buf; continue; }
      if(buf.p<min_box.p || buf.q<min_box.q || buf.p>max_box.p || buf.q>max_box.q){
        ans[i]=buf.swap()+ans1;
        continue;
      }
      if(buf.q<min_box.p || buf.p<min_box.q || buf.q>max_box.p || buf.p>max_box.q){
        ans[i]=buf+ans1;
        continue;
      }

      int rev=-1; for(int j=i+1; j<NN; j++) if(R[i]==R[j]) rev=j;
      Node res=query(i+1,corner_p+1); if(corner_val.q!=0) res=res.swap();
      ans[i]=res+corner_val;
      Node revres = (ans[i]-ans1).swap()+ans1;
      ans[rev]=revres;
    }

  cout<<"!";
  rep(i,NN) cout<<" "<<(ans[i].val()+N);
  cout<<endl;

  return 0;
}