結果

問題 No.144 エラトステネスのざる
ユーザー yuki2006yuki2006
提出日時 2015-02-05 23:48:09
言語 PyPy2
(7.3.15)
結果
WA  
実行時間 -
コード長 492 bytes
コンパイル時間 345 ms
コンパイル使用メモリ 76,928 KB
実行使用メモリ 100,352 KB
最終ジャッジ日時 2024-06-23 08:53:56
合計ジャッジ時間 7,817 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 TLE -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

# -*- 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