結果
| 問題 |
No.458 異なる素数の和
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-03-26 14:11:43 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 116 ms / 2,000 ms |
| コード長 | 930 bytes |
| コンパイル時間 | 196 ms |
| コンパイル使用メモリ | 82,440 KB |
| 実行使用メモリ | 61,952 KB |
| 最終ジャッジ日時 | 2024-07-06 06:27:17 |
| 合計ジャッジ時間 | 2,516 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 28 |
ソースコード
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))