【算法】珂朵莉树板子

 https://www.luogu.com.cn/problem/CF896C <---例子

#include<bits/stdc++.h>
using namespace std;
#define SET0(a) memset(a,0,sizeof(a))
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define DWN(i,a,b) for(int i=(a);i>=(b);i--)
typedef long long ll;
#define IT set<node>::iterator
int n,m;
ll seed,vmax;
const int maxn=1e5+5;
const int mod=1e9+7;
ll a[maxn];
ll rd(){
    ll ret=seed;
    seed=(seed*7+13)%mod;
    return ret;
}
ll fpow(ll b,ll p,ll mo){
    ll res=1; b%=mo;
    for(;p;p>>=1,b=(b*b)%mo){
        if(p&1) res=(res*b+mo)%mo;
		res=(res+mo)%mo;
    }
    return (res+mo)%mo;
}

struct node{
    int l,r;
    mutable ll v;
    node(int L,int R=-1,ll V=0):l(L),r(R),v(V){};//构造函数
    bool operator<(const node&o)const{//运算符重载
        return l<o.l;
    }
};
set<node>s;
//分裂
IT split(int p){
    IT it=s.lower_bound(node(p));
    if(it->l==p) return it; //如果p在it对应的区间里面
    it--; //否则找上一个it（p必然在相应区间上）

    //模拟分裂区间
    int L=it->l; int R=it->r; ll V=it->v;
    s.erase(it);
    s.insert(node(L,p-1,V));
    return s.insert(node(p,R,V)).first; //insert会返回一个pair类型 first是迭代器 second是bool类型
}
//推平
void assign(int l,int r,ll v=0){
    IT itr=split(r+1); IT itl=split(l);
    s.erase(itl,itr);
    s.insert(node(l,r,v));
}
//sum up
void add(int l,int r,ll v=1){
    IT itr=split(r+1); IT itl=split(l);
    for(;itl!=itr;itl++){
        itl->v+=v;
    }
}
//rank
ll rk(int l,int r,int k){
    IT itr=split(r+1); IT itl=split(l);
    //建立一个vector取出区间 first存区间值 second存区间长
    vector<pair<ll,int> >v;
    v.clear();

    for(;itl!=itr;itl++){
        v.push_back(pair<ll,int>(itl->v,itl->r-itl->l+1)); //
    }
    sort(v.begin(),v.end());
    vector<pair<ll,int> >::iterator vit;
    for(vit=v.begin();vit!=v.end();vit++){
        k-=vit->second;
        if(k<=0) return vit->first;
    }
    return -1;
}
//幂次和
ll sum_pow(int l,int r,int k,int mo){
    ll res=0;
    IT itr=split(r+1); IT itl=split(l);
    for(;itl!=itr;itl++){
        res=((itl->r-itl->l+1)*fpow(itl->v,ll(k),ll(mo))+res)%mo;
        res%=mo;
    }
    return res;
}
int main(){
    cin>>n>>m>>seed>>vmax;
    FOR(i,1,n){
        a[i]=rd()%vmax+1;
        s.insert(node(i,i,a[i])); //初始化
    }
    s.insert(node(n+1,n+1,0));

    FOR(i,1,m){
        int op=rd()%4 +1;
        int l,r;
        l=rd()%n +1;
        r=rd()%n +1;
        if(l>r) swap(l,r);

        int x,y;

        if(op==3) x=rd()%(r-l+1)+1;
        else x=rd()%vmax+1;

        if(op==4) y=rd()%vmax+1;

        switch (op)
        {
        case 1:
            add(l,r,x);
            break;
        
        case 2:
            assign(l,r,x);
            break;
        
        case 3:
            cout<<rk(l,r,x)<<endl;
            break;
        
        case 4:
            cout<<sum_pow(l,r,x,y)<<endl;
            break;

        default:
            break;
        }
    }
    return 0;
}