test/library_checker/data_structure/ordered_set.test.cpp

• View this file on GitHub
• Last update: 2024-12-15 12:50:17+09:00
• Problem: https://judge.yosupo.jp/problem/ordered_set

Depends on

• Static Ordered Set (lib/data_structure/static_ordered_set.hpp)

Code

#define PROBLEM "https://judge.yosupo.jp/problem/ordered_set"

#include <iostream>
#include <vector>

#include "../../../lib/data_structure/static_ordered_set.hpp"

using namespace std;

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    int n, q; cin >> n >> q;
    StaticOrderedSet<long long> ps;
    vector<int> a(n);
    for(int i = 0; i < n; ++i){
        cin >> a[i];
        ps.add(a[i]);
    }
    using T = pair<int, int>;
    vector<T> query(q);
    for(int i = 0; i < q; ++i){
        cin >> query[i].first >> query[i].second;
        if(query[i].first != 2){
            ps.add(query[i].second);
        }
    }
    ps.build();

    for(int i = 0; i < n; ++i){
        ps.insert(a[i]);
    }
    for(int i = 0; i < q; i++){
        auto [type, x] = query[i];
        if(type == 0){
            if(!ps.exist(x)) ps.insert(x);
        } else if(type == 1){
            if(ps.exist(x)) ps.erase(x);
        } else if(type == 2){
            cout << ps.kth_small_element(x) << "\n";
        } else if(type == 3){
            cout << ps.count_less(x + 1) << "\n";
        } else if(type == 4){
            cout << ps.reverse_lower_bound(x) << "\n";
        } else if(type == 5){
            cout << ps.lower_bound(x) << "\n";
        }
    }
}

#line 1 "test/library_checker/data_structure/ordered_set.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/ordered_set"

#include <iostream>
#include <vector>

#line 2 "lib/data_structure/static_ordered_set.hpp"

/**
 * @brief Static Ordered Set
 */

#line 8 "lib/data_structure/static_ordered_set.hpp"
#include <algorithm>

template <typename T>
struct StaticOrderedSet{
    struct compress{
        std::vector<T> sorted, compressed;

        compress(){}

        void init(const std::vector<T> &vec){
            int n = vec.size();
            compressed.resize(n);
            for(T x : vec){
                sorted.emplace_back(x);
            }
            std::sort(sorted.begin(), sorted.end());
            sorted.erase(unique(sorted.begin(), sorted.end()), sorted.end());
            for(int i = 0; i < n; ++i){
                compressed[i] = std::lower_bound(sorted.begin(), sorted.end(), vec[i]) - sorted.begin();
            }
        }

        int get(const T &x) const{
            return std::lower_bound(sorted.begin(), sorted.end(), x) - sorted.begin();
        }

        T inv(const T &x) const{
            return sorted[x];
        }

        size_t size() const{
            return sorted.size();
        }

        std::vector<T> getCompressed() const{
            return compressed;
        }
    };

    struct BinaryIndexedTree{
        int N;
        std::vector<T> BIT;
        BinaryIndexedTree() {}

        void init(int size){
            N = size;
            BIT.assign(N + 1, 0);
        }

        T get(int i){
            return sum(i + 1) - sum(i);
        }

        void add(int i, T x){
            i++;
            while(i <= N){
                BIT[i] += x;
                i += i & -i;
            }
        }

        T sum(int i) const {
            T ans = 0;
            while(i > 0){
                ans += BIT[i];
                i -= i & -i;
            }
            return ans;
        }

        T sum(int L, int R) const {
            return sum(R) - sum(L);
        }

        int lower_bound(T x) const {
            if(x <= 0){
                return 0;
            } else{
                int v = 0, r = 1;
                while(r < N) r = r << 1;
                for(int len = r; len > 0; len = len >> 1){
                    if(v + len < N && BIT[v + len] < x){
                        x -= BIT[v + len];
                        v += len;
                    }
                }
                return v;
            }
        }

