SegmentTree/old/update-rmin.cpp

• View this file on GitHub
• Last update: 2022-02-01 17:39:26+09:00

Code

#include <bits/stdc++.h>
using namespace std;

//一つの値を更新、区間の最小値/最大値

template<class T, const T&(*op)(const T&,const T&)>
struct segmentTree {
  segmentTree(int x){
    while(x > n) n *= 2;
    inf = -op(inf-1, -inf+1);
    d.resize(2*n - 1, inf);
  }
  void update(int k, T x){
    assert(0 <= k && k < n);
    k += n - 1;
    d[k] = x;
    while(k > 0){
      k = (k - 1) / 2;
      d[k] = op(d[k*2+1], d[k*2+2]);
    }
  }
  void set(int i, T x){
    assert(0 <= i && i < n);
    d[i + n - 1] = x;
  }
  T &operator[](int i){
    assert(0 <= i && i < n);
    return d[i + n - 1];
  }
  void build(){
    for(int k = n - 2; k >= 0; k--)
      d[k] = op(d[2*k+1], d[2*k+2]);
  }
  T query(int l, int r){
    assert(0 <= l && l <= r && r <= n);
    return query_sub(l, r, 0, 0, n);
  }
  int min_right(int l, int r, T x){
    assert(0 <= l && l <= r && r <= n);
    return min_right_sub(l, r, x, 0, 0, n);
  }
  int min_left(int l, int r, T x){
    assert(0 <= l && l <= r && r <= n);
    return min_left_sub(l, r, x, 0, 0, n);
  }
  void debug(int l = 0, int r = -1){
    if(r == -1) r = n;
    cerr << "[" << l << "," << r << "): ";
    for(int i = l; i < r; i++) cerr << d[i] << " ";
    cerr << "\n";
  }
  private:
  T inf = numeric_limits<T>::max();
  int n = 1;
  vector<T> d;
  T query_sub(int a, int b, int k, int l, int r){
    if(r <= a || b <= l) return inf;
    if(a <= l && r <= b) return d[k];
    T vlef = query_sub(a, b, k*2+1, l, (r+l)/2);
    T vrig = query_sub(a, b, k*2+2, (r+l)/2, r);
    return op(vlef, vrig);
  }
  //範囲外であればreturn a-1
  int min_right_sub(int a, int b, T x, int k, int l, int r){
    if(op(d[k]+op(1,-1),x)==x || r <= a || b <= l) return a - 1;
    if(k >= n - 1) return k - (n - 1);

    int vrig = min_right_sub(a, b, x, 2*k+2, (l+r)/2, r);
    if(vrig != a - 1) return vrig;
    return min_right_sub(a, b, x, 2*k+1, l, (l+r)/2);
  }
  //範囲外であればreturn b
  int min_left_sub(int a, int b, T x, int k, int l, int r){
    if(op(d[k]+op(1,-1),x)==x || r <= a || b <= l) return b;
    else if(k >= n - 1) return k - (n - 1);

    int vlef = min_left_sub(a, b, x, 2*k+1, l, (l+r)/2);
    if(vlef != b) return vlef;
    return min_left_sub(a, b, x, 2*k+2, (l+r)/2, r);
  }
};

#line 1 "SegmentTree/old/update-rmin.cpp"
#include <bits/stdc++.h>
using namespace std;

//一つの値を更新、区間の最小値/最大値

template<class T, const T&(*op)(const T&,const T&)>
struct segmentTree {
  segmentTree(int x){
    while(x > n) n *= 2;
    inf = -op(inf-1, -inf+1);
    d.resize(2*n - 1, inf);
  }
  void update(int k, T x){
    assert(0 <= k && k < n);
    k += n - 1;
    d[k] = x;
    while(k > 0){
      k = (k - 1) / 2;
      d[k] = op(d[k*2+1], d[k*2+2]);
    }
  }
  void set(int i, T x){
    assert(0 <= i && i < n);
    d[i + n - 1] = x;
  }
  T &operator[](int i){
    assert(0 <= i && i < n);
    return d[i + n - 1];
  }
  void build(){
    for(int k = n - 2; k >= 0; k--)
      d[k] = op(d[2*k+1], d[2*k+2]);
  }
  T query(int l, int r){
    assert(0 <= l && l <= r && r <= n);
    return query_sub(l, r, 0, 0, n);
  }
  int min_right(int l, int r, T x){
    assert(0 <= l && l <= r && r <= n);
    return min_right_sub(l, r, x, 0, 0, n);
  }
  int min_left(int l, int r, T x){
    assert(0 <= l && l <= r && r <= n);
    return min_left_sub(l, r, x, 0, 0, n);
  }
  void debug(int l = 0, int r = -1){
    if(r == -1) r = n;
    cerr << "[" << l << "," << r << "): ";
    for(int i = l; i < r; i++) cerr << d[i] << " ";
    cerr << "\n";
  }
  private:
  T inf = numeric_limits<T>::max();
  int n = 1;
  vector<T> d;
  T query_sub(int a, int b, int k, int l, int r){
    if(r <= a || b <= l) return inf;
    if(a <= l && r <= b) return d[k];
    T vlef = query_sub(a, b, k*2+1, l, (r+l)/2);
    T vrig = query_sub(a, b, k*2+2, (r+l)/2, r);
    return op(vlef, vrig);
  }
  //範囲外であればreturn a-1
  int min_right_sub(int a, int b, T x, int k, int l, int r){
    if(op(d[k]+op(1,-1),x)==x || r <= a || b <= l) return a - 1;
    if(k >= n - 1) return k - (n - 1);

    int vrig = min_right_sub(a, b, x, 2*k+2, (l+r)/2, r);
    if(vrig != a - 1) return vrig;
    return min_right_sub(a, b, x, 2*k+1, l, (l+r)/2);
  }
  //範囲外であればreturn b
  int min_left_sub(int a, int b, T x, int k, int l, int r){
    if(op(d[k]+op(1,-1),x)==x || r <= a || b <= l) return b;
    else if(k >= n - 1) return k - (n - 1);

    int vlef = min_left_sub(a, b, x, 2*k+1, l, (l+r)/2);
    if(vlef != b) return vlef;
    return min_left_sub(a, b, x, 2*k+2, (l+r)/2, r);
  }
};

Back to top page