結果

問題 No.3461 Min GCD
コンテスト
ユーザー MM
提出日時 2026-02-28 15:01:43
言語 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,712 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 6,758 ms
コンパイル使用メモリ 382,784 KB
実行使用メモリ 19,656 KB
最終ジャッジ日時 2026-02-28 15:01:57
合計ジャッジ時間 13,495 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample TLE * 1 -- * 2
other -- * 21
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#include<atcoder/all>
#define chmin(x,y) (x) = min((x),(y))
#define chmax(x,y) (x) = max((x),(y))
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define vec vector
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
#define pb push_back 
#define eb emplace_back

using namespace std;
using namespace atcoder;
using ll = long long;
using ld = long double;
const ll mod = 998244353;
using mint = modint998244353;
const vector<int> dx = {1,0,-1,0}, dy = {0,1,0,-1};
// using Graph = vector<vector<pair<int,ll>>>;
using Graph = vector<vector<int>>;

vector<int> divisors(int n){
  vector<int> res;
  for(int i = 1; i * i <= n; i++){
    if(n % i == 0){
      res.pb(i);
      if(i * i != n)
        res.pb(n / i);
    }
  }
  sort(all(res));
  return res;
}

int main(){
  // input + prep
  int N; ll K;
  cin >> N >> K;
  vector<int> A(N),B(N);
  vector<vector<int>> D(N);
  rep(i,N){
    cin >> A[i];
    D[i] = divisors(A[i]);
  }
  rep(i,N) cin >> B[i];

  // prep
  vector<vector<int>> E(N),F(N);
  rep(i,N){
    int M = (int)D[i].size();
    E[i].reserve(M * sizeof(int));
    F[i].reserve(M * sizeof(int));
    rep(j,M){
      E[i][j] = (B[i] % D[i][j] == 0 ? 0 : D[i][j] - B[i] % D[i][j]);
      F[i][j] = (j > 0 ? min(F[i][j-1],E[i][j]) : E[i][j]);
    }
  }

  // solve
  int ok = 1, ng = 100001;
  while(ng - ok > 1){
    int mid = (ok + ng);
    bool can = 1;
    ll cnt = 0;

    rep(i,N){
      int pos = lower_bound(all(D[i]),mid) - D[i].begin();
      if(pos == (int)D[i].size())
        can = 0;
      else
        cnt += F[i][pos];
    }

    if(cnt > K) can = 0;

    if(can) ok = mid;
    else ng = mid;
  }
  
  // output
  cout << ok << endl;
}
0