        int upper_bound(T x) const {
            if(x < 0){
                return 0;
            } else{
                int v = 0, r = 1;
                while(r <= N) r = r << 1;
                for(int len = r; len > 0; len = len >> 1){
                    if(v + len <= N && BIT[v + len] <= x){
                        x -= BIT[v + len];
                        v += len;
                    }
                }
                return v;
            }
        }

        T operator [](int i) const {
            return sum(i, i + 1);
        }
    };

    std::vector<T> a;
    compress comp;
    T siz, sum;
    BinaryIndexedTree cnt, val;

    StaticOrderedSet(){ }

    void add(T x){
        a.push_back(x);
    }

    void build(){
        comp.init(a);
        siz = T(0);
        sum = T(0);
        cnt.init(comp.size());
        val.init(comp.size());
    }

    T size(){
        return siz;
    }

    T all_sum(){
        return sum;
    }

    void insert(T x, T count = 1){
        int idx = comp.get(x);
        cnt.add(idx, count);
        val.add(idx, count * x);
        siz += count;
        sum += count * x;
    }

    void erase(T x, T count = 1){
        T idx = comp.get(x);
        if(cnt.get(idx) < count){
            count = cnt.get(idx);
        }
        cnt.add(idx, -count);
        val.add(idx, -count * x);
        siz -= count;
        sum -= count * x;
    }

    // 1-indexed
    T kth_small_element(T k){
        if(k <= 0 || size() < k){
            return -1;
        }
        T idx = cnt.lower_bound(k);
        return comp.inv(idx);
    }

    T kth_large_element(T k){
        if(k <= 0 || size() < k){
            return -1;
        }
        T rev_k = size() - k + 1;
        return kth_small_element(rev_k);
    }

    // 1-indexed
    T kth_small_sum(T k){
        if(size() < k){
            return val.sum((int) comp.size());
        }
        T idx = cnt.lower_bound(k);
        T sum = val.sum(idx);
        sum += comp.inv(idx) * (k - cnt.sum(idx));
        return sum;
    }

    T kth_large_sum(T k){
        if(size() < k){
            return val.sum((int) comp.size());
        }
        T rev_k = size() - k;
        return val.sum((int) comp.size()) - kth_small_sum(rev_k);
    }

    // x 未満
    T count_less(T x){
        return cnt.sum(comp.get(x));
    }

    // x 以上
    T count_moreq(T x){
        return size() - count_less(x);
    }

    T val_less(T x){
        return val.sum(comp.get(x));
    }

    T val_moreq(T x){
        return all_sum() - val_less(x);
    }

    T lower_bound(T x){
        return kth_small_element(count_less(x) + 1);
    }

    T upper_bound(T x){
        return kth_small_element(count_less(x + 1) + 1);
    }

    // x 以下で最大
    T reverse_lower_bound(T x){
        return kth_small_element(count_less(x + 1));
    }

    bool exist(T x){
        return count_less(x) != count_less(x + 1);
    }
};
#line 7 "test/library_checker/data_structure/ordered_set.test.cpp"

using namespace std;

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    int n, q; cin >> n >> q;
    StaticOrderedSet<long long> ps;
    vector<int> a(n);
    for(int i = 0; i < n; ++i){
        cin >> a[i];
        ps.add(a[i]);
    }
    using T = pair<int, int>;
    vector<T> query(q);
    for(int i = 0; i < q; ++i){
        cin >> query[i].first >> query[i].second;
        if(query[i].first != 2){
            ps.add(query[i].second);
        }
    }
    ps.build();

    for(int i = 0; i < n; ++i){
        ps.insert(a[i]);
    }
    for(int i = 0; i < q; i++){
        auto [type, x] = query[i];
        if(type == 0){
            if(!ps.exist(x)) ps.insert(x);
        } else if(type == 1){
            if(ps.exist(x)) ps.erase(x);
        } else if(type == 2){
            cout << ps.kth_small_element(x) << "\n";
        } else if(type == 3){
            cout << ps.count_less(x + 1) << "\n";
        } else if(type == 4){
            cout << ps.reverse_lower_bound(x) << "\n";
        } else if(type == 5){
            cout << ps.lower_bound(x) << "\n";
        }
    }
}

Back to top page