結果

問題 No.144 エラトステネスのざる
ユーザー FromBooska
提出日時 2023-02-15 13:21:21
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 411 bytes
コンパイル時間 137 ms
コンパイル使用メモリ 82,492 KB
実行使用メモリ 98,108 KB
最終ジャッジ日時 2024-07-17 23:06:56
合計ジャッジ時間 4,259 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 1 -- * 6
権限があれば一括ダウンロードができます

ソースコード

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] = ((1-exe_prob)**(divisor_count(i)-2))
            
#print(is_prime)
print(sum(is_prime))
0