結果

問題 No.1234 典型RMQ
コンテスト
ユーザー intfans
提出日時 2024-05-06 22:58:57
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,834 bytes
コンパイル時間 2,120 ms
コンパイル使用メモリ 200,748 KB
最終ジャッジ日時 2025-02-21 11:28:20
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 11 WA * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

using ll = long long;

/*
link:https://yukicoder.me/problems/no/1234
长度为n的数目，q次操作，两种类型
1. a[l,r]的元素 +x
2. 求 a[l,r]区间最小值
*/

template<int B_size = -1>
struct SqrtDecomposition{
  using S = long long; // 待分块数据类型
  using T = long long; // 答案数据类型
  vector<S> A;
  vector<long long> lazy, mns;
  int B; // 分块大小
  int n; // 数组长度
  SqrtDecomposition(const int _n) : n(_n) {
    A.resize(_n);
    B = B_size == -1 ? sqrt(n) : B_size;
  }
  SqrtDecomposition(const vector<S>& v_) {
    A = v_;
    n = (int)(A.size());
    B = B_size == -1 ? sqrt(n) : B_size;
    init();
  }
  int get_start(int idx) const {
    return idx - idx % B;
  }
  int get_end(int bucket_start) const {
    return bucket_start + B < n ? bucket_start + B : n;
  }
  void set(int idx, S val){
    assert(idx >=0 && idx < n);
  }
  void init(){
    int m = (n + B - 1) / B;
    lazy.assign(m, 0);
    mns.resize(m);
    for (int i = 0, j = 0; i < n; i += B, j++) {
    	mns[j] = *min_element(A.begin() + i, A.begin() + get_end(i));
    }
  }
  T get_part(int l, int r) {
    T res = 1e18;
    int b = l / B;
    for(int i=l;i<r;++i){
    	if(A[i]+lazy[b]<res) res=A[i]+lazy[b];
    }
    return res;
  }
  T get_all(int l, int r) {
    return mns[l/B];
  }
  void build(int l, int r) { // [l, r)
  }
  void update_part(int l, int r, S x) {
  	int b = l / B;
  	for (int i = l; i < r; ++i) {
  		A[i] += x;
  	}
  	mns[b] = 1e18;
  	for (int i = get_start(l); i < get_end(b*B); ++i) {
  		mns[b] = min(mns[b], A[i] + lazy[b]);
  	}
    build(l,r);
  }
  void update_all(int l, int r, S x) {
  	lazy[l / B] += x;
  	mns[l / B] += x;
  }
  void prod(int l, int r, S val) {
    int x = l / B, y = r / B;
    if (x == y) {
      update_part(l, r, val);
    } else {
      update_part(get_start(l), x * B + B, val);
      for (int i = x + 1; i < y; ++i) {
        update_all(i * B, i * B + B, val);
      }
      update_part(y * B, r, val);
    }
  }
  T get(int l, int r) {
    int x = l / B, y = r / B;
    T res = 1e18;
    if (x == y) {
      res = min(res,get_part(l, r));
    } else {
      res += get_part(l, x * B + B);
      for (int i = x + 1; i < y; ++i) {
        res = min(res,get_all(i * B, i * B + B));
      }
      res = min(res, get_part(y * B, r));
    }
    return res;
  }
};

int main() {
    ios::sync_with_stdio(false); cin.tie(nullptr);
    
    int n, q;
    cin >> n;
    SqrtDecomposition<350> s(n);
    for (auto &x : s.A){
    	cin >> x;
    }
    s.init();
    cin >> q;
    for (int _i = 0; _i < q; ++_i) {
    	int k, l, r, c;
    	cin >> k >> l >> r >> c;
    	l--;
    	// continue;
    	if (k == 1) {
    		s.prod(l, r, c);
    	} else {
    		cout << s.get(l, r) << '\n';
    	}
    }

    return 0;
}
0