結果

問題 No.458 異なる素数の和
コンテスト
ユーザー doikimihiro
提出日時 2018-06-10 01:11:51
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 674 ms / 2,000 ms
コード長 1,553 bytes
コンパイル時間 928 ms
コンパイル使用メモリ 83,168 KB
実行使用メモリ 180,864 KB
最終ジャッジ日時 2024-06-30 13:07:19
合計ジャッジ時間 5,278 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <iomanip>
#include <iostream>
#include <vector>
using namespace std;
#define MAX_N 20001
int N;
// int dp[MAX_N + 10][MAX_N + 10];
vector<int> sieve(int n) {
    vector<int> prime;
    vector<bool> is_prime(n + 1, true);
    is_prime[0] = is_prime[1] = false;
    for (int i = 2; i <= n; ++i) {
        if (is_prime[i]) {
            prime.push_back(i);
            for (int j = 2 * i; j <= n; j += i) { is_prime[j] = false; }
        }
    }
    return prime;
}

int main() {
    cin >> N;
    // if (N == 1) {
    //     cout << -1 << endl;
    //     return 0;
    // }
    vector<int> prime   = sieve(N);
    int number_of_prime = prime.size();

    vector<vector<int>> dp(N + 1, vector<int>(number_of_prime + 1, -1));
    for (int i = 0; i <= 0; i++) {
        for (int j = 0; j <= number_of_prime; j++) { dp[i][j] = 0; }
    }

    for (int i = 0; i <= N; i++) {
        for (int j = 0; j < number_of_prime; j++) {
            dp[i][j + 1] = dp[i][j];
            if (i - prime[j] >= 0 && dp[i - prime[j]][j] != -1) {
                dp[i][j + 1] = max(dp[i][j + 1], dp[i - prime[j]][j] + 1);
            }
        }
    }
    // for (int i = 0; i <= N; i++) {
    //     for (int j = 0; j < number_of_prime; j++) {
    //         cout << setw(3) << dp[i][j];
    //     }
    //     cout << endl;
    // }

    cout << dp[N][number_of_prime] << endl;

    return 0;
}
// dp[i][j]は和がiになるもので、j番目の要素まで使った場合の最大の素数の数
// dp[i][j]=max(dp[i][j-1],dp[i-xj][j-1]+1);
0