結果

問題 No.144 エラトステネスのざる
ユーザー sue_charosue_charo
提出日時 2017-04-17 12:40:43
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 1,629 bytes
コンパイル時間 175 ms
コンパイル使用メモリ 10,936 KB
実行使用メモリ 14,500 KB
最終ジャッジ日時 2023-09-26 10:35:36
合計ジャッジ時間 4,911 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 30 ms
14,444 KB
testcase_01 AC 30 ms
9,964 KB
testcase_02 AC 31 ms
9,896 KB
testcase_03 AC 30 ms
10,068 KB
testcase_04 AC 30 ms
10,076 KB
testcase_05 AC 30 ms
9,928 KB
testcase_06 AC 35 ms
9,912 KB
testcase_07 AC 35 ms
9,896 KB
testcase_08 AC 34 ms
10,004 KB
testcase_09 AC 34 ms
10,072 KB
testcase_10 AC 34 ms
9,960 KB
testcase_11 AC 35 ms
10,088 KB
testcase_12 AC 35 ms
9,956 KB
testcase_13 TLE -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

# coding: utf-8
import array, bisect, collections, heapq, itertools, math, random, re, string, sys, time
sys.setrecursionlimit(10 ** 7)
inf = 10 ** 20
mod = 10 ** 9 + 7

def II(): return int(input())
def ILI(): return list(map(int, input().split()))
def IAI(LINE): return [ILI() for __ in range(LINE)]
def IDI(): return {key: value for key, value in ILI()}


def make_prime_list(num):
    if num < 2:
        return []

    # 0のものは素数じゃないとする
    prime_list = [i for i in range(num + 1)]
    prime_list[1] = 0   # 1は素数ではない
    num_sqrt = math.sqrt(num)

    for prime in prime_list:
        if prime == 0:
            continue
        if prime > num_sqrt:
            break

        for non_prime in range(2 * prime, num, prime):
            prime_list[non_prime] = 0

    return [prime for prime in prime_list if prime != 0]


def search_divisor_num_1(num):
    if num < 0:
        return None
    elif num == 1:
        return 1
    else:
        num_sqrt = math.floor(math.sqrt(num))
        prime_list = make_prime_list(num_sqrt)

        divisor_num = 1
        for prime in prime_list:
            count = 1
            while num % prime == 0:
                num //= prime
                count += 1
            divisor_num *= count

        if num != 1:
            divisor_num *= 2

        return divisor_num


def solve(N, p):
    ans = 0.0
    for n in range(2, N + 1):
        ans += (1.0 - p) ** (search_divisor_num_1(n) - 2)
    return ans


def main():
    N, p = input().split()
    N = int(N)
    p = float(p)
    print(solve(N, p))


if __name__ == "__main__":
    main()
0