結果
| 問題 |
No.458 異なる素数の和
|
| ユーザー |
|
| 提出日時 | 2017-03-26 14:12:44 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
TLE
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 930 bytes |
| コンパイル時間 | 157 ms |
| コンパイル使用メモリ | 12,800 KB |
| 実行使用メモリ | 11,392 KB |
| 最終ジャッジ日時 | 2024-07-06 06:27:32 |
| 合計ジャッジ時間 | 13,670 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 26 TLE * 2 |
ソースコード
def primes2(limit):
''' returns a list of prime numbers upto limit.
source: Rossetta code: Sieve of Eratosthenes
http://rosettacode.org/wiki/Sieve_of_Eratosthenes#Odds-only_version_of_the_array_sieve_above
'''
if limit < 2: return []
if limit < 3: return [2]
lmtbf = (limit - 3) // 2
buf = [True] * (lmtbf + 1)
for i in range((int(limit ** 0.5) - 3) // 2 + 1):
if buf[i]:
p = i + i + 3
s = p * (i + 1) + i
buf[s::p] = [False] * ((lmtbf - s) // p + 1)
return [2] + [i + i + 3 for i, v in enumerate(buf) if v]
def solve(N):
primes = primes2(N)
dp = [0] * (N + 1)
dp[0] = 1
cumsum = 0
for p in primes:
cumsum += p
for i in range(min(cumsum, N), p - 1, -1):
tmp = dp[i - p]
if tmp and tmp + 1 > dp[i]:
dp[i] = tmp + 1
return dp[N] - 1
N = int(input())
print(solve(N))