結果

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

ソースコード

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;
        if(n==1){
            cout << "!" << " " << 1 << " " << 1 << endl;
            return 0;
        }
        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[j]/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);
            }  
            dbg(vs);
            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;
    }
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0