結果
問題 |
No.826 連絡網
|
ユーザー |
![]() |
提出日時 | 2023-05-12 21:17:26 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 525 ms / 2,000 ms |
コード長 | 520 bytes |
コンパイル時間 | 177 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 44,588 KB |
最終ジャッジ日時 | 2024-11-28 16:59:55 |
合計ジャッジ時間 | 18,429 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 30 |
ソースコード
import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines import numpy as np N, P = map(int, read().split()) U = N + 1 is_prime = np.zeros(U, np.bool_) is_prime[2] = 1 is_prime[3::2] = 1 for p in range(3, U, 2): if p * p >= U: break if is_prime[p]: is_prime[p*p::p+p] = 0 M = N // 2 + 1 n_single = np.count_nonzero(is_prime[M:]) + 1 if P == 1 or (P >= M and is_prime[P]): answer = 1 else: answer = N - n_single print(answer)