結果
| 問題 | No.826 連絡網 | 
| コンテスト | |
| ユーザー |  maspy | 
| 提出日時 | 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)
            
            
            
        