結果

問題 No.3461 Min GCD
コンテスト
ユーザー koheijkt
提出日時 2026-02-28 15:10:27
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
TLE  
実行時間 -
コード長 1,969 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 219 ms
コンパイル使用メモリ 78,172 KB
実行使用メモリ 203,360 KB
最終ジャッジ日時 2026-02-28 15:10:50
合計ジャッジ時間 13,685 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys, math
sys.setrecursionlimit(10**8)
sys.set_int_max_str_digits(0)
INF = 1e18
MOD = 998244353
from bisect import bisect_left, bisect_right
from collections import deque, defaultdict, Counter
from itertools import product, combinations, permutations, groupby, accumulate
from heapq import heapify, heappop, heappush
input = sys.stdin.readline
def I():   return input().rstrip()
def II():  return int(input().rstrip())
def IS():  return input().rstrip().split()
def MII(): return map(int, input().rstrip().split())
def LI():  return list(input().rstrip())
def TII(): return tuple(map(int, input().rstrip().split()))
def LII(): return list(map(int, input().rstrip().split()))
def LSI(): return list(map(str, input().rstrip().split()))
def GMI(): return list(map(lambda x: int(x) - 1, input().rstrip().split()))
def kiriage(a, b): return (a+b-1)//b

N, K = MII()
A = LII()

def make_divisors(n):
    lower_divisors , upper_divisors = [], []
    i = 1
    while i*i <= 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]

A_divs = []
for a in A:
    A_divs.append(make_divisors(a)[::-1])
B = LII()

def solve(m): # 全ペアの GCD を m以上にできますか？
    total_cnt = 0
    for i in range(N):
        cnt = K + 1

        # a の約数を大きいものから見る
        for div in A_divs[i]:
            if div < m:
                break
            # B[i] を div の約数になるまでいくらコストが必要か？
            cost = B[i]%div
            if cost != 0:
                cost = div - cost
            # 最小か？            
            if cost < cnt:
                cnt = cost

        total_cnt += cnt

    return total_cnt <= K

ok = 1
ng = 10**5 + 1
while abs(ok - ng) != 1:
    mid = (ok + ng)//2
    if solve(mid):
        ok = mid
    else:
        ng = mid

print(ok)
0