結果
問題 | No.843 Triple Primes |
ユーザー |
![]() |
提出日時 | 2019-06-28 21:40:59 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 32 ms / 2,000 ms |
コード長 | 1,009 bytes |
コンパイル時間 | 1,855 ms |
コンパイル使用メモリ | 193,040 KB |
最終ジャッジ日時 | 2025-01-07 05:27:10 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 42 |
ソースコード
#include <bits/stdc++.h> #define rep(i, a, n) for(int i = a; i < n; i++) #define int long long using namespace std; typedef pair<int, int> P; const int mod = 1000000007; const int INF = 1e18; const int MAX_N = 500010; int prime[MAX_N]; bool is_prime[MAX_N + 1]; int sieve(int n){ int p = 0; for(int i = 0; i <= n; i++) is_prime[i] = true; is_prime[0] = is_prime[1] = false; for(int i = 2; i <= n; i++){ if(is_prime[i]){ prime[p++] = i; for(int j = 2*i; j <= n; j += i) is_prime[j] = false; } } return p; } signed main(){ cin.tie(nullptr); ios::sync_with_stdio(false); int n; cin >> n; int p = sieve(n + 1); int ans = 0; rep(i, 0, p){ int tmp = prime[i] * prime[i]; if(tmp > 1000000) break; rep(j, 0, p){ int tmp2 = tmp - prime[j]; if(tmp2 <= 0) break; if(tmp2 > n) continue; if(is_prime[tmp2]) ans++; } } cout << ans << endl; }