[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=100010;
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=1; a%=mod;
    for (; b; a=a*a%mod,b>>=1)
        if (b & 1) res=res*a%mod;
    return res;
}

sit split(int pos){
    sit it=S.lower_bound((P){pos,0,0});
    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-1,v});
    return S.insert((P){pos,r,v}).first;
}

void assign(int l,int r,ll v){
    sit it2=split(r+1),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+1),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+1),it1=split(l); int tot=0;
    for (sit it=it1; it!=it2; it++) ve[++tot]=pair<ll,int>(it->v,(it->r)-(it->l)+1);
    sort(ve+1,ve+tot+1);
    rep(i,1,tot){
        k-=ve[i].second;
        if (k<=0) return ve[i].first;
    }
    return 0;
}

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

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

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