結果
問題 | No.106 素数が嫌い!2 |
ユーザー | H3PO4 |
提出日時 | 2020-04-26 14:58:29 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 1,395 ms / 5,000 ms |
コード長 | 468 bytes |
コンパイル時間 | 427 ms |
コンパイル使用メモリ | 12,416 KB |
実行使用メモリ | 74,792 KB |
最終ジャッジ日時 | 2024-11-17 07:07:52 |
合計ジャッジ時間 | 17,554 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 849 ms
59,560 KB |
testcase_01 | AC | 479 ms
43,980 KB |
testcase_02 | AC | 471 ms
43,940 KB |
testcase_03 | AC | 464 ms
44,204 KB |
testcase_04 | AC | 824 ms
59,560 KB |
testcase_05 | AC | 1,272 ms
74,668 KB |
testcase_06 | AC | 1,395 ms
74,792 KB |
testcase_07 | AC | 1,393 ms
74,660 KB |
testcase_08 | AC | 529 ms
44,584 KB |
testcase_09 | AC | 526 ms
44,720 KB |
testcase_10 | AC | 1,246 ms
74,792 KB |
testcase_11 | AC | 782 ms
57,508 KB |
testcase_12 | AC | 1,213 ms
73,768 KB |
testcase_13 | AC | 466 ms
44,076 KB |
testcase_14 | AC | 478 ms
44,076 KB |
testcase_15 | AC | 479 ms
44,324 KB |
ソースコード
import numpy as np N, K = map(int, input().split()) def primes(n): is_prime = [True] * (n + 1) is_prime[0] = False is_prime[1] = False for i in range(2, int(n ** 0.5) + 1): if not is_prime[i]: continue for j in range(i * 2, n + 1, i): is_prime[j] = False return (i for i in range(n + 1) if is_prime[i]) A = np.zeros(N + 1, dtype=int) for p in primes(N): A[p::p] += 1 print(np.count_nonzero(A >= K))