結果
| 問題 |
No.407 鴨等素数間隔列の数え上げ
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-02-23 05:47:04 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 342 ms / 1,000 ms |
| コード長 | 889 bytes |
| コンパイル時間 | 168 ms |
| コンパイル使用メモリ | 12,672 KB |
| 実行使用メモリ | 49,172 KB |
| 最終ジャッジ日時 | 2024-12-31 10:25:05 |
| 合計ジャッジ時間 | 3,635 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 5 |
| other | AC * 31 |
ソースコード
import sys
from math import sqrt
def solve():
N, L = map(int, input().split())
limit = L // (N - 1) + 1
primes = rwh_primes1(limit)
ans = 0
for p in primes:
if L - (N - 1) * p + 1 <= 0:
break
ans += L - (N - 1) * p + 1
print(ans)
def rwh_primes1(n):
# http://stackoverflow.com/questions/2068372/fastest-way-to-list-all-primes-below-n-in-python/3035188#3035188
""" Returns a list of primes < n """
sieve = [True] * (n//2)
for i in range(3,int(n**0.5)+1,2):
if sieve[i//2]:
sieve[i*i//2::i] = [False] * ((n-i*i-1)//(2*i)+1)
return [2] + [2*i+1 for i in range(1,n//2) if sieve[i]]
def debug(x, table):
for name, val in table.items():
if x is val:
print('DEBUG:{} -> {}'.format(name, val), file=sys.stderr)
return None
if __name__ == '__main__':
solve()