結果

問題 No.3461 Min GCD
コンテスト
ユーザー くらげ
提出日時 2026-01-01 10:48:36
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 802 ms / 2,000 ms
コード長 1,531 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,463 ms
コンパイル使用メモリ 228,292 KB
実行使用メモリ 219,264 KB
最終ジャッジ日時 2026-02-28 13:02:20
合計ジャッジ時間 7,632 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge7
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using P = pair<ll,ll>;
#define rep(i,n) for(int (i) = 0; i<(n); i++)

const int N = 100000;
const ll inf = 1e18;

int main(){
    int n; ll k;
    cin >> n >> k;
    vector<ll> a(n), b(n);
    rep(i,n) cin >> a[i];
    rep(i,n) cin >> b[i];

    // 制約
    assert(1<=n && n<=N);
    assert(1<=k && k<=inf);
    rep(i,n) assert(1<=a[i] && a[i]<=N);
    rep(i,n) assert(1<=b[i] && b[i]<=N);

    // div[i] : i の約数
    vector<vector<int>> div(N+10);
    for(int i=1; i<=N+5; i++){
        int j = i;
        while(j<=N+5){
            div[j].push_back(i);
            j += i;
        }
    }

    // gcd を div の値にするのに必要な操作回数
    vector<vector<P>> min_op(n);
    rep(i,n){
        ll min_cost = inf;
        for(int j = div[a[i]].size()-1; 0<=j; j--){
            ll p = div[a[i]][j];
            ll cost = (b[i]%p==0? 0 : p-b[i]%p);
            if(cost<min_cost){
                min_op[i].push_back(P(p,cost));
                min_cost = cost;
            }
        }
        sort(min_op[i].begin(), min_op[i].end());
    }

    // 二分探索
    int ok = 1, ng = *min_element(a.begin(),a.end())+1;
    while(1<abs(ok-ng)){
        int mid = (ok+ng)/2;
        ll sum = 0; // 加算が必要な回数
        rep(i,n) sum += (*lower_bound(min_op[i].begin(), min_op[i].end(), P(mid,0))).second;
        if(sum<=k) ok = mid;
        else ng = mid;
    }

    cout << ok << endl;

    //cerr << clock()/1000 << "ms" << endl;

}
0