結果
問題 | No.2849 Birthday Donuts |
ユーザー | 👑 AngrySadEight |
提出日時 | 2024-06-11 01:35:57 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,450 bytes |
コンパイル時間 | 170 ms |
コンパイル使用メモリ | 82,292 KB |
実行使用メモリ | 85,436 KB |
最終ジャッジ日時 | 2024-07-01 01:43:16 |
合計ジャッジ時間 | 7,434 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | AC | 138 ms
83,140 KB |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | RE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | RE | - |
testcase_17 | RE | - |
testcase_18 | RE | - |
testcase_19 | RE | - |
testcase_20 | RE | - |
testcase_21 | WA | - |
ソースコード
def my_gcd(a, b): ret = 0 if a < b: a, b = b, a if b == 0: ret = a else: ret = my_gcd(b, a % b) return ret divisor = [-1 for _ in range(300002)] for i in range(2, 300001): for j in range(2, 300000 // i + 1): divisor[j * i] = j phi = [0 for _ in range(300002)] for i in range(2, 300001): now = i div = -1 tot = i while True: if divisor[now] == -1: if now != div: tot = tot * (now - 1) tot //= now break if div != divisor[now]: div = divisor[now] tot = tot * (div - 1) tot //= div now //= div phi[i] = tot phi_sum = [0 for _ in range(300002)] for i in range(2, 300001): phi_sum[i] = phi_sum[i - 1] + phi[i] T = int(input()) for _ in range(T): L, R = map(int, input().split()) ans = 0 for i in range(2, 401): k1 = (L - 1) // i k2 = R // i if k1 < k2: ans += phi[i] left = 0 right = 1 prev = R while True: if R < 400: break nl = (L - 1) // (left + 1) nr = R // (right + 1) val = max(nl, nr) if left != right: ans += (phi_sum[prev] - phi_sum[max(val, 400)]) if val <= 400: break prev = val if nl >= nr: left += 1 else: right += 1 print(ans)