結果
問題 | No.843 Triple Primes |
ユーザー |
![]() |
提出日時 | 2019-06-28 21:46:52 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 128 ms / 2,000 ms |
コード長 | 1,707 bytes |
コンパイル時間 | 1,449 ms |
コンパイル使用メモリ | 166,552 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-19 13:59:59 |
合計ジャッジ時間 | 5,362 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 42 |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll = long long; template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; } #define FOR(i,a,b) for(ll i=(a);i<(b);++i) #define ALL(v) (v).begin(), (v).end() #define p(s) cout<<(s)<<endl #define p2(s, t) cout << (s) << " " << (t) << endl #define br() p("") #define pn(s) cout << (#s) << " " << (s) << endl #define p_yes() p("Yes") #define p_no() p("No") const ll mod = 1e9 + 7; const ll inf = 1e18; const int N_MAX = 1e6 + 10; bool is_prime[N_MAX]; // エラトステネスのふるい void Eratosthenes(){ FOR(i, 0, N_MAX){ is_prime[i] = true; } is_prime[0] = false; is_prime[1] = false; for(ll i=2; i*i<=N_MAX; i++){ if(is_prime[i]){ int p = i; int step_p = p; // それの倍数を消していく while(p < N_MAX){ p += step_p; if(p<N_MAX){ is_prime[p] = 0; } } } } } int main(){ cin.tie(0); ios::sync_with_stdio(false); // input ll N; cin >> N; Eratosthenes(); ll limit = min(N+1, (ll)sqrt(500000)+100); ll count = 0; FOR(r, 2, limit){ if(!is_prime[r]) continue; FOR(p, 2, N+1){ if(!is_prime[p]) continue; ll q = r*r - p; if(q<=1) continue; if(q>N) continue; if(is_prime[q]){ // cout << "pqr " << p << ' ' << q << ' ' << r << endl; count++; } } } p(count); return 0; }