結果

問題 No.144 エラトステネスのざる
ユーザー FromBooskaFromBooska
提出日時 2023-02-15 13:24:06
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 466 bytes
コンパイル時間 216 ms
コンパイル使用メモリ 82,532 KB
実行使用メモリ 97,984 KB
最終ジャッジ日時 2024-07-17 23:08:59
合計ジャッジ時間 4,351 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 33 ms
57,600 KB
testcase_01 AC 34 ms
51,712 KB
testcase_02 AC 33 ms
51,840 KB
testcase_03 AC 34 ms
51,712 KB
testcase_04 AC 35 ms
51,968 KB
testcase_05 AC 35 ms
52,352 KB
testcase_06 AC 40 ms
59,008 KB
testcase_07 AC 39 ms
59,264 KB
testcase_08 AC 39 ms
59,648 KB
testcase_09 AC 39 ms
59,392 KB
testcase_10 AC 43 ms
59,264 KB
testcase_11 AC 39 ms
59,904 KB
testcase_12 AC 38 ms
59,240 KB
testcase_13 TLE -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

N, exe_prob = map(float, input().split())
N = int(N)

def divisor_count(n):
    cnt = 0
    i = 1
    while i*i < n:
        if n % i == 0:
            cnt += 2
        i += 1
    if i*i == n:
        cnt += 1
    return cnt

is_prime = [1] * (N+1)
is_prime[0], is_prime[1] = 0, 0

for i in range(2, N+1):
    is_prime[i] = pow(1-exe_prob, divisor_count(i)-2)
    #is_prime[i] = ((1-exe_prob)**(divisor_count(i)-2))
            
#print(is_prime)
print(sum(is_prime))
0