結果

問題 No.144 エラトステネスのざる
ユーザー FromBooska
提出日時 2023-02-15 13:35:29
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 175 ms / 2,000 ms
コード長 1,124 bytes
コンパイル時間 188 ms
コンパイル使用メモリ 82,372 KB
実行使用メモリ 84,236 KB
最終ジャッジ日時 2024-07-17 23:16:43
合計ジャッジ時間 2,752 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 落とし穴がある。たとえば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