結果
問題 | No.144 エラトステネスのざる |
ユーザー | sue_charo |
提出日時 | 2017-04-17 13:35:08 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,628 bytes |
コンパイル時間 | 288 ms |
コンパイル使用メモリ | 86,936 KB |
実行使用メモリ | 87,852 KB |
最終ジャッジ日時 | 2023-09-26 10:37:38 |
合計ジャッジ時間 | 6,783 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 167 ms
84,528 KB |
testcase_01 | AC | 167 ms
79,760 KB |
testcase_02 | AC | 165 ms
79,692 KB |
testcase_03 | AC | 166 ms
79,812 KB |
testcase_04 | AC | 168 ms
79,616 KB |
testcase_05 | AC | 169 ms
79,768 KB |
testcase_06 | AC | 195 ms
82,636 KB |
testcase_07 | AC | 195 ms
82,724 KB |
testcase_08 | AC | 202 ms
82,724 KB |
testcase_09 | AC | 190 ms
82,476 KB |
testcase_10 | AC | 194 ms
82,528 KB |
testcase_11 | AC | 193 ms
82,756 KB |
testcase_12 | AC | 192 ms
82,468 KB |
testcase_13 | TLE | - |
testcase_14 | -- | - |
testcase_15 | -- | - |
testcase_16 | -- | - |
testcase_17 | -- | - |
testcase_18 | -- | - |
testcase_19 | -- | - |
ソースコード
# 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()