結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー yuki2006
提出日時 2015-02-05 23:48:09
言語 PyPy2
(7.3.15)
結果
WA  
実行時間 -
コード長 492 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 345 ms
コンパイル使用メモリ 76,928 KB
実行使用メモリ 100,352 KB
最終ジャッジ日時 2024-06-23 08:53:56
合計ジャッジ時間 7,817 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample WA * 3
other WA * 3 TLE * 1 -- * 13
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

# -*- coding: utf-8 -*-

import random
import time


N, p = map(float, raw_input().split())
N = int(N)

QQ = 10 ** 6
total = 0
q = 0

# モンテカルロ法

for q in xrange(QQ):

    prime = [True] * (N + 1)
    prime[0] = prime[1] = False

    for v in xrange(2, N + 1):
        if not prime[v]:
            continue
        if random.random() <= p:
            for i in xrange(v + v, N + 1, v):
                prime[i] = False

    total += sum(prime)
    q += 1
print q, 1.0 * total / q
0