結果
問題 | No.843 Triple Primes |
ユーザー |
|
提出日時 | 2019-06-28 21:42:03 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 16 ms / 2,000 ms |
コード長 | 839 bytes |
コンパイル時間 | 1,615 ms |
コンパイル使用メモリ | 169,112 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-19 13:57:30 |
合計ジャッジ時間 | 2,924 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 42 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define REP(i,n) for (int i=0;i<(n);i++) #define REP2(i,m,n) for (int i=m;i<(n);i++) typedef long long ll; ll N; bool isP[1010101]; void solve() { cin >> N; isP[0] = isP[1] = false; REP2(i, 2, 1010101) isP[i] = true; vector<ll> P; ll pmax = 0; for (ll i = 2; i < 1010101; ++i) if (isP[i]) { if (i <= N) { P.push_back(i); pmax = max(pmax, i); } for (ll j = i + i; j < 1010101; j += i) { isP[j] = false; } } ll ans = 0; for (auto r: P) if (r * r <= 2 * pmax) { for (auto q: P) if (r * r - q >= 2 && r * r - q <= N && isP[r * r - q]) { ans += 1; } } cout << ans << endl; } int main() { cin.tie(0); ios::sync_with_stdio(false); solve(); }