結果

問題 No.144 エラトステネスのざる
ユーザー buey_tbuey_t
提出日時 2023-04-29 14:44:41
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 314 ms / 2,000 ms
コード長 1,642 bytes
コンパイル時間 250 ms
コンパイル使用メモリ 82,124 KB
実行使用メモリ 163,712 KB
最終ジャッジ日時 2024-04-29 08:53:39
合計ジャッジ時間 6,836 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 254 ms
163,200 KB
testcase_01 AC 252 ms
163,200 KB
testcase_02 AC 250 ms
163,372 KB
testcase_03 AC 248 ms
163,200 KB
testcase_04 AC 250 ms
163,328 KB
testcase_05 AC 314 ms
163,276 KB
testcase_06 AC 251 ms
163,272 KB
testcase_07 AC 251 ms
163,200 KB
testcase_08 AC 241 ms
163,200 KB
testcase_09 AC 244 ms
163,328 KB
testcase_10 AC 242 ms
163,328 KB
testcase_11 AC 243 ms
163,328 KB
testcase_12 AC 248 ms
163,276 KB
testcase_13 AC 264 ms
163,456 KB
testcase_14 AC 306 ms
163,584 KB
testcase_15 AC 303 ms
163,712 KB
testcase_16 AC 296 ms
163,712 KB
testcase_17 AC 298 ms
163,712 KB
testcase_18 AC 303 ms
163,712 KB
testcase_19 AC 286 ms
163,712 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