結果
問題 | No.458 異なる素数の和 |
ユーザー |
![]() |
提出日時 | 2020-09-04 12:48:13 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 240 ms / 2,000 ms |
コード長 | 1,287 bytes |
コンパイル時間 | 2,124 ms |
コンパイル使用メモリ | 199,124 KB |
最終ジャッジ日時 | 2025-01-14 04:26:10 |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include <bits/stdc++.h> using namespace std; class SieveOfEratosthenes { private: std::vector<int> sieve; public: SieveOfEratosthenes(int n) : sieve(n + 1) { std::iota(sieve.begin(), sieve.end(), 0); for (int i = 2; i * i < n + 1; i++) { if (sieve[i] != i) continue; for (int j = i * i; j < n + 1; j += i) { if (sieve[j] == j) sieve[j] = i; } } } bool is_prime(int x) const { return x != 0 && x != 1 && sieve[x] == x; } std::vector<int> factor(int x) const { std::vector<int> res; while (x > 1) { res.emplace_back(sieve[x]); x /= sieve[x]; } return res; } }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; SieveOfEratosthenes soe(N); vector<int> primes; for (int i = 0; i < N + 1; i++) { if (soe.is_prime(i)) primes.emplace_back(i); } vector<vector<int>> dp(primes.size() + 1, vector<int>(N + 1, 0)); dp[0][0] = 1; for (int i = 0; i < (int)primes.size(); i++) { for (int j = 0; j < N + 1; j++) { dp[i + 1][j] = max(dp[i + 1][j], dp[i][j]); if (j + primes[i] <= N && dp[i][j] > 0) dp[i + 1][j + primes[i]] = max(dp[i + 1][j + primes[i]], dp[i][j] + 1); } } cout << dp[primes.size()][N] - 1 << '\n'; return 0; }