-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathstarry_sky_tree.cpp
137 lines (119 loc) · 3 KB
/
starry_sky_tree.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
#define _CRT_SECURE_NO_WARNINGS
#include <cstdio>
#include <iostream>
#include <stack>
#include <queue>
#include <algorithm>
#include <functional>
#include <set>
#include <map>
#include <string>
#include <vector>
#include <cmath>
#include<sstream>
#include<list>
#include<iomanip>
#include <cstdlib>
#include <cstring>
#include <stack>
#include <bitset>
#include <cassert>
#include <stdlib.h>
#include <stdio.h>
using namespace std;
const int INF = 100000000;
const long long LINF = 3e18 + 7;
const int MAX_N = 200010;
const int MAX_W = 10002;
const int MAX_ARRAYK = 100000;
double PI = 3.14159265358979323846;
//using ll = long long;
// RMQ and RAQ (遅延評価セグメント木) Starry Sky Tree
// AOJ DSL_2_H https://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=DSL_2_H&lang=jp
// https://qiita.com/ageprocpp/items/9ea58ac181d31cfdfe02
// https://algo-logic.info/segment-tree/
int N;
long long dat[2 * MAX_N - 1];
long long lazy[2 * MAX_N - 1];
// initialization
void init(int n_) {
N = 1;
while (N < n_) N *= 2;
for (int i = 0; i < 2 * N - 1; i++) {
dat[i] = LINF;
lazy[i] = LINF;
}
}
void eval(int k) {
if (lazy[k] == LINF) {
return;
}
if (k < N - 1) {
lazy[k * 2 + 1] = (lazy[k * 2 + 1] == LINF) ? lazy[k] : lazy[k * 2 + 1] + lazy[k];
lazy[k * 2 + 2] = (lazy[k * 2 + 2] == LINF) ? lazy[k] : lazy[k * 2 + 2] + lazy[k];
}
dat[k] = (dat[k] == LINF) ? lazy[k] : dat[k] + lazy[k];
lazy[k] = LINF;
}
void update(int a, int b, long long x, int k, int l, int r) {
eval(k);
if (a <= l && r <= b) { // 完全に内側の時
lazy[k] = (lazy[k] == LINF) ? x : lazy[k] + x;
eval(k);
}
else if (a < r && l < b) { // 一部区間が被る時
update(a, b, x, k * 2 + 1, l, (l + r) / 2); // 左の子
update(a, b, x, k * 2 + 2, (l + r) / 2, r); // 右の子
dat[k] = min(dat[k * 2 + 1], dat[k * 2 + 2]);
}
}
// [a, b) の最小値を求める
// k は節点の番号、l, rはその節点が[l, r)に対応づいていることを表す。
// したがって、外からquery(a, b, 0, 0, n)をよぶ。
// (k(0)番目の節点、l(0) ~ r(n)の間を初めは指定)
long long query(int a, int b, int k, int l, int r) {
eval(k);
// [a, b)と[l, r)が交差しなければ INF
if (r <= a || b <= l) {
return LINF;
}
// [a, b) が[l, r)を完全に含んでいれば、この節点の値
if (a <= l && r <= b) {
return dat[k];
}
else {
long long v1 = query(a, b, k * 2 + 1, l, (l + r) / 2);
long long v2 = query(a, b, k * 2 + 2, (l + r) / 2, r);
return min(v1, v2);
}
}
int Q;
long long L[MAX_N], R[MAX_N];
long long X[MAX_N];
int T[MAX_N];
long long A[MAX_N];
int main() {
cin >> N >> Q;
for (int q = 0; q < Q; q++) {
cin >> T[q];
if (T[q] == 1) {
cin >> L[q] >> R[q];
}
else {
cin >> L[q] >> R[q] >> X[q];
}
}
init(N);
for (int i = 0; i < 2 * N - 1; i++) {
dat[i] = 0;
}
for (int q = 0; q < Q; q++) {
if (T[q] == 1) {
cout << query(L[q], R[q] + 1, 0, 0, N + 1) << endl;
}
else {
update(L[q], R[q] + 1, X[q], 0, 0, N + 1);
}
}
return 0;
}