結果

問題 No.3461 Min GCD
コンテスト
ユーザー 👑 loop0919
提出日時 2026-02-28 15:19:33
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
MLE  
実行時間 -
コード長 1,146 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 332 ms
コンパイル使用メモリ 78,012 KB
実行使用メモリ 556,732 KB
最終ジャッジ日時 2026-02-28 15:20:31
合計ジャッジ時間 11,295 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge7
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18 MLE * 1 -- * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import math
from bisect import bisect_left
from functools import cache

LIMIT = 10**5 + 1
INF = 10**9


@cache
def factorize(n):
    factorized = []
    for i in range(1, math.isqrt(n) + 1):
        if n % i == 0:
            factorized.append(i)
            factorized.append(n // i)

    if math.isqrt(n) ** 2 == n:
        factorized.pop()

    return factorized


N, K = [int(s) for s in input().split()]
A = [int(s) for s in input().split()]
B = [int(s) for s in input().split()]

need_op = [dict() for _ in range(N)]

for i, (a, b) in enumerate(zip(A, B)):
    for fac in factorize(a):
        need_op[i][fac] = (fac - b) % fac

    min_cnt = INF

    for key in sorted(need_op[i], reverse=True):
        min_cnt = min(need_op[i][key], min_cnt)
        need_op[i][key] = min_cnt

keys = [list(sorted(need_op[i])) for i in range(N)]


def is_ok(mid):
    cnt = 0

    for i, key in enumerate(keys):
        idx = bisect_left(key, mid)
        cnt += need_op[i][key[idx]]

    return cnt <= K


ok, ng = 0, min(A) + 1
while abs(ok - ng) > 1:
    mid = (ok + ng) // 2

    if is_ok(mid):
        ok = mid
    else:
        ng = mid

print(ok)
0