結果
問題 | No.3101 Range Eratosthenes Query |
ユーザー |
|
提出日時 | 2025-04-11 21:58:49 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 323 ms / 3,000 ms |
コード長 | 4,574 bytes |
コンパイル時間 | 10,953 ms |
コンパイル使用メモリ | 547,908 KB |
実行使用メモリ | 57,260 KB |
最終ジャッジ日時 | 2025-04-11 21:59:12 |
合計ジャッジ時間 | 19,138 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 24 |
ソースコード
#pragma GCC optimize("O3") #include "bits/stdc++.h" #include <stdio.h> #include "atcoder/all" #include <boost/multiprecision/cpp_int.hpp> using namespace boost::multiprecision; using namespace std; using namespace atcoder; using ll = long long; using ull = unsigned long long; using mint = modint998244353; using mint2 = modint1000000007; const double eps = 1e-9; #define OVERLOAD3(_1, _2, _3, call,...) call #define REP0(n) for (ll i = 0; i < ll(n); ++i) #define REP1(i, n) for (ll i = 0; i < ll(n); ++i) #define REP2(i, a, b) for (ll i = ll(a); i < ll(b); ++i) #define REP(...) OVERLOAD3(__VA_ARGS__, REP2, REP1, REP0)(__VA_ARGS__) #define RREP0(n) for (ll i = ll(n) - 1; i >= 0; --i) #define RREP1(i, n) for (ll i = ll(n) - 1; i >= 0; --i) #define RREP2(i, a, b) for (ll i = ll(a) - 1; i >= b; --i) #define RREP(...) OVERLOAD3(__VA_ARGS__, RREP2, RREP1, RREP0)(__VA_ARGS__) #define REPE(i, c) for (auto&& i : c) #define TCASE() ll _ttt; cin >> _ttt; while(_ttt--) template<class... T> void COMMON_IN(T&... a) { (cin >> ... >> a); } #define LI(...) ll __VA_ARGS__; COMMON_IN(__VA_ARGS__); using vl = vector<ll>; using vvl = vector<vl>; using vvvl = vector<vvl>; #define VLI(a, n) vector<ll> a; a.assign(n, 0); for(auto& x : a) cin >> x; #define DLI(a, n) deque<ll> a; a.assign(n, 0); for(auto& x : a) cin >> x; #define VSI(a, n) vector<string> a; a.assign(n, ""); for(auto& x : a) cin >> x; #define VVLI(a, n, m) vector<vector<ll>> a(n, vector<ll>(m)); for(auto& x : a) for(auto& y : x) cin >> y; using vd = vector<double>; using pl = pair<ll, ll>; using vpl = vector<pl>; struct uv { ll u; ll v; ll c = 1; }; #define VUVI(a, n) vector<uv> a; a.assign(n, {0, 0, 1}); for(auto& x : a) {cin >> x.u >> x.v; x.u--; x.v--;} #define VUVCI(a, n) vector<uv> a; a.assign(n, {0, 0, 0}); for(auto& x : a) {cin >> x.u >> x.v >> x.c; x.u--; x.v--;} vvl to_edge(const ll n, const vector<uv>& v) { vvl ret(n); for (auto& x : v) ret[x.u].push_back(x.v); return ret; } vvl to_edge_d(const ll n, const vector<uv>& v) { vvl ret(n); for (auto& x : v) { ret[x.u].push_back(x.v); ret[x.v].push_back(x.u); } return ret; } template <class T = long long> using pql = priority_queue<T>; template <class T = long long> using pqg = priority_queue<T, vector<T>, greater<T>>; using vm = vector<mint>; using vvm = vector<vm>; ostream& operator<<(ostream& os, const mint& v) { os << v.val(); return os; } void yn(bool f) { std::cout << (f ? "Yes" : "No") << "\n"; }; template<class T> void ov(const T& v) { for (auto it = v.begin(); it != v.end(); it++) { if (it != v.begin()) cout << " "; cout << *it; } }; template<> void ov(const vm& v) { for (auto it = v.begin(); it != v.end(); it++) { if (it != v.begin()) cout << " "; cout << it->val(); } }; const ll llhuge = 1LL << 61; //十分でかいが多少足しても溢れない数値 auto gmax(const auto& a, const auto& b) { return a > b ? a : b; }; auto gmin(const auto& a, const auto& b) { return a < b ? a : b; }; bool chmin(auto& a, const auto& b) { const auto olda = a; a = gmin(a, b); return olda != a; } bool chmax(auto& a, const auto& b) { const auto olda = a; a = gmax(a, b); return olda != a; } void sort(auto& v) { std::ranges::sort(v); }; void rev(auto& v) { std::ranges::reverse(v); }; template<class T> auto p_sum(const T& v) { T ret(v.size()); std::partial_sum(begin(v), end(v), ret.begin()); return ret; }; template<class T, class K = T::key_type, class V = T::mapped_type> [[nodiscard]] V myat(const T& mp, const K& key) { return mp.find(key) == mp.end() ? V() : mp.at(key); } long long ceil(const long long a, const long long b) { return (a + b - 1) / b; } //セグ木(sum) ll op(ll a, ll b) { return a + b; } ll e() { return 0; } using st = segtree<ll, op, e>; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); vl del(1000002, 0); REP(i, 1, ssize(del)) { for (ll mul = 2; i * mul < ssize(del); mul++) { chmax(del[i * mul], i); } } vpl mp; mp.reserve(ssize(del)); REP(i, 1, ssize(del)) { if (del[i] != 0) mp.push_back({ del[i], i }); } sort(mp.begin(), mp.end()); LI(q); vector<tuple<ll, ll, ll>> qli(q); REP(i, q) { LI(l, r); qli[i] = { l, r, i }; } sort(qli.begin(), qli.end(), [](const auto& a, const auto& b) { return get<0>(a) < get<0>(b); }); rev(qli); ll mpix = ssize(mp) - 1; vl seginit(ssize(del), 1); st seg(seginit); vl ans(q); REP(i, q) { auto [l, r, idx] = qli[i]; while (mpix >= 0 && mp[mpix].first >= l) { seg.set(mp[mpix].second, 0); mpix--; } ans[idx] = seg.prod(l, r + 1); } REP(i, q) { cout << ans[i] << "\n"; } }