結果

問題 No.3054 Modulo Inequalities
コンテスト
ユーザー suisen
提出日時 2025-01-18 16:17:59
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 534 ms / 2,500 ms
コード長 787 bytes
コンパイル時間 391 ms
コンパイル使用メモリ 82,320 KB
実行使用メモリ 108,196 KB
最終ジャッジ日時 2025-02-02 13:29:26
合計ジャッジ時間 11,187 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

# correct

M = 10 ** 6 + 10

n = int(input())

W = [0] * (2 * M + 1)

for _ in range(n):
    a, b = map(int, input().split())
    W[M + a + 1] += 1
    W[M + b + 1] -= 1
    W[M + a - b + 1] -= 1

# accumulate
for i in range(2 * M):
    W[i + 1] += W[i]

S = [0] * (M + 1)
S[0] = W[M]
for i in range(1, M + 1):
    S[i] = W[M - i] + W[M + i]

del W

# zeta transform
is_prime = [True] * (M + 1)
for p in range(2, M + 1):
    if not is_prime[p]:
        continue
    q = M // p
    while q >= p:
        is_prime[p * q] = False  # sieve
        S[q] += S[p * q]
        q -= 1
    while q >= 1:
        S[q] += S[p * q]
        q -= 1

max_f, argmax = -1, 0

for m in range(1, M + 1):
    f_m = -n * (M // m) - S[0] - S[m]
    if f_m > max_f:
        max_f, argmax = f_m, m
print(argmax)
0