結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー FromBooska
提出日時 2023-02-15 13:35:29
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 108 ms / 2,000 ms
コード長 1,124 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 161 ms
コンパイル使用メモリ 85,428 KB
実行使用メモリ 177,992 KB
最終ジャッジ日時 2026-04-01 14:23:16
合計ジャッジ時間 2,046 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

# 落とし穴がある。たとえば8
# 2でも4でも倍数としてひっかからなければ8は生き残る
# 2で8がひっかからない確率は1-exe_prob
# 4で8がひっかからないかどうかを考えるときに4が生きているか死んでいるかの場合分けをしたくなる
# しかし4が死んでいれば8も死んでいる、つまり完全相関
# つまり4まで来たときに8が生き残っていれば4も生き残っている
# 8が生き残るかは2と4の2回のテストをクリアすればいい、4の場合分けは必要ない
# つまり各数の約数の個数だけが必要となる
# (1-exe_prob)**約数個数は時間かかるのでpow
# さらに約数個数は下から積み上げる、エラトステネスの容量、そうしないとTLE

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

div_count = [0]*(N+1)

for p in range(2, N+1):
    for q in range(p*2, N+1, p):
        div_count[q] += 1

ans = 0
for i in range(2, N+1):
    if div_count[i] == 0:
        ans += 1
    else:
        ans += pow(1-exe_prob, div_count[i])

print(ans)

#print(div_count)
0