結果
問題 | No.1657 Sum is Prime (Easy Version) |
ユーザー |
![]() |
提出日時 | 2021-08-27 21:51:50 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 318 ms / 2,000 ms |
コード長 | 2,581 bytes |
コンパイル時間 | 4,586 ms |
コンパイル使用メモリ | 275,892 KB |
最終ジャッジ日時 | 2025-01-24 02:54:25 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#pragma region Macros // #pragma GCC target("avx2") #pragma GCC optimize("O3") #include <bits/stdc++.h> #include <atcoder/all> #define rep(i, n) for (long long i = 0; i < (n); i++) #define rrep(i, n) for (long long i = (n - 1); i >= 0; i--) #define ALL(v) v.begin(), v.end() #define endl "\n" #define fi first #define se second #define popcount(bit) __builtin_popcount(bit) #define popcountll(bit) __builtin_popcountll(bit) #define pb push_back #define eb emplace_back using namespace std; using namespace atcoder; using P = pair<int, int>; using PL = pair<long long, long long>; using Graph = vector<vector<int>>; typedef long long ll; template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } const int dx[4] = {1, 0, -1, 0}; const int dy[4] = {0, 1, 0, -1}; const int fx[8] = {0, 1, 1, 1, 0, -1, -1, -1}; const int fy[8] = {1, 1, 0, -1, -1, -1, 0, 1}; template <typename T> const auto INF = numeric_limits<T>::max()/2; struct Sieve_prime{ vector<int> f, primes; int n; Sieve_prime(int n = 1): n(n), f(n +1){ f[0] = f[1] = -1; for (ll i = 2; i <= n; i++){ if(f[i])continue; primes.emplace_back(i); f[i] = i; for (ll j = i * i; j <= n; j += i){ //この条件をつけると最小素因数が記録される //つけないと最大素因数が記録される。 if(!f[j]) f[j] = i; } } } bool isPrime(int x) { return f[x] == x; } vector<int> factorlist(int x){ vector<int> res; while (x != 1){ res.emplace_back(f[x]); x /= f[x]; } return res; } vector<P> factor(int x) { vector<int> fl = factorlist(x); if (fl.size() == 0) return {}; vector<P> res(1, P(fl[0], 0)); for (int p : fl) { if (res.back().first == p) { res.back().second++; } else { res.emplace_back(p, 1); } } return res; } vector<int> divisors(int x){ assert(x > 0); vector<int> res = {1}; for(auto [p, e] : factor(x)){ int n = res.size() * e; for(int i=0; i<n; ++i){ res.push_back(res[i] * p); } } sort(ALL(res)); return res; } }; Sieve_prime Sieve(20000000 + 2); int main() { ll l, r; cin >> l >> r; deque<ll> q; ll ans = 0; ll sum = 0; for (ll i = l; i <= r * 2; i++){ if(Sieve.isPrime(i)){ if(l <= i && i <= r)ans++; if(i == 2)continue; ll c = i / 2; if(l <= c && c + 1 <= r)ans++; } } cout << ans << endl; }