summaryrefslogtreecommitdiff
path: root/test/GNUmakefile
diff options
context:
space:
mode:
Diffstat (limited to 'test/GNUmakefile')
0 files changed, 0 insertions, 0 deletions
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
struct SegTree {
	using T = ll; using U = ll;
	int n;
	static constexpr T E = 0; // Neutral element for combine
	static constexpr U UF = inf; // Unused value by updates
	vector<T> tree;
	int h;
	vector<U> lazy;
	vector<int> k; // size of segments (optional)

	SegTree(const vector<T>& a) : n(sz(a) + 1), tree(2 * n, E),
	//SegTree(int size, T def = E) : n(size + 1), tree(2 * n, def),
			h(__lg(2 * n)), lazy(n, UF), k(2 * n, 1) {
		copy(all(a), tree.begin() + n);
		for (int i = n - 1; i > 0; i--) {
			k[i] = 2 * k[2 * i];
			tree[i] = comb(tree[2 * i], tree[2 * i + 1]);
	}}

	T comb(T a, T b) {return a + b;} // Modify this + E

	void apply(int i, U val) { // And this + UF
		tree[i] = val * k[i];
		if (i < n) lazy[i] = val; // Don't forget this
	}

	void push_down(int i) {
		if (lazy[i] != UF) {
			apply(2 * i, lazy[i]);
			apply(2 * i + 1, lazy[i]);
			lazy[i] = UF;
	}}

	void push(int i) {
		for (int s = h; s > 0; s--) push_down(i >> s);
	}

	void build(int i) {
		while (i /= 2) {
			push_down(i);
			tree[i] = comb(tree[2 * i], tree[2 * i + 1]);
	}}

	void update(int l, int r, U val) {
		l += n, r += n;
		int l0 = l, r0 = r;
		push(l0), push(r0 - 1);
		for (; l < r; l /= 2, r /= 2) {
			if (l&1) apply(l++, val);
			if (r&1) apply(--r, val);
		}
		build(l0), build(r0 - 1);
	}

	T query(int l, int r) {
		l += n, r += n;
		push(l), push(r - 1);
		T resL = E, resR = E;
		for (; l < r; l /= 2, r /= 2) {
			if (l&1) resL = comb(resL, tree[l++]);
			if (r&1) resR = comb(tree[--r], resR);
		}
		return comb(resL, resR);
	}

	// Optional:
	int lower_bound(int l, int r, T x) {
		l += n, r += n;
		push(l), push(r - 1);
		int a[64] = {}, lp = 0, rp = 64;
		for (; l < r; l /= 2, r /= 2) {
			if (l&1) a[lp++] = l++;
			if (r&1) a[--rp] = --r;
		}
		for (int i : a) if (i != 0 && tree[i] >= x) { // Modify this
			while (i < n) {
				push_down(i);
				if (tree[2 * i] >= x) i = 2 * i; // And this
				else i = 2 * i + 1;
			}
			return i - n;
		}
		return -1;
	}
};