結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー buey_t
提出日時 2023-04-29 14:44:41
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 300 ms / 2,000 ms
コード長 1,642 bytes
コンパイル時間 270 ms
コンパイル使用メモリ 82,376 KB
実行使用メモリ 163,832 KB
最終ジャッジ日時 2024-11-18 10:58:11
合計ジャッジ時間 6,662 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

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