結果

問題 No.144 エラトステネスのざる
ユーザー mkawa2
提出日時 2020-01-03 00:44:50
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
TLE  
実行時間 -
コード長 1,009 bytes
コンパイル時間 234 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 124,212 KB
最終ジャッジ日時 2024-11-22 18:49:11
合計ジャッジ時間 29,001 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 7
権限があれば一括ダウンロードができます

ソースコード

diff # 

from itertools import *
from bisect import *
#from math import *
from collections import *
from heapq import *
from random import *
from decimal import *
import numpy as np
import sys

sys.setrecursionlimit(10 ** 6)
int1 = lambda x: int(x) - 1
p2D = lambda x: print(*x, sep="\n")
def II(): return int(sys.stdin.readline())
def MI(): return map(int, sys.stdin.readline().split())
def MI1(): return map(int1, sys.stdin.readline().split())
def MF(): return map(float, sys.stdin.readline().split())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LI1(): return list(map(int1, sys.stdin.readline().split()))
def LF(): return list(map(float, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
dij = [(1, 0), (0, 1), (-1, 0), (0, -1)]

def main():
    n,p=input().split()
    n=int(n)
    p=float(p)
    ev=[1]*(n+1)
    ans=0
    for x in range(2,n+1):
        ans+=ev[x]
        for y in range(x*2,n+1,x):
            ev[y]*=1-p
    print(ans)

main()
0