結果

問題 No.843 Triple Primes
ユーザー kibunakibuna
提出日時 2019-06-28 21:46:49
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 110 ms / 2,000 ms
コード長 1,729 bytes
コンパイル時間 1,664 ms
コンパイル使用メモリ 172,484 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-19 17:40:37
合計ジャッジ時間 5,060 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 110 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 2 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 3 ms
4,348 KB
testcase_07 AC 83 ms
4,348 KB
testcase_08 AC 92 ms
4,348 KB
testcase_09 AC 109 ms
4,348 KB
testcase_10 AC 87 ms
4,348 KB
testcase_11 AC 98 ms
4,348 KB
testcase_12 AC 106 ms
4,348 KB
testcase_13 AC 103 ms
4,348 KB
testcase_14 AC 103 ms
4,348 KB
testcase_15 AC 85 ms
4,348 KB
testcase_16 AC 86 ms
4,348 KB
testcase_17 AC 2 ms
4,348 KB
testcase_18 AC 2 ms
4,348 KB
testcase_19 AC 2 ms
4,348 KB
testcase_20 AC 20 ms
4,348 KB
testcase_21 AC 10 ms
4,348 KB
testcase_22 AC 51 ms
4,348 KB
testcase_23 AC 52 ms
4,348 KB
testcase_24 AC 21 ms
4,348 KB
testcase_25 AC 15 ms
4,348 KB
testcase_26 AC 108 ms
4,348 KB
testcase_27 AC 4 ms
4,348 KB
testcase_28 AC 104 ms
4,348 KB
testcase_29 AC 23 ms
4,348 KB
testcase_30 AC 107 ms
4,348 KB
testcase_31 AC 6 ms
4,348 KB
testcase_32 AC 3 ms
4,348 KB
testcase_33 AC 26 ms
4,348 KB
testcase_34 AC 44 ms
4,348 KB
testcase_35 AC 110 ms
4,348 KB
testcase_36 AC 13 ms
4,348 KB
testcase_37 AC 79 ms
4,348 KB
testcase_38 AC 61 ms
4,348 KB
testcase_39 AC 105 ms
4,348 KB
testcase_40 AC 2 ms
4,348 KB
testcase_41 AC 2 ms
4,348 KB
testcase_42 AC 90 ms
4,348 KB
testcase_43 AC 98 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using vi = vector<int>;
using vl = vector<ll>;
using vvi = vector<vi>;
using vvl = vector<vl>;
const int INF = 1 << 28;
const ll MOD = 1000000007;
template <class T>
bool chmax(T &a, const T &b) {
    return (a < b) ? (a = b, 1) : 0;
}
template <class T>
bool chmin(T &a, const T &b) {
    return (b < a) ? (a = b, 1) : 0;
}

struct Sieve {
    vector<bool> isPrime;
    Sieve(size_t max) : isPrime(max + 1, true) {
        isPrime[0] = false; // 0は素数ではない
        isPrime[1] = false; // 1は素数ではない

        for (size_t i = 2; i * i <= max; ++i)         // 0からsqrt(max)まで調べる
            if (isPrime[i])                           // iが素数ならば
                for (size_t j = 2; i * j <= max; ++j) // (max以下の)iの倍数は
                    isPrime[i * j] = false;           // 素数ではない
    }
    bool operator()(size_t n) { return isPrime[n]; }
};

int main() {
    int n;
    cin >> n;
    int m = int(sqrt(double(2 * n)));
    Sieve si(n);
    vi primesr;
    vi primespq;
    primesr.reserve(n);
    primespq.reserve(n);
    for (int i = 2; i <= m; ++i) {
        if (si(i))
            primesr.push_back(i);
    }
    for (int i = 2; i <= n; ++i) {
        if (si(i))
            primespq.push_back(i);
    }
    int ret = 0;
    for (auto &r : primesr) {
        for (auto &p : primespq) {
            if (p > r * r / 2)
                break;
            if (binary_search(primespq.begin(), primespq.end(), r * r - p))
                ret += 2;
        }
    }
    if (n > 1)
        ret -= 1;
    cout << ret << "\n";

    return 0;
}
0