結果
問題 | No.106 素数が嫌い!2 |
ユーザー | maspy |
提出日時 | 2023-05-12 21:17:02 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 1,052 ms / 5,000 ms |
コード長 | 589 bytes |
コンパイル時間 | 107 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 56,112 KB |
最終ジャッジ日時 | 2024-05-06 10:44:28 |
合計ジャッジ時間 | 18,283 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1,044 ms
54,580 KB |
testcase_01 | AC | 1,052 ms
53,432 KB |
testcase_02 | AC | 1,026 ms
53,420 KB |
testcase_03 | AC | 1,022 ms
53,680 KB |
testcase_04 | AC | 1,017 ms
54,960 KB |
testcase_05 | AC | 1,012 ms
55,468 KB |
testcase_06 | AC | 1,034 ms
55,856 KB |
testcase_07 | AC | 1,016 ms
56,112 KB |
testcase_08 | AC | 1,039 ms
53,424 KB |
testcase_09 | AC | 996 ms
53,936 KB |
testcase_10 | AC | 1,028 ms
55,988 KB |
testcase_11 | AC | 1,030 ms
53,548 KB |
testcase_12 | AC | 1,021 ms
55,476 KB |
testcase_13 | AC | 1,020 ms
53,944 KB |
testcase_14 | AC | 1,038 ms
53,556 KB |
testcase_15 | AC | 999 ms
53,548 KB |
ソースコード
import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines import numpy as np N,K = map(int,read().split()) def make_prime(U): is_prime = np.zeros(U,np.bool_) is_prime[2] = 1 is_prime[3::2] = 1 M = int(U**.5)+1 for p in range(3,M,2): if is_prime[p]: is_prime[p*p::p+p] = 0 return is_prime, is_prime.nonzero()[0] U = 2 * 10 ** 6 + 10 _,primes = make_prime(U) cnt = np.zeros(U,np.int32) for p in primes: cnt[p::p] += 1 x = cnt[2:N+1] answer = np.count_nonzero(x >= K) print(answer)