結果

問題 No.3093 Please GCD
ユーザー eggkid6eggkid6
提出日時 2022-04-01 22:18:16
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 3,834 bytes
コンパイル時間 2,704 ms
コンパイル使用メモリ 214,344 KB
実行使用メモリ 25,332 KB
平均クエリ数 227.00
最終ジャッジ日時 2024-11-20 09:51:04
合計ジャッジ時間 9,242 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using V = vector<ll>;
using VV = vector<V>;
using VVV = vector<VV>;
using VVVV = vector<VVV>;
using VVVVV = vector<VVVV>;
using VVVVVV = vector<VVVVV>;
using VS = vector<string>;
using VB = vector<bool>;
using VVB = vector<VB>;
using P = pair<ll,ll>;
using M = map<ll,ll>;
using Q = queue<ll>;
using PQ = priority_queue<ll>;
using PQG = priority_queue<ll,V,greater<ll>>;
using S = set<ll>;
using VP = vector<P>;
const ll MOD = 1000000007;
const ll mod = 998244353;
const ll INF = 1LL << 60;
#define rep(i,n) for(ll i = 0; i < n; i++)
#define rep2(i,s,n) for(ll i = s; i < n; i++)
#define per(i,n) for(ll i = n-1; i >= 0; i--)
#define per2(i,s,n) for(ll i = n-1; i >= s; i--)
#define all(x) (x).begin(),(x).end()
#define rall(x) (x).rbegin(),(x).rend()
#define fi first
#define se second
#define pb push_back
#define pf push_front
#define ppb pop_back
#define ppf pop_front
#define eb emplace_back
#define lb lower_bound
#define ub upper_bound
template<class T>bool chmin(T&a, const T&b){if(a>b){a=b;return true;}return false;}
template<class T>bool chmax(T&a, const T&b){if(a<b){a=b;return true;}return false;}
template<class T>void Vin(vector<T>&a){rep(i,(ll)a.size())cin>>a[i];}
template<class T>void VVin(vector<vector<T>>&a){rep(i,(ll)a.size())Vin(a[i]);}
template<class T>void Vout(const vector<T>&a){rep(i,(int)a.size()-1)cout<<a[i]<<" ";cout<<a.back()<<'\n';}
template<class T>void Voutl(const vector<T>&a){rep(i,(int)a.size())cout<<a[i]<<'\n';}
ll power(ll a,ll b,const ll&M){a%=M;ll res=1;while(b>0){if(b&1)res=res*a%M;a=a*a%M;b>>=1;}return res;}
const ll H[4] = {0,1,0,-1}, W[4] = {1,0,-1,0};

struct Sieve {
    long long n;
    vector<long long> f, primes;
    Sieve(long long n): n(n), f(n+1) {
        f[0] = f[1] = -1;
        for(long long i = 2; i <= n; i++) {
            if(f[i]) continue;
            primes.push_back(i);
            f[i] = i;
            for(long long j = i*i; j <= n; j += i) if(!f[j]) f[j] = i;
        }
    }
    bool isPrime(long long x) {
        return f[x] == x;
    }
    vector<long long> factorList(long long x) {
        if(!x) return {};
        vector<long long> res;
        while(x != 1) {
            res.push_back(f[x]);
            x /= f[x];
        }
        return res;
    }
    vector<pair<long long,long long>> factor(long long x) {
        vector<long long> fl = factorList(x);
        if(fl.size() == 0) return {};
        vector<pair<long long,long long>> res = {{fl[0],0}};
        for(long long p : fl) {
            if(p == res.back().fi) res.back().se++;
            else res.emplace_back(p,1);
        }
        return res;
    }
    vector<pair<long long,long long>> bigfactor(long long x) {
        if(!x) return {};
        vector<pair<long long,long long>> res;
        for(long long p : primes) {
            long long y = 0;
            while(x%p == 0) x /= p, y++;
            if(y) res.emplace_back(p,y);
            if(x == 1) break;
        }
        if(x != 1) res.emplace_back(x,1);
        return res;
    }
};

int main() {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    int n;
    cin >> n;
    int ans = 1;
    Sieve s(n + 1);
    ll sz = s.primes.size();
    ll y = 1;
    S st;
    per(i, sz) {
        y *= s.primes[i];
        if(y > 100000000000000 || i == 0) {
            cout << "? " << y << endl;
            ll z;
            cin >> z;
            VP p = s.bigfactor(z);
            for(auto q : p) {
                st.insert(q.fi);
                ans *= q.fi;
            }
            y = 1;
        }
    }
    for(ll x : st) {
        if(ans * x > n) break;
        ll z = 1;
        while(z <= n) z *= x;
        z /= x;
        cout << "? " << z << endl;
        ll t;
        cin >> t;
        ans *= t / x;
    }
    cout << "! " << ans << endl;
}
0