結果

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

ソースコード

diff #
raw source code 

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

import random
import time

start = time.clock()

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

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

# モンテカルロ法

while time.clock() - start < 0.8:

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

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

    total += sum(prime)
    q += 1

print 1.0 * total / q
0