結果

問題 No.144 エラトステネスのざる
ユーザー buey_tbuey_t
提出日時 2023-04-29 14:44:41
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 427 ms / 2,000 ms
コード長 1,642 bytes
コンパイル時間 628 ms
コンパイル使用メモリ 86,908 KB
実行使用メモリ 157,624 KB
最終ジャッジ日時 2023-08-11 17:13:19
合計ジャッジ時間 9,877 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 336 ms
157,000 KB
testcase_01 AC 347 ms
157,036 KB
testcase_02 AC 364 ms
157,044 KB
testcase_03 AC 336 ms
157,232 KB
testcase_04 AC 336 ms
156,956 KB
testcase_05 AC 339 ms
156,992 KB
testcase_06 AC 334 ms
157,152 KB
testcase_07 AC 355 ms
156,956 KB
testcase_08 AC 347 ms
157,108 KB
testcase_09 AC 353 ms
156,952 KB
testcase_10 AC 345 ms
157,224 KB
testcase_11 AC 355 ms
157,008 KB
testcase_12 AC 347 ms
157,052 KB
testcase_13 AC 381 ms
157,104 KB
testcase_14 AC 427 ms
157,592 KB
testcase_15 AC 425 ms
157,624 KB
testcase_16 AC 411 ms
157,552 KB
testcase_17 AC 425 ms
157,536 KB
testcase_18 AC 420 ms
157,360 KB
testcase_19 AC 387 ms
157,312 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import sqrt,sin,cos,tan,ceil,radians,floor,gcd,exp,log,log10,log2,factorial,fsum
from heapq import heapify, heappop, heappush
from bisect import bisect_left, bisect_right
from copy import deepcopy
import copy
import random
from random import randrange
from collections import deque,Counter,defaultdict
from itertools import permutations,combinations
from decimal import Decimal,ROUND_HALF_UP
#tmp = Decimal(mid).quantize(Decimal('0'), rounding=ROUND_HALF_UP)
from functools import lru_cache, reduce
#@lru_cache(maxsize=None)
from operator import add,sub,mul,xor,and_,or_,itemgetter
import sys
input = sys.stdin.readline
# .rstrip()
INF = 10**18
mod1 = 10**9+7
mod2 = 998244353

#DecimalならPython
#再帰ならPython!!!!!!!!!!!!!!!!!!!!!!!!!!


'''
処理が行われた場合どうなる
処理が行われなかった場合は
処理が行われなかった場合、素数の数が一つ増えるし、その先でも処理が行われるかどうかの選択がある
行われた場合、その倍数はすべて候補から外れる
逆なんじゃね?
約数の数を数えておく。それで、p^約数の数をしのげば素数の加わる
これだな

'''

def sieve(n):
    is_prime = [0 for _ in range(n+1)]
    is_prime[0] = 0
    is_prime[1] = 0

    for i in range(2, n+1):
        j = 2 * i
        while j <= n:
            is_prime[j] += 1
            j += i
    table = [ i for i in range(1, n+1) if is_prime[i]]
    return is_prime, table

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

prime,t = sieve(10**6+1)

ans = 0
for i in range(2,N+1):
    ans += pow(1-P,prime[i])

print(ans)
0