結果

問題 No.3461 Min GCD
コンテスト
ユーザー i_taku
提出日時 2026-03-03 08:36:57
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 716 ms / 2,000 ms
コード長 1,071 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 490 ms
コンパイル使用メモリ 85,444 KB
実行使用メモリ 302,384 KB
最終ジャッジ日時 2026-03-03 08:37:04
合計ジャッジ時間 6,052 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

from bisect import bisect_left

def make_divisors(n):
    lower_divisors, upper_divisors = [], []
    i = 1
    while i**2 <= n:
        if n % i == 0:
            lower_divisors.append(i)
            if i != n // i:
                upper_divisors.append(n // i)
        i += 1
    return lower_divisors + upper_divisors[::-1]

N, K = map(int, input().split())
A = list(map(int, input().split()))
B = list(map(int, input().split()))
divs = [make_divisors(a) for a in A]
cums = [[1 << 60] * len(divs[i]) for i in range(N)]
for i in range(N):
    div = divs[i]
    cum = cums[i]
    for k in range(len(div)):
        d = div[k]
        cum[k] = (B[i] + d - 1) // d * d - B[i]
    for k in range(len(div) - 1, 0, -1):
        cum[k - 1] = min(cum[k], cum[k - 1])

ok = 1
ng = 10**11
while ng - ok > 1:
    mid = (ok + ng) // 2
    cnt = 0
    for i in range(N):
        idx = bisect_left(divs[i], mid)
        if idx == len(divs[i]):
            cnt = K + 1
            break
        cnt += cums[i][idx]
    if cnt <= K:
        ok = mid
    else:
        ng = mid
print(ok)
0