結果

問題 No.2496 LCM between Permutations
ユーザー maeshunmaeshun
提出日時 2023-10-06 22:34:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 3,438 bytes
コンパイル時間 4,090 ms
コンパイル使用メモリ 234,484 KB
実行使用メモリ 25,836 KB
平均クエリ数 928.07
最終ジャッジ日時 2024-07-26 16:40:22
合計ジャッジ時間 7,884 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 5 WA * 21 RE * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
#define rep(i, n) for(int i=0;i<(n);++i)
#define rep1(i, n) for(int i=1;i<=(n);i++)
#define ll long long
using mint = modint998244353;
using P = pair<ll,ll>;
using lb = long double;
using T = tuple<ll, ll, ll>;
#ifdef LOCAL
#  include <debug_print.hpp>
#  define dbg(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__)
#else
#  define dbg(...) (static_cast<void>(0))
#endif

// Sieve of Eratosthenes
// https://youtu.be/UTVg7wzMWQc?t=2774
//long longだとおそい
struct Sieve {
  int n;
  vector<int> f, primes;
  Sieve(int n=1):n(n), f(n+1) {
    f[0] = f[1] = -1;
    for (ll i = 2; i <= n; ++i) {
      if (f[i]) continue;
      primes.push_back(i);
      f[i] = i;
      for (ll j = i*i; j <= n; j += i) {
        if (!f[j]) f[j] = i;
      }
    }
  }
  bool isPrime(int x) { return f[x] == x;}
  vector<int> factorList(int x) {
    vector<int> res;
    while (x != 1) {
      res.push_back(f[x]);
      x /= f[x];
    }
    return res;
  }
  vector<P> factor(int x) {
    vector<int> fl = factorList(x);
    if (fl.size() == 0) return {};
    vector<P> res(1, P(fl[0], 0));
    for (int p : fl) {
      if (res.back().first == p) {
        res.back().second++;
      } else {
        res.emplace_back(p, 1);
      }
    }
    return res;
  }
  vector<pair<ll,int>> factor(ll x) {
    vector<pair<ll,int>> res;
    for (int p : primes) {
      int y = 0;
      while (x%p == 0) x /= p, ++y;
      if (y != 0) res.emplace_back(p,y);
    }
    if (x != 1) res.emplace_back(x,1);
    return res;
  }
};

int query(int i, int j) {
    cout << "? " << i+1 << " " << j+1 << endl;
    int k;
    cin >> k;
    return k; 
};

int main()
{
    int n;
    cin >> n;
    Sieve s(n);
    int p = 0;
    for(int i=1;i<=n;i++){
        if(s.isPrime(i)) p = i;
    }
    dbg(p);
    vector<int> ans(n);
    rep(i,n){
        ans[i] = query(0,i);
    }
    int cnt = 0;
    rep(i,n){
        if(ans[i]%p==0) cnt++;
    }
    vector<int> a(n), b(n);
    if(cnt==1){
        int j = -1;
        rep(i,n){  
            if(ans[i]%p==0) {
                b[i] = p;
                j = i;
            }
        }
        a[0] = ans[0]/p;
        for(int i=1;i<n;i++){
            int ret = query(i,j);
            a[i] = ret/p;
        }
        vector<int> vs;
        for(int i=0;i<n;i++){
            if(a[i]==1) vs.push_back(i);
        }  
        int ret = query(vs[0], (j+1)%n);
        int k = -1;
        if(ret%p==0) {
            a[vs[0]] = p;
            a[vs[1]] = 1;
            k = vs[1];
        }
        else{
            a[vs[0]] = 1;
            a[vs[1]] = p;
            k = vs[0];
        }
        for(int i=0;i<n;i++){
            int ret = query(k,i);
            b[i] = ret;
        }
    }
    else{
        rep(i,n){
            b[i] = ans[i]/p;
        }
        vector<int> vs;
        rep(i,n){
            if(b[i]==1) vs.push_back(i);
        }
        int ret = query(1, vs[0]);
        int k = 0;
        if(ret%p==0){
            b[vs[0]] = p;
            b[vs[1]] = 1;
            k = vs[1];
        }
        else{
            b[vs[0]] = 1;
            b[vs[1]] = p;
            k = vs[0];
        }
        rep(i,n){
            a[i] = query(i, k);
        }
    }
    cout << "! ";
    rep(i,n){
        cout << a[i] << " ";
    }
    rep(i,n){
        cout << b[i] << " ";
    }
    cout << endl;
    return 0;
}
0