#if __INCLUDE_LEVEL__ == 0 #include __BASE_FILE__ namespace { void solve() { int n; cin >> n; auto query = [&](int i, int j) -> int { print('?', i, j); int x; cin >> x; return x; }; vector> pre(n * n + 1); for (int i : rep1(n)) { if (!atcoder::internal::is_prime_constexpr(i)) { continue; } for (int j : rep1(i + 1, n)) { if (i * j <= n) { continue; } if (!atcoder::internal::is_prime_constexpr(j)) { continue; } pre[i * j] = {i, j}; } } vector ans(n + n, 1); if (3 <= n) { while (true) { int i = randint(1, n); int j = randint(1, n); int x = query(i, j); if (pre[x].first) { for (int k : rep(n + n)) { if (k < n) { ans[k] = query(k + 1, j); } else { ans[k] = query(i, k - n + 1); } } for (int z : rep(2)) { auto rng = ans | views::drop(n * z) | views::take(n); int mn = ranges::min(rng); for (int& e : rng) { e /= mn; } } break; } } } if (2 <= n) { vector indices; for (int k : rep(n + n)) { if (ans[k] == 1) { indices.push_back(k); } } assert(len(indices) == 4); for (int za : rep(2)) { for (int zb : rep(2, 4)) { int x = query(indices[za] + 1, indices[zb] - n + 1); if (x == 1) { ans[indices[za ^ 1]] = query(indices[za ^ 1] + 1, indices[zb] - n + 1); ans[indices[zb ^ 1]] = query(indices[za] + 1, indices[zb ^ 1] - n + 1); } } } } print('!', ans); } } // namespace int main() { ios::sync_with_stdio(false); solve(); } #else // __INCLUDE_LEVEL__ #include using namespace std; #include template bool chmin(T& x, U&& y) { return y < x && (x = forward(y), true); } template bool chmax(T& x, U&& y) { return x < y && (x = forward(y), true); } namespace std { template istream& operator>>(istream& is, pair& p) { return is >> p.first >> p.second; } template istream& operator>>(istream& is, tuple& t) { return apply([&is](auto&... xs) -> istream& { return (is >> ... >> xs); }, t); } template istream& operator>>(istream& is, tuple&& t) { return is >> t; } template >* = nullptr> auto operator>>(istream& is, R&& r) -> decltype(is >> *begin(r)) { for (auto&& e : r) { is >> e; } return is; } template ostream& operator<<(ostream& os, const pair& p) { return os << p.first << ' ' << p.second; } template ostream& operator<<(ostream& os, const tuple& t) { auto f = [&os](const auto&... xs) -> ostream& { [[maybe_unused]] auto sep = ""; ((os << exchange(sep, " ") << xs), ...); return os; }; return apply(f, t); } template >* = nullptr> auto operator<<(ostream& os, R&& r) -> decltype(os << *begin(r)) { auto sep = ""; for (auto&& e : r) { os << exchange(sep, " ") << e; } return os; } } // namespace std template void print(Ts&&... xs) { cout << tie(xs...) << '\n'; } inline mt19937_64 mt_engine( chrono::steady_clock::now().time_since_epoch().count()); template double uniform(double a, double b) { assert(a <= b); if constexpr (Log) { assert(0 < a); return exp(uniform(log(a), log(b))); } else { return uniform_real_distribution(a, b)(mt_engine); } } template int randint(int a, int b) { assert(a <= b); if constexpr (Log) { assert(0 < a); return round(uniform(a - 0.5, b + 0.5)); } else { return uniform_int_distribution(a, b)(mt_engine); } } inline bool bernoulli(double p = 0.5) { assert(0 <= p && p <= 1); return bernoulli_distribution(p)(mt_engine); } inline auto rep(int l, int r) { return views::iota(min(l, r), r); } inline auto rep(int n) { return rep(0, n); } inline auto rep1(int l, int r) { return rep(l, r + 1); } inline auto rep1(int n) { return rep(1, n + 1); } inline auto per(int l, int r) { return rep(l, r) | views::reverse; } inline auto per(int n) { return per(0, n); } inline auto per1(int l, int r) { return per(l, r + 1); } inline auto per1(int n) { return per(1, n + 1); } inline auto len = ranges::ssize; #endif // __INCLUDE_LEVEL__