結果

問題 No.1653 Squarefree
ユーザー rlangevinrlangevin
提出日時 2023-10-25 08:57:17
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 594 ms / 2,000 ms
コード長 762 bytes
コンパイル時間 217 ms
コンパイル使用メモリ 81,940 KB
実行使用メモリ 92,288 KB
最終ジャッジ日時 2024-09-25 04:05:17
合計ジャッジ時間 16,707 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 594 ms
91,752 KB
testcase_01 AC 84 ms
72,448 KB
testcase_02 AC 86 ms
72,448 KB
testcase_03 AC 88 ms
74,880 KB
testcase_04 AC 83 ms
72,448 KB
testcase_05 AC 81 ms
72,704 KB
testcase_06 AC 83 ms
72,576 KB
testcase_07 AC 81 ms
72,832 KB
testcase_08 AC 83 ms
72,448 KB
testcase_09 AC 85 ms
72,576 KB
testcase_10 AC 82 ms
72,704 KB
testcase_11 AC 276 ms
87,552 KB
testcase_12 AC 279 ms
87,552 KB
testcase_13 AC 529 ms
92,032 KB
testcase_14 AC 506 ms
91,904 KB
testcase_15 AC 563 ms
92,224 KB
testcase_16 AC 558 ms
92,032 KB
testcase_17 AC 546 ms
91,776 KB
testcase_18 AC 581 ms
91,904 KB
testcase_19 AC 551 ms
91,904 KB
testcase_20 AC 518 ms
91,904 KB
testcase_21 AC 564 ms
91,648 KB
testcase_22 AC 541 ms
92,160 KB
testcase_23 AC 543 ms
91,520 KB
testcase_24 AC 517 ms
92,288 KB
testcase_25 AC 519 ms
91,908 KB
testcase_26 AC 522 ms
91,820 KB
testcase_27 AC 464 ms
91,516 KB
testcase_28 AC 471 ms
91,680 KB
testcase_29 AC 494 ms
91,672 KB
testcase_30 AC 512 ms
91,840 KB
testcase_31 AC 483 ms
91,648 KB
testcase_32 AC 515 ms
92,088 KB
testcase_33 AC 484 ms
91,676 KB
testcase_34 AC 486 ms
91,740 KB
testcase_35 AC 520 ms
92,084 KB
testcase_36 AC 79 ms
73,448 KB
testcase_37 AC 81 ms
72,448 KB
testcase_38 AC 85 ms
72,848 KB
testcase_39 AC 78 ms
74,056 KB
testcase_40 AC 79 ms
73,880 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def isqrt(n):
    yes = 1
    no = 10 ** 9 + 5
    while no - yes != 1:
        mid = (yes + no)//2
        if mid * mid <= n:
            yes = mid
        else:
            no = mid
    return yes * yes == n

L, R = map(int, input().split())
D = list(range(L, R + 1))
M = 10 ** 6 + 5
Prime = [1] * M
for i in range(2, M):
    if not Prime[i]:
        continue
    for j in range(2 * i, M, i):
        Prime[j] = 0
    start = (L + i - 1)//i * i
    for j in range(start - L, R - L + 1, i):
        if D[j] % (i * i) == 0:
            D[j] = 0
        elif D[j] % i == 0:
            D[j] //= i

ans = 0
for i in range(R - L + 1):
    if D[i] == 1:
        ans += 1
    elif D[i] == 0:
        continue    
    elif not isqrt(D[i]):
        ans += 1

print(ans)
0