#include using namespace std; #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 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; } const int sign_offset=100; int main(){ cin>>N; int NN=N*N-N; Node R[2500]; R[0]={0,0}; for(int i=2; i<=NN; i++){ cout<<"? "<>R[i-1].p>>R[i-1].q; } 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); } } assert(ans1.p!=-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; assert(corner_p!=-1); 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.pmax_box.p || buf.q>max_box.q){ swap(buf.p,buf.q); ans[i]=buf+ans1; continue; } if(buf.qmax_box.p || buf.p>max_box.q){ ans[i]=buf+ans1; continue; } int rev=-1; for(int j=i+1; j>res.p>>res.q; if(corner_val.q!=0) swap(res.p,res.q); ans[i]=res+corner_val; Node revres = ans[i]-ans1; swap(revres.p,revres.q); revres=revres+ans1; ans[rev]=revres; } cout<<"!"; rep(i,NN) cout<<" "<<(ans[i].val()+N); cout<