結果

問題 No.458 異なる素数の和
ユーザー bombrary_skbombrary_sk
提出日時 2019-05-19 16:40:08
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 520 ms / 2,000 ms
コード長 1,809 bytes
コンパイル時間 1,022 ms
コンパイル使用メモリ 61,528 KB
実行使用メモリ 495,880 KB
最終ジャッジ日時 2023-10-17 07:55:41
合計ジャッジ時間 18,989 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 514 ms
495,880 KB
testcase_01 AC 516 ms
495,880 KB
testcase_02 AC 514 ms
495,880 KB
testcase_03 AC 514 ms
495,880 KB
testcase_04 AC 512 ms
495,880 KB
testcase_05 AC 513 ms
495,880 KB
testcase_06 AC 516 ms
495,880 KB
testcase_07 AC 508 ms
495,880 KB
testcase_08 AC 520 ms
495,880 KB
testcase_09 AC 515 ms
495,880 KB
testcase_10 AC 513 ms
495,880 KB
testcase_11 AC 510 ms
495,880 KB
testcase_12 AC 512 ms
495,880 KB
testcase_13 AC 514 ms
495,880 KB
testcase_14 AC 512 ms
495,880 KB
testcase_15 AC 513 ms
495,880 KB
testcase_16 AC 511 ms
495,880 KB
testcase_17 AC 513 ms
495,880 KB
testcase_18 AC 515 ms
495,880 KB
testcase_19 AC 514 ms
495,880 KB
testcase_20 AC 515 ms
495,880 KB
testcase_21 AC 512 ms
495,880 KB
testcase_22 AC 512 ms
495,880 KB
testcase_23 AC 513 ms
495,880 KB
testcase_24 AC 513 ms
495,880 KB
testcase_25 AC 513 ms
495,880 KB
testcase_26 AC 511 ms
495,880 KB
testcase_27 AC 511 ms
495,880 KB
testcase_28 AC 515 ms
495,880 KB
testcase_29 AC 506 ms
495,880 KB
testcase_30 AC 514 ms
495,880 KB
権限があれば一括ダウンロードができます

ソースコード

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 dp1[3000][21000];
void solve1(int N)
{
  for (int i = 0; i < 3000; i++) {
    for (int j = 0; j < 21000; j++) {
      dp1[i][j] = -1;
    }
  }
  dp1[0][0] = 0;
  for (int i = 0; i < prime.size(); i++) {
    for (int j = 0; j <= N; j++) {
      if (dp1[i][j] < 0) continue;
      dp1[i + 1][j] = max(dp1[i + 1][j], dp1[i][j]);
      if (j + prime[i] > N) continue;
      dp1[i + 1][j + prime[i]] = max(dp1[i + 1][j + prime[i]], dp1[i][j] + 1);
    }
  }
}

int dp2[3000][21000];
void solve2(int N)
{
  for (int i = 0; i < 3000; i++) {
    for (int j = 0; j < 21000; j++) {
      dp2[i][j] = -1;
    }
  }
  dp2[0][0] = 0;
  for (int i = 0; i < prime.size(); i++) {
    for (int j = 0; j <= N; j++) {
      if (dp2[i][j] < 0) continue;
      dp2[i + 1][j] = max(dp2[i + 1][j], dp2[i][j]);
      dp2[i + 1][j + prime[i]] = max(dp2[i + 1][j + prime[i]], dp2[i][j] + 1);
    }
  }
}

int main()
{
  sieve(20000);


  solve1(20000);
  solve2(20000);

  int N; cin >> N;
  cout << dp1[prime.size()][N] << endl;

  //for (int i = 10000; i <= 20000; i++) {
  //  int a1 = dp1[prime.size()][i];
  //  int a2 = dp2[prime.size()][i];
  //  cout << i << ' ' << a1 << ' ' << a2 << endl;
  //  if (a1 != a2) {
  //    return 0;
  //  }
  //}

  return 0;
}
0