結果

問題 No.3101 Range Eratosthenes Query
コンテスト
ユーザー syndrome
提出日時 2025-04-11 22:27:21
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,016 ms / 3,000 ms
コード長 3,066 bytes
コンパイル時間 2,496 ms
コンパイル使用メモリ 220,476 KB
実行使用メモリ 41,160 KB
最終ジャッジ日時 2025-04-11 22:28:15
合計ジャッジ時間 24,722 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 24
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// (⁠◕⁠ᴗ⁠◕⁠✿⁠)

// #pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
#define rep(i, n) for (ll i = 0; i < (n); i++)
#define srep(i, s, n) for (ll i = s; i < (n); i++)
#define len(x) ((int)(x).size())
#define all(x) (x).begin(), (x).end()
using namespace std;
template<typename T> using vc = vector<T>;
template<typename T> using vv = vc<vc<T>>;
template<typename T> using vvv = vv<vc<T>>;
using vi = vc<int>;using vvi = vv<int>; using vvvi = vv<vi>;
using ll = long long;using vl = vc<ll>;using vvl = vv<ll>; using vvvl = vv<vl>;
using ld = long double; using vld = vc<ld>; using vvld = vc<vld>; using vvvld = vc<vvld>;
using uint = unsigned int;
using ull = unsigned long long;
const ld pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
// const ll mod = 1000000007;
const ll mod = 998244353;
inline bool inside(ll y, ll x, ll H, ll W) {return 0 <= (y) and (y) < (H) and 0 <= (x) and (x) < (W); }

#define debug(var)  do{std::cout << #var << " : \n";view(var);}while(0)
template<typename T> void view(T e){cout << e << endl;}
template<typename T> void view(const vc<T>& v){for(const auto& e : v){ cout << e << " "; } cout << endl;}
template<typename T> void view(const vv<T>& vv){ for(const auto& v : vv){ view(v); } }

int op(int a, int b){return a + b;}
int e(){return 0;}
template<class T, T (*op)(T, T), T (*e)()> struct segtree{
    public : 
        segtree(int n) : _n(n){
            depth = 0;
            while ((1 << depth) < _n) depth++;
            size = 1 << depth;
            tree = vc<T>(2 * size, e());
        }

        void set(int p, T x){
            assert (0 <= p && p < _n);
            p += size;
            tree[p] = x;
            for (int i = 0; i < depth; i++){
                p >>= 1;
                tree[p] = op(tree[2 * p], tree[2 * p + 1]);
            }
            return;
        }

        T get(int p){
            assert(0 <= p && p < _n);
            return tree[p + size];
        }

        T prod(int l, int r){
            assert (0 <= l && l <= r && r <= _n);
            if (l == r) return e();
            T nl = e(), nr = e();
            l += size;
            r += size;
            while(l < r){
                if (l & 1) nl = op(nl, tree[l++]);
                if (r & 1) nr = op(nr, tree[--r]);
                l >>= 1;
                r >>= 1;
            }
            return op(nl, nr);
        }
        
        T all_prod(){return tree[1];}
    private : 
        int _n, depth, size;
        vc<T> tree;
};

vc<pair<int, int>> querys[1001001];
int main(){
    int Q; cin >> Q;
    rep(i, Q){
        int l, r; cin >> l >> r; r++;
        querys[l].push_back({r, i});
    }
    segtree<int, op, e> seg(1001001);
    rep(i, 1001001) seg.set(i, 1);
    vi ans(Q, -1);
    for (int i = 1001000; i >= 1; i--){
        for (int j = 2 * i; j < 1001001; j += i) seg.set(j, 0);
        for (auto [r, id] : querys[i]) ans[id] = seg.prod(i, r);
    }
    for (auto x : ans) cout << x << endl;
}
0