#283654 (C++14) No.732 3PrimeCounting

• 問題: No.732 3PrimeCounting
• ユーザー: memmemmi
• 提出日時: 2018-09-08 01:10:05
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 109 ms / 3,000 ms
• コード長: 1,413 bytes
• コンパイル時間: 872 ms
• コンパイル使用メモリ: 109,664 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-05-07 14:03:23
• 合計ジャッジ時間: 3,940 ms
• ジャッジサーバーID （参考情報）: judge4 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 5,248 KB
testcase_01	AC	2 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	1 ms; 5,376 KB
testcase_04	AC	2 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	AC	2 ms; 5,376 KB
testcase_08	AC	2 ms; 5,376 KB
testcase_09	AC	1 ms; 5,376 KB
testcase_10	AC	2 ms; 5,376 KB
testcase_11	AC	2 ms; 5,376 KB
testcase_12	AC	2 ms; 5,376 KB
testcase_13	AC	1 ms; 5,376 KB
testcase_14	AC	1 ms; 5,376 KB
testcase_15	AC	1 ms; 5,376 KB
testcase_16	AC	2 ms; 5,376 KB
testcase_17	AC	1 ms; 5,376 KB
testcase_18	AC	1 ms; 5,376 KB
testcase_19	AC	1 ms; 5,376 KB
testcase_20	AC	2 ms; 5,376 KB
testcase_21	AC	2 ms; 5,376 KB
testcase_22	AC	2 ms; 5,376 KB
testcase_23	AC	1 ms; 5,376 KB
testcase_24	AC	1 ms; 5,376 KB
testcase_25	AC	2 ms; 5,376 KB
testcase_26	AC	2 ms; 5,376 KB
testcase_27	AC	2 ms; 5,376 KB
testcase_28	AC	1 ms; 5,376 KB
testcase_29	AC	2 ms; 5,376 KB
testcase_30	AC	2 ms; 5,376 KB
testcase_31	AC	2 ms; 5,376 KB
testcase_32	AC	2 ms; 5,376 KB
testcase_33	AC	2 ms; 5,376 KB
testcase_34	AC	2 ms; 5,376 KB
testcase_35	AC	2 ms; 5,376 KB
testcase_36	AC	2 ms; 5,376 KB
testcase_37	AC	2 ms; 5,376 KB
testcase_38	AC	1 ms; 5,376 KB
testcase_39	AC	2 ms; 5,376 KB
testcase_40	AC	2 ms; 5,376 KB
testcase_41	AC	2 ms; 5,376 KB
testcase_42	AC	2 ms; 5,376 KB
testcase_43	AC	2 ms; 5,376 KB
testcase_44	AC	1 ms; 5,376 KB
testcase_45	AC	1 ms; 5,376 KB
testcase_46	AC	2 ms; 5,376 KB
testcase_47	AC	1 ms; 5,376 KB
testcase_48	AC	2 ms; 5,376 KB
testcase_49	AC	2 ms; 5,376 KB
testcase_50	AC	1 ms; 5,376 KB
testcase_51	AC	2 ms; 5,376 KB
testcase_52	AC	1 ms; 5,376 KB
testcase_53	AC	3 ms; 5,376 KB
testcase_54	AC	34 ms; 5,376 KB
testcase_55	AC	36 ms; 5,376 KB
testcase_56	AC	35 ms; 5,376 KB
testcase_57	AC	9 ms; 5,376 KB
testcase_58	AC	8 ms; 5,376 KB
testcase_59	AC	4 ms; 5,376 KB
testcase_60	AC	11 ms; 5,376 KB
testcase_61	AC	12 ms; 5,376 KB
testcase_62	AC	46 ms; 5,376 KB
testcase_63	AC	20 ms; 5,376 KB
testcase_64	AC	15 ms; 5,376 KB
testcase_65	AC	14 ms; 5,376 KB
testcase_66	AC	2 ms; 5,376 KB
testcase_67	AC	2 ms; 5,376 KB
testcase_68	AC	44 ms; 5,376 KB
testcase_69	AC	41 ms; 5,376 KB
testcase_70	AC	36 ms; 5,376 KB
testcase_71	AC	36 ms; 5,376 KB
testcase_72	AC	19 ms; 5,376 KB
testcase_73	AC	79 ms; 5,376 KB
testcase_74	AC	78 ms; 5,376 KB
testcase_75	AC	2 ms; 5,376 KB
testcase_76	AC	41 ms; 5,376 KB
testcase_77	AC	8 ms; 5,376 KB
testcase_78	AC	58 ms; 5,376 KB
testcase_79	AC	40 ms; 5,376 KB
testcase_80	AC	53 ms; 5,376 KB
testcase_81	AC	38 ms; 5,376 KB
testcase_82	AC	2 ms; 5,376 KB
testcase_83	AC	8 ms; 5,376 KB
testcase_84	AC	9 ms; 5,376 KB
testcase_85	AC	32 ms; 5,376 KB
testcase_86	AC	54 ms; 5,376 KB
testcase_87	AC	107 ms; 5,376 KB
testcase_88	AC	109 ms; 5,376 KB

ソースコード

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cfloat>
#include <map>
#include <utility>
#include <set>
#include <iostream>
#include <memory>
#include <string>
#include <vector>
#include <algorithm>
#include <functional>
#include <sstream>
#include <complex>
#include <stack>
#include <queue>
#include <iomanip>
#include <array>
#include <numeric>
#include <regex>
#include <bitset>
#include <deque>

using namespace std;
typedef long long ll;
typedef pair<int, int> p_ii;
#define REP(i,m,n) for(int i=(int)m ; i < (int) n ; ++i )
#define rep(i,n) REP(i,0,n)

const int INF = 1e9;
const double PI = acos(-1.0);
const ll MOD = 1e9 + 7;

bool is_prime[300100];
int prime[300100];
int cnt[300100];
int ex_max = -1;

void sieve(int n) {
	int p = 0;
	rep(i, n + 1)is_prime[i] = true;
	is_prime[0] = is_prime[1] = false;
	REP(i, 2, n + 1) {
		if (is_prime[i]) {
			prime[p++] = i;
			for (int j = 2 * i; j <= n; j += i) is_prime[j] = false;
		}
	}
}

int main() {
	int n;
	cin >> n;

	ll res = 0;
	sieve(3 * n+10);
	cnt[prime[0] + prime[1]]++;
	ex_max = max(ex_max, prime[0] + prime[1]);

	for (int i = 2; prime[i] <= n; i++) {
		for (int j = i + 1; prime[j] - prime[i] <= ex_max; j++) {
			res += cnt[prime[j] - prime[i]];
		}
		for (int j = 0; j < i; j++) {
			cnt[prime[j] + prime[i]]++;
		}
		ex_max = max(ex_max, prime[i] + prime[i - 1]);
	}

	cout << res << endl;

	return 0;
}