結果

問題 No.458 異なる素数の和
コンテスト
ユーザー onakaT_Titai
提出日時 2018-10-10 22:58:02
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 83 ms / 2,000 ms
コード長 1,798 bytes
コンパイル時間 816 ms
コンパイル使用メモリ 98,860 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-12 17:24:44
合計ジャッジ時間 2,649 ms
ジャッジサーバーID
（参考情報） 		judge / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <bitset>
#include <fstream>
#include <string>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <vector>
#include <algorithm>
#include <numeric>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <deque>
#include <functional>
#include <cctype>
#include <list>
#include <limits>
//#include <boost/multiprecision/cpp_int.hpp>

const double EPS = (1e-10);


using namespace std;
using Int = long long;
//using namespace boost::multiprecision;

const Int MOD = 1000000007;

Int mod_pow(Int x, Int n) {
    Int res = 1;
    while(n > 0) {
        if(n & 1) res = (res * x) % MOD; //ビット演算(最下位ビットが1のとき)
        x = (x * x) % MOD;
        n >>= 1; //右シフト(n = n >> 1)
    }
    return res;
}

// エラトステネスの篩
// [0,n) の範囲の i について、primes[i] != 0 ⇔ i は素数
vector<int> erat(int n) {
    vector<int> primes(n);
    for (int i = 2; i < n; ++i) primes[i] = i;
    for (int i = 2; i*i < n; ++i)
        if (primes[i])
            for (int j = i*i; j < n; j+=i) primes[j] = 0;
    // 素数のみをvectorに格納する場合は以下の行を追加 (0 である要素を全削除)
    // primes.erase(remove(primes.begin(), primes.end(), 0), primes.end());
    return primes;
}

int dp[20001];

int main(){
    cin.tie(0);

    int N; cin >> N;
    vector<int> e = erat(N+1);
    vector<int> p;
    for (int i = 0; i < N+1; i++){
        if (e[i] > 0) p.push_back(i);
    }
    
    for (auto i : p){
        for (int j = 20001; j >= 0; j--){
            if ((j == 0 || dp[j]) && j + i < 20001){
                dp[j+i] = max(dp[j+i], dp[j] + 1);
            } 
        }
    }

    if (dp[N] > 0) cout << dp[N] << endl;
    else cout << -1 << endl;
}
0