結果

問題 No.843 Triple Primes
ユーザー lumc_
提出日時 2019-06-28 21:29:36
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 1,182 bytes
コンパイル時間 1,087 ms
コンパイル使用メモリ 104,284 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-19 13:48:42
合計ジャッジ時間 3,147 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff #

// includes {{{
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<tuple>
#include<cmath>
#include<random>
#include<cassert>
#include<bitset>
#include<cstdlib>
// #include<deque>
// #include<multiset>
// #include<cstring>
// #include<bits/stdc++.h>
// }}}
using namespace std;
using ll = long long;

// O(N log log N)
/// --- primes {{{ ///
#include <vector>
vector< int > primes(int n) {
  vector< int > res;
  for(int i = 2; i <= n; ++i) {
    bool isp = 1;
    for(int p : res) {
      if((ll) p * p > i) break;
      if(i % p == 0) {
        isp = 0;
        break;
      }
    }
    if(isp) res.emplace_back(i);
  }
  return res;
}
/// }}}--- ///

bool isprime[1123456];

int main() {
  std::ios::sync_with_stdio(false), std::cin.tie(0);
  int n;
  cin >> n;

  auto ps = primes(n);
  for(int p : ps) isprime[p] = 1;
  int ans = 0;

  for(int ri = 0; ri < ps.size() && ps[ri] * ps[ri] <= n * 2; ri++) {
    for(int pi = 0; pi < ps.size(); pi++) {
      int q = ps[ri] * ps[ri] - ps[pi];
      if(0 <= q && q <= n && isprime[q]) ans++;
    }
  }



  cout << ans << endl;
  return 0;
}
0