結果
問題 |
No.732 3PrimeCounting
|
ユーザー |
|
提出日時 | 2021-03-25 19:30:44 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 840 ms / 3,000 ms |
コード長 | 1,007 bytes |
コンパイル時間 | 3,164 ms |
コンパイル使用メモリ | 197,040 KB |
最終ジャッジ日時 | 2025-01-19 21:41:49 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 89 |
ソースコード
#include <bits/stdc++.h> using namespace std; vector<bool> psieve(int n) { vector<bool> pr(n); for (int i = 2; i < n; i++) pr[i] = true; for (int i = 2; i * i < n; i++) { if (!pr[i]) continue; for (int j = i * i; j < n; j += i) pr[j] = false; } return pr; } int main() { ios_base::sync_with_stdio(0); cin.tie(0); int n; cin >> n; auto pr = psieve(3 * n + 1); vector<int> p, q; for (int i = 0; i <= n; i++) if (pr[i]) p.push_back(i); for (int i = 0; i <= 3 * n; i++) if (pr[i]) q.push_back(i); vector<int> suff(3 * n + 1); for (int x: q) { for (int y: p) { if (x - y < 0) break; suff[x - y]++; } } long long ret = 0; // a + b = x(prime) - c for (int i = 0; i < p.size(); i++) { // b // remove c = b for (int x: q) if (x - p[i] > 0) --suff[x - p[i]]; for (int j = 0; j < i; j++) { // a ret += suff[p[j] + p[i]]; } } cout << ret << endl; return 0; }