結果

問題 No.2849 Birthday Donuts
ユーザー ecotteaecottea
提出日時 2024-06-20 02:05:46
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 3,626 ms / 6,000 ms
コード長 1,704 bytes
コンパイル時間 228 ms
コンパイル使用メモリ 82,464 KB
実行使用メモリ 83,744 KB
最終ジャッジ日時 2024-07-01 01:48:22
合計ジャッジ時間 70,343 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 97 ms
74,616 KB
testcase_01 AC 97 ms
74,452 KB
testcase_02 AC 3,512 ms
83,260 KB
testcase_03 AC 3,566 ms
82,908 KB
testcase_04 AC 3,602 ms
83,600 KB
testcase_05 AC 3,603 ms
83,688 KB
testcase_06 AC 3,626 ms
83,376 KB
testcase_07 AC 3,593 ms
83,608 KB
testcase_08 AC 3,570 ms
83,744 KB
testcase_09 AC 3,598 ms
83,096 KB
testcase_10 AC 3,585 ms
83,124 KB
testcase_11 AC 3,589 ms
83,156 KB
testcase_12 AC 3,256 ms
83,648 KB
testcase_13 AC 3,394 ms
83,688 KB
testcase_14 AC 3,417 ms
83,604 KB
testcase_15 AC 3,414 ms
83,688 KB
testcase_16 AC 3,580 ms
83,516 KB
testcase_17 AC 3,266 ms
83,304 KB
testcase_18 AC 3,569 ms
83,192 KB
testcase_19 AC 3,354 ms
83,084 KB
testcase_20 AC 3,395 ms
83,196 KB
testcase_21 AC 599 ms
83,592 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math

n = 2 * 10**5 + 10

pf = [0] * n
for i in range(n):
    pf[i] = i

for i in range(2, n):
    if pf[i] == i:
        j = 2 * i
        while j < n:
            pf[j] = i
            j = j + i
# print(pf)

phi = [1] * n
for i in range(2, n):
    x = i
    pp = -1
    while x > 1:
        p = pf[x]
        if p == pp:
            phi[i] *= p
        else:
            phi[i] *= p - 1
        pp = p
        x //= p
# print(phi)

acc = [0] * (n + 1)
for i in range(2, n):
    acc[i] = acc[i - 1] + phi[i]


import re
import sys

content = input()
pattern = r'^(\d+)$'
result = re.match(pattern, content)

if not result:
    sys.exit("format error.")

T = int(content)

prv = 0

for t in range(T):
    content = input()
    pattern = r'^(\d+) (\d+)$'
    result = re.match(pattern, content)

    if not result:
        sys.exit("format error.")

    l, r = map(int, content.split())

    if not (0 <= l and l <= 10**18):
        sys.exit("alpha")
    if not (0 <= r and r <= 10**18):
        sys.exit("beta")

    l ^= prv
    r ^= prv

    if not (2 <= l and l <= r and r <= 2 * 10**5):
        sys.exit("L, R")
    
    l -= 1

    res = 0
    
    sqrtR = int(math.sqrt(l + r))

    for i in range(2, sqrtR + 1):
        if l // i != r // i:
            res += phi[i]
    
    i1 = 0
    i2 = 1
    q_prv = r
    while q_prv > sqrtR:
        ql = l // (i1 + 1)
        qr = r // (i2 + 1)
        q_max = max(ql, qr)
        if i1 != i2:
            res += acc[q_prv] - acc[max(q_max, sqrtR)]
            # print(i1, i2, max(q_max, sqrtR), q_prv)
        
        q_prv = q_max
        if ql >= qr:
            i1 += 1
        else:
            i2 += 1
    
    print(res)
    
    prv = res
0