結果
問題 |
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))