datastructures/sparseTable.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

struct SparseTable {
	vector<vector<int>> st;
	vector<ll> *a;

	int better(int lidx, int ridx) {
		return a->at(lidx) <= a->at(ridx) ? lidx : ridx;
	}

	void init(vector<ll> *vec) {
		a = vec;
		st.assign(__lg(sz(*a)) + 1, vector<int>(sz(*a)));
		iota(all(st[0]), 0);
		for (int j = 0; (2 << j) <= sz(*a); j++) {
			for (int i = 0; i + (2 << j) <= sz(*a); i++) {
				st[j + 1][i] = better(st[j][i] , st[j][i + (1 << j)]);
	}}}

	int queryIdempotent(int l, int r) {
		int j = __lg(r - l); //31 - builtin_clz(r - l);
		return better(st[j][l] , st[j][r - (1 << j)]);
	}
}