結果
問題 |
No.826 連絡網
|
ユーザー |
![]() |
提出日時 | 2020-02-06 12:51:46 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
RE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 519 bytes |
コンパイル時間 | 84 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 44,592 KB |
最終ジャッジ日時 | 2024-09-25 06:44:16 |
合計ジャッジ時間 | 17,329 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | RE * 2 |
other | RE * 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)