結果

問題 No.458 異なる素数の和
ユーザー bombrary_sk
提出日時 2019-05-19 13:39:52
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 151 ms / 2,000 ms
コード長 1,090 bytes
コンパイル時間 577 ms
コンパイル使用メモリ 60,348 KB
実行使用メモリ 249,644 KB
最終ジャッジ日時 2024-09-17 06:27:56
合計ジャッジ時間 4,129 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>

using namespace std;

// エラトステネスの篩  O(nloglogn)
// create p (1<=p<=n and p is a prime nubmer)
#define MAX_N 100000
vector<long long> prime;
bool is_prime[MAX_N + 1];
long long sieve(long long n)
{
  long long p = 0;
  for (long long i = 0; i <= n; i++) is_prime[i] = true;
  is_prime[0] = is_prime[1] = false;
  for (long long i = 2; i <= n; i++) {
    if (!is_prime[i]) continue;
    prime.push_back(i);
    for (long long j = 2 * i; j <= n; j += i) {
      is_prime[j] = false;
    }
  }
  return p;
}

int dp[3000][21000];
int main()
{
  for (int i = 0; i < 3000; i++) {
    for (int j = 0; j < 21000; j++) {
      dp[i][j] = -1;
    }
  }
  int N;
  cin >> N;
  sieve(N);

  dp[0][0] = 0;
  for (int i = 0; i < prime.size(); i++) {
    for (int j = 0; j <= N; j++) {
      if (dp[i][j] < 0) continue;
      dp[i + 1][j] = max(dp[i + 1][j], dp[i][j]);
      if (j + prime[i] > N) continue;
      dp[i + 1][j + prime[i]] = max(dp[i + 1][j + prime[i]], dp[i][j] + 1);
    }
  }
  cout << dp[prime.size()][N] << endl;

  return 0;
}
0