#include <bits/stdc++.h>
//#include <chrono>
//#pragma GCC optimize("O2")
using namespace std;
#define reps(i,s,n) for(int i = s; i < n; i++)
#define rep(i,n) reps(i,0,n)
#define Rreps(i,n,e) for(int i = n - 1; i >= e; --i)
#define Rrep(i,n) Rreps(i,n,0)
#define ALL(a) a.begin(), a.end()

using ll = long long;
using vec = vector<ll>;
using mat = vector<vec>;

ll N,M,H,W,K,A,B;
string S;
using P = pair<ll, ll>;
const ll INF = (1LL<<60);

template<class T> bool chmin(T &a, const T &b){
    if(a > b) {a = b; return true;}
    else return false;
}
template<class T> bool chmax(T &a, const T &b){
    if(a < b) {a = b; return true;}
    else return false;
}
template<class T> void print_vector(std::vector<T> v,bool endline = true){
    if(!v.empty()){
        for(std::size_t i{}; i<v.size()-1; ++i) std::cout<<v[i]<<" ";
        std::cout<<v.back();
    }
    if(endline) std::cout<<std::endl;
}

P Q(int i, int j){
    cout<<"? "<<i + 1<<' '<<j + 1<<endl;
    int p, q; cin>>p>>q;
    return make_pair(p, q);
}

vec calc(P p, int n, int o){
    vec res;
    rep(_, 2){
        int nn = p.first, oo = p.second;
        nn += n; oo += o;
        if(nn >= 1 && nn < N && oo >= 0 && oo < N) res.push_back(nn * N + oo);
        swap(p.first, p.second);
    }
    return res;
}

int main(){
    //cin.tie(nullptr);
    //ios::sync_with_stdio(false);
    cin>>N;
    map<P, vector<int> > id;
    reps(i, 1, N * N - N){
        id[Q(i, 0)].push_back(i);
    }
    P big = id.rbegin()->first, small = id.begin()->first;
    if(N == 2){
        if(small.first == -1) cout<<"! 3 2"<<endl;
        else cout<<"! 2 3"<<endl;
        return 0;
    }
    int n_bit, one_bit, small_id = id.begin()->second[0], big_id = id.rbegin()->second[0];
    if(small.first == 0){
        one_bit = 0; n_bit = 1;
    }else if(big.first < 0){
        one_bit = n_bit = N - 1;
    }else if(big.first - small.first == N) {
        one_bit = -small.first; n_bit = 1 - small.second;
    }else{
        one_bit = -small.second; n_bit = 1 - small.first;
    }
    vec ans(N * N - N);
    ans[0] = n_bit * N + one_bit;
    //cout<<n_bit<<' '<<one_bit<<endl;
    for(pair<P, vector<int> > pv : id){
        cout<<pv.first.first<<' '<<pv.first.second<<endl;
        //print_vector(pv.second);
        vec kouho = calc(pv.first, n_bit, one_bit);
        assert(!kouho.empty());
        if(pv.second.size() == 1){
            ans[pv.second[0]] = kouho[0];
        }else {
            vec nx_kouho;
            if(n_bit == one_bit) nx_kouho = calc(Q(pv.second[0], small_id), 1, 0);
            else nx_kouho = calc(Q(pv.second[0], big_id), N - 1, N - 1);
            assert(!nx_kouho.empty());
            rep(i, 2) {
                if (kouho[i] == nx_kouho[0] || (nx_kouho.size() == 2 && kouho[i] == nx_kouho[1])) ans[pv.second[0]] = kouho[i];
                else ans[pv.second[1]] = kouho[i];
            }
        }
    }
    cout<<"! ";
    print_vector(ans);
}