結果
問題 | No.843 Triple Primes |
ユーザー |
|
提出日時 | 2019-06-28 21:42:30 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 7 ms / 2,000 ms |
コード長 | 2,244 bytes |
コンパイル時間 | 1,827 ms |
コンパイル使用メモリ | 195,580 KB |
最終ジャッジ日時 | 2025-01-07 05:27:33 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 42 |
ソースコード
#include <bits/stdc++.h> #define LLI long long int #define FOR(v, a, b) for(LLI v = (a); v < (b); ++v) #define FORE(v, a, b) for(LLI v = (a); v <= (b); ++v) #define REP(v, n) FOR(v, 0, n) #define REPE(v, n) FORE(v, 0, n) #define REV(v, a, b) for(LLI v = (a); v >= (b); --v) #define ALL(x) (x).begin(), (x).end() #define RALL(x) (x).rbegin(), (x).rend() #define ITR(it, c) for(auto it = (c).begin(); it != (c).end(); ++it) #define RITR(it, c) for(auto it = (c).rbegin(); it != (c).rend(); ++it) #define EXIST(c,x) ((c).find(x) != (c).end()) #define fst first #define snd second #define popcount __builtin_popcount #define UNIQ(v) (v).erase(unique(ALL(v)), (v).end()) #define bit(i) (1LL<<(i)) #define sz(v) ((LLI)(v).size()) #ifdef DEBUG #include <misc/C++/Debug.cpp> #else #define dump(...) ((void)0) #endif #define gcd __gcd using namespace std; template <class T> constexpr T lcm(T m, T n){return m/gcd(m,n)*n;} template <typename I> void join(ostream &ost, I s, I t, string d=" "){for(auto i=s; i!=t; ++i){if(i!=s)ost<<d; ost<<*i;}ost<<endl;} template <typename T> istream& operator>>(istream &is, vector<T> &v){for(auto &a : v) is >> a; return is;} template <typename T, typename U> istream& operator>>(istream &is, pair<T,U> &p){is >> p.first >> p.second; return is;} template <typename T, typename U> bool chmin(T &a, const U &b){return (a>b ? a=b, true : false);} template <typename T, typename U> bool chmax(T &a, const U &b){return (a<b ? a=b, true : false);} template <typename T, size_t N, typename U> void fill_array(T (&a)[N], const U &v){fill((U*)a, (U*)(a+N), v);} vector<bool> prime_table(int n){ vector<bool> res(n+1, true); res[0] = res[1] = false; FOR(i,2,n) if(res[i]) for(int j=2*i; j<=n; j+=i) res[j] = false; return res; } int main(){ cin.tie(0); ios::sync_with_stdio(false); auto is_prime = prime_table(500000); int N; while(cin >> N){ LLI ans = 0; if(N == 1){ ans = 0; }else{ // r = 2, p = 2, q = 2 ans += 1; // p = 2, q is odd, r is odd // q + 2 = r^2 int temp = 0; FORE(r,3,N){ LLI q = r*r-2; if(q<=N and is_prime[r] and is_prime[q]) ++temp; } ans += temp*2; } cout << ans << endl; } return 0; }