Codeforces896C - Willem, Chtholly and Seniorious 珂朵莉树

        珂朵莉树是用来做区间推平跟一些奇怪的查询，比如本题区间幂次和，基于 set 来维护三元组 [ L，R，w ]，在数据随机的情况下，复杂的是能得到保证的。在 split 操作里要先 split [ R + 1 ] ，在 split [ L ]，反过来的话有可能会使得前一个迭代器失效。

        本题需要实现：区间加，区间更改，区间幂次和，区间第 K 小。

#include<bits/stdc++.h>
constexpr auto Inf = 0x3F3F3F3F;
#define Loc freopen("baka.int", "r", stdin)
typedef long long LL;
using namespace std;

/* 1000 1000 658073485 946088556 */
const int Mod = 1E9 + 7;
struct Node {
	int L, R; mutable LL w;
	Node(int L, int R = -1, LL w = 0) : L(L), R(R), w(w) {};
	bool operator < (const Node& Tmp) const {
		return L < Tmp.L;
	}
}; set<Node> Set; vector<pair<LL, int> > vec; LL Seed, vmax;

LL pow(LL w, LL N, LL M) {
	LL ans = 1; w %= M;
	while (N) {
		if (N & 1) ans = ans * w % M; w = w * w % M; N >>= 1;
	} return ans;
}

LL Get() {
	LL res = Seed; Seed = (Seed * 7 + 13) % Mod; return res;
}

set<Node>::iterator split(int pos) {
	auto iter = Set.upper_bound(Node(pos)); --iter;
	if (iter->L == pos)
		return iter;
	if (iter->R < pos)
		return Set.end();
	LL w = iter->w; int L = iter->L, R = iter->R;
	Set.erase(iter);
	Set.insert(Node(L, pos - 1, w));
	return Set.insert(Node(pos, R, w)).first;
}

void add(int _L, int _R, LL _w) {
	set<Node>::iterator R = split(_R + 1), L = split(_L);
	for (; L != R; ++L)
		L->w += _w;
}

void cov(int _L, int _R, LL _w) {
	set<Node>::iterator R = split(_R + 1), L = split(_L);
	Set.erase(L, R);
	Set.insert(Node(_L, _R, _w));
}

LL rak(int _L, int _R, int K) {
	vec.clear();
	set<Node>::iterator R = split(_R + 1), L = split(_L);
	for (; L != R; ++L)
		vec.push_back(make_pair(L->w, L->R - L->L + 1));
	sort(vec.begin(), vec.end());
	for (auto iter : vec) {
		K -= iter.second;
		if (K <= 0) return iter.first;
	} return -1;
}

LL psm(int _L, int _R, int Ex, int M) {
	LL ans = 0;
	set<Node>::iterator R = split(_R + 1), L = split(_L);
	for (; L != R; ++L)
		ans = (ans + (L->R - L->L + 1) * pow(L->w, Ex, M)) % M;
	return ans;
}

int main() {
	int N, M;
	cin >> N >> M >> Seed >> vmax;
	for (int pos = 1; pos <= N; pos++)
		Set.insert(Node(pos, pos, Get() % vmax + 1));
	for (int pos = 1, L, R, X, Y, o; pos <= M; pos++) {
		o = Get() % 4 + 1;
		L = Get() % N + 1;
		R = Get() % N + 1;
		if (L > R) swap(L, R);

		X = o == 3 ? Get() % (R - L + 1) + 1 : Get() % vmax + 1;
		if (o == 4) Y = Get() % vmax + 1;
		if (o == 1) add(L, R, X);
		else if (o == 2) cov(L, R, X);
		else if (o == 3) cout << rak(L, R, X) << endl;
		else cout << psm(L, R, X, Y) << endl;
	}
}