[CF896C]Willem, Chtholly and Seniorious(珂朵莉树)

https://www.cnblogs.com/WAMonster/p/10181214.html

主要用于支持含有较难维护的区间操作与查询的问题，要求其中区间赋值操作（assign()）是纯随机的。

注意要先split(r+1)再split(l)，最好最后设一个点(n+1,n+1,0)

 #include<set>
 #include<cstdio>
 #include<algorithm>
 #include<iostream>
 #define rep(i,l,r) for (int i=(l); i<=(r); i++)
 typedef long long ll;
 using namespace std;

 const int N=;
 int n,m,x,y,mx; ll seed;
 struct P{ int l,r; mutable ll v; };
 bool operator <(const P &a,const P &b){ return a.l<b.l; }
 set<P>S;
 typedef set<P>::iterator sit;
 pair<ll,int>ve[N];

 ll ksm(ll a,ll b,ll mod){
     ll res=; a%=mod;
     for (; b; a=a*a%mod,b>>=)
         if (b & ) res=res*a%mod;
     return res;
 }

 sit split(int pos){
     sit it=S.lower_bound((P){pos,,});
     if (it!=S.end() && it->l==pos) return it;
     it--; int l=it->l,r=it->r; ll v=it->v;
     S.erase(it); S.insert((P){l,pos-,v});
     return S.insert((P){pos,r,v}).first;
 }

 void assign(int l,int r,ll v){
     sit it2=split(r+),it1=split(l);
     S.erase(it1,it2); S.insert((P){l,r,v});
 }

 void add(int l,int r,ll v){
     sit it2=split(r+),it1=split(l);
     for (sit it=it1; it!=it2; it++) it->v+=v;
 }

 ll kth(int l,int r,int k){
     sit it2=split(r+),it1=split(l); int tot=;
     for (sit it=it1; it!=it2; it++) ve[++tot]=pair<ll,int>(it->v,(it->r)-(it->l)+);
     sort(ve+,ve+tot+);
     rep(i,,tot){
         k-=ve[i].second;
         if (k<=) return ve[i].first;
     }
     return ;
 }

 ll que(int l,int r,int x,ll y){
     sit it2=split(r+),it1=split(l);
     ll res=;
     for (sit it=it1; it!=it2; it++) res=(res+((it->r)-(it->l)+)*ksm(it->v,x,y))%y;
     return res;
 }

 int rnd(){ int ret=(int)seed; seed=(seed*+)%; return ret; }

 int main(){
     freopen("CF896C.in","r",stdin);
     freopen("CF896C.out","w",stdout);
     cin>>n>>m>>seed>>mx;
     rep(i,,n) S.insert((P){i,i,rnd()%mx+});
     S.insert((P){n+,n+,});
     rep(i,,m){
         int op=rnd()%+,l=rnd()%n+,r=rnd()%n+;
         if (l>r) swap(l,r);
         if (op==) x=rnd()%(r-l+)+; else x=rnd()%mx+;
         if (op==) y=rnd()%mx+;
         if (op==) add(l,r,x);
         if (op==) assign(l,r,x);
         if (op==) cout<<kth(l,r,x)<<endl;
         if (op==) cout<<que(l,r,x,y)<<endl;
     }
     return ;
 }