[ODT]CF896C

写一种数据结构，支持：

• $1$ $l$ $r$ $x$ ：将$[l,r]$ 区间所有数加上$x$
• $2$ $l$ $r$ $x$ ：将$[l,r]$ 区间所有数改成$x$
• $3$ $l$ $r$ $x$ ：输出将$[l,r]$ 区间从小到大排序后的第$x$ 个数是的多少(即区间第$x$ 小，数字大小相同算多次，保证 $1\le$ $x$ $\le$ $r-l+1$ )
• $4$ $l$ $r$ $x$ $y$ ：输出$[l,r]$ 区间每个数字的$x$ 次方的和模$y$ 的值(即(${\sum }_{i=l}^r{a}_i^x$ ) $\mathrm{mod}y$ )

其中数据在本机上根据程序生成（生成器是随机数据）， $n,m\le 10^5,x,y\le 10^9$ 。

---

珂朵莉树的元老题。首先声明珂朵莉树只有在随机数据下保证复杂度正确，大多数情况下只能用来骗分。然后用这一题来写一下ODT的基本流程。

首先珂朵莉树并不是树，只是借用了STL::set来解决区间问题，其中set中存储的结构体 $(l,r,v)$ 表示 $[l,r]$ 区间存储的值为 $v$ 。同时支持合并 $assign$ 和分割 $split$ 操作，分别表示将 $[l,r]$ 区间的值全部附为 $v$ ，将 $[l,r]$ 区间改为 $[l,x),[x,r]$ 两段，并返回后一段的指针。然后其他操作直接在STL上暴力 $++it$ 即可。

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

namespace QuickMul{
	ll Qmul(ll a,ll b,ll p){
		ll ret=1;
		while(b){
			if(b&1)
				ret=(ret*a)%p;
			b>>=1;
			a=((a%p)*(a%p))%p;
		}
		return ret;
	}
} 

namespace ChthollyTree{
	struct node_c{
		ll l,r;
		mutable ll v;
		node_c(ll il,ll ir,ll iv){
			l=il,r=ir,v=iv;
		}
		inline bool operator <(const node_c o)const{
			return l<o.l;
		}
	};
	struct rank{
		ll num,cnt;
		bool operator < (const rank a)const{
			return num<a.num;
		}
		rank(ll inum,ll icnt){
			num=inum;
			cnt=icnt;
		}
		rank(){
			num=0;
			cnt=0;
		}
	};
	#define iter set<node_c>::iterator
	set<node_c> odt;
	iter split(ll x){
		iter it=odt.lower_bound(node_c(x,0,0));
		if(it!=odt.end() && it->l == x){
			return it;
		}
		it--;
		if(it->r < x)
			return odt.end();
		ll l=it->l;
		ll r=it->r;
		ll v=it->v;
		odt.erase(it);
		odt.insert(node_c(l,x-1,v));
		return odt.insert(node_c(x,r,v)).first; 
	}
	void assign(ll l,ll r,ll v){
		iter itr=split(r+1),itl=split(l);
		odt.erase(itl,itr);
		odt.insert(node_c(l,r,v));
	}
	void add(int l,int r,int x){
		iter itr=split(r+1),itl=split(l);
		for(;itl!=itr;itl++){
			itl->v+=x;
		}
	}
	ll getkth(int l,int r,int k){
		iter itr=split(r+1),itl=split(l);
		vector<rank> v;
		for(;itl!=itr;itl++){
			v.push_back(rank(itl->v,(itl->r)-(itl->l)+1));
		}
		sort(v.begin(),v.end());
		for(int i=0;i<(int)v.size();++i){
			if(v[i].cnt>=k)
				return v[i].num;
			else
				k-=v[i].cnt;
		}
	}
	ll getpow(ll l,ll r,ll x,ll y){
		iter itr=split(r+1),itl=split(l);
		ll ret=0;
		for(;itl!=itr;++itl){
			ret=(ret+((((itl->r)-(itl->l)+1)%y)*(QuickMul::Qmul(itl->v,x,y))%y))%y;
		}
		return ret;
	}
}
using namespace ChthollyTree;
ll n,m,seed,vmax,a[100005];

ll rnd(){
	ll ret=seed;
	seed=(seed*7+13)%1000000007;
	return ret;
}

int main(){
	cin>>n>>m>>seed>>vmax;
	for(int i=1;i<=n;++i)
		a[i]=(rnd()%vmax)+1;
	for(int i=1;i<=n;++i){
		odt.insert(node_c(i,i,a[i]));
	}
	for(int i=1;i<=m;++i){
		ll op=(rnd()%4)+1;
		ll l=(rnd()%n)+1;
		ll r=(rnd()%n)+1;
		ll x,y;
		if(l>r)
			swap(l,r);
		if(op==3){
			x=(rnd()%(r-l+1))+1;
			printf("%lld\n",getkth(l,r,x));
		}
		else{
			x=(rnd()%vmax)+1;
		}
		if(op==4){
			y=(rnd()%vmax)+1;
			printf("%lld\n",getpow(l,r,x,y));
		}
		if(op==1)
			add(l,r,x);
		if(op==2)
			assign(l,r,x); 
	}
	
}