結果

問題 No.843 Triple Primes
コンテスト
ユーザー Kichico
提出日時 2022-05-12 13:45:26
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 335 ms / 2,000 ms
コード長 2,890 bytes
コンパイル時間 2,057 ms
コンパイル使用メモリ 207,760 KB
最終ジャッジ日時 2025-01-29 06:19:44
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
using ll = int64_t;
using ld = long double;
using ull = unsigned long long;
#define ALL(x) x.begin(),x.end()
#define rep(iter,from,to) for(ll iter=from;iter<to;++iter)
#define fore(variable,container) for(auto& variable:container)
#define forc(variable,container) for(auto& variable:container) cout<<variable<<endl;
const ll MOD = 1e9 + 7;
const ll INF = 1e17;
const vector<ll> dx{ 1,0,-1,0 }, dy{ 0,1,0,-1 };
//#######################################################################
void op(vector<ll> vec) {
    ll size = (ll)vec.size();
    for (ll i = 0; i < size - 1; ++i) cout << vec[i] << " ";
    cout << vec.back() << endl;
}

void op(vector<vector<ll>> vec) {
    ll height = (ll)vec.size();
    ll width = (ll)vec[0].size();
    for (ll i = 0; i < height; ++i) {
        for (ll j = 0; j < width - 1; ++j) cout << vec[i][j] << " ";
        cout << vec[i].back() << endl;
    }
}

void twoText(bool identifier, string outTrue, string outFalse) {
    if (identifier) cout << outTrue << endl;
    else cout << outFalse << endl;
}

void twoText(bool identifier) {
    if (identifier) cout << "Yes" << endl;
    else cout << "No" << endl;
}

template <class T>
T vecsum(vector<T>& vec) {
    return accumulate(ALL(vec), (T)0);
}

template<class T, ll>
T vecsum(vector<T>& vec, ll K) {
    ll ret = 0;
    rep(i, 0, K) ret += vec[i];
    return ret;
}

template <class T>
struct grid {
    vector<vector<T>> field;
    grid(ll height, ll width) { field = vector<vector<T>>(height, vector<T>(width, (T)0)); }
    void input() { rep(i, 0, field.size()) rep(j, 0, field[i].size()) cin >> field[i][j]; }
};

//#########################################################################
bool isPrime(ll Num) {
    ll root = sqrt(Num);
    if (Num == 0 || Num == 1) return false;
    for (ll i = 2; i <= root; ++i) if (Num % i == 0) return false;
    return true;
}

vector<ll> enumeratePrime(ll maxx) {
    vector<bool> notp(maxx + 1, false);
    vector<ll> ret;
    rep(i, 2, maxx + 1) {
        if (notp[i]) continue;
        if (isPrime(i)) ret.emplace_back(i);
        ll init = 1;
        while (init * i <= maxx) {
            notp[init * i] = true;
            init++;
        }
    }
    return ret;
}

void solve() {
    ll N; cin >> N;
    auto primes = enumeratePrime(N);
    unordered_set<ll> primeSet;
    fore(x, primes) primeSet.emplace(x);
    ll ans = 0;
    fore(r, primes) {
        if (r * r > 2 * N) break;
        fore(p, primes) {
            //if (2 * p > r * r) break;
            ll q = r * r - p;
            if (q < 0) break;
            if (primeSet.find(q) != primeSet.end()) {
                ans++;
            }
        }
    }
    cout << ans << endl;
}


int main(void) {
    std::cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << std::fixed << std::setprecision(15);
    solve();
}
0