結果

問題 No.3461 Min GCD
コンテスト
ユーザー Guran08
提出日時 2026-02-28 16:07:44
言語 C++23
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
TLE  
実行時間 -
コード長 1,981 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 7,856 ms
コンパイル使用メモリ 392,180 KB
実行使用メモリ 7,848 KB
最終ジャッジ日時 2026-02-28 16:07:56
合計ジャッジ時間 11,200 ms
ジャッジサーバーID
（参考情報） 		judge7 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 3 TLE * 1 -- * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
//using mint = modint1000000007;
typedef long long ll;
typedef pair<ll, ll> P;
typedef tuple<ll, ll, ll> T;
template<typename T>bool chmax(T& a, const T& b) { if (a < b) { a = b;return true; } else { return false; } }
template<typename T>bool chmin(T& a, const T& b) { if (a > b) { a = b;return true; } else { return false; } }
template<class... Args> void dbg(Args&&... args) { ((cout << args << ' '), ...);cout << '\n'; }

const int di[] = { -1,0,1,0 };
const int dj[] = { 0,-1,0,1 };
const long long INF = 6000000000000000000;
const int inf = 1001001001;


int main(void) {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int n; cin >> n;
    ll k; cin >> k;
    vector<int> a(n, 0), b(n, 0);
    for (int i = 0; i < n; i++) {
        cin >> a[i];
    }
    for (int i = 0; i < n; i++) {
        cin >> b[i];
    }
    unordered_map<ll, set<int>> c;
    for (int j = 0; j < n; j++) {
        if (c.contains(a[j]))continue;
        ll g = gcd(a[j], b[j]);
        for (long long i = g; i * i <= a[j]; ++i) {
            if (a[j] % i == 0) {
                c[j].insert(i);
                // 重複しないならば i の相方である N/i も push
                if (a[j] / i != i)c[j].insert(a[j] / i);
            }
        }
    }

    priority_queue<P, vector<P>, greater<P>> pq;
    for (int i = 0; i < n; i++) {
        pq.emplace(gcd(a[i], b[i]), i);
    }
    int ans = inf;
    while (pq.size()) {
        auto [g, p] = pq.top();pq.pop();
        // dbg(g, p, k);
        if (g == a[p])continue;
        auto it = c[p].find(g);
        it++;
        int x = *it;
        int cnt = b[p] % x;
        if (k - (x - cnt) >= 0)k -= x - cnt;
        else break;

        b[p] += x - cnt;
        pq.emplace(gcd(a[p], b[p]), p);
    }
    for (int i = 0; i < n; i++) {
        chmin(ans, gcd(a[i], b[i]));
    }
    cout << ans << "\n";
}
0