CF896C Willem, Chtholly and Seniorious 珂朵莉树

问题描述

CF896C

LG-CF896C

---

题解

我expect就是T飞，从这里跳下去，也不碰和珂朵莉相关的任何东西。

珂朵莉树真好使。

珂朵莉树模板。

---

\(\mathrm{Code}\)

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

#define int long long
#define IT set<node>::iterator

template <typename Tp>
void read(Tp &x){
    x=0;char ch=1;int fh;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-'){fh=-1;ch=getchar(); }
    else fh=1;
    while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
    x*=fh;
}

struct node{
    int l,r;
    mutable int v;
    node(int L,int R=-1,int V=0):l(L),r(R),v(V){}
    bool operator<(node a) const{
        return l<a.l;
    }
};

int ksm(int x,int p,int mod){
    int res=1;x%=mod;
    while(p){
        if(p&1) res=res*x%mod;p>>=1;
        x=x*x%mod;
    }
    return res;
}

set<node>s;

IT split(int pos){
    IT it=s.lower_bound(node(pos));
    if(it!=s.end()&&it->l==pos) return it;
    --it;
    int L=it->l,R=it->r,V=it->v;
    s.erase(it);s.insert(node(L,pos-1,V));
    return s.insert(node(pos,R,V)).first;
}

void add(int l,int r,int val){
    IT rr=split(r+1),ll=split(l);
    for(;ll!=rr;++ll) ll->v+=val;
}

void assign(int l,int r,int val){
    IT rr=split(r+1),ll=split(l);
    s.erase(ll,rr);
    s.insert(node(l,r,val));
}

int rk(int l,int r,int k){
    vector <pair <int,int> > v;v.clear();
    IT rr=split(r+1),ll=split(l);
    for(;ll!=rr;ll++){
        v.push_back( pair<int,int> (ll->v,ll->r-ll->l+1));
    }
    sort(v.begin(),v.end());
    for( vector< pair<int,int> >::iterator it=v.begin();it!=v.end();it++){
        k-=it->second;
        if(k<=0) return it->first;
    }
}

int sum(int l,int r,int pos,int mod){
    IT rr=split(r+1),ll=split(l);
    int res=0;
    for(;ll!=rr;++ll) res=(res+(ll->r-ll->l+1)*ksm(ll->v,pos,mod))%mod;
    return res;
}

const int maxn=100007;

int n,m,seed,mx;
int a[maxn];
int random(){
    int res=seed;
    seed=(seed*7+13)%1000000007;
    return res;
}

signed main(){
    read(n);read(m);read(seed);read(mx);
    for(int i=1;i<=n;i++){
        a[i]=random()%mx+1;
        s.insert(node(i,i,a[i]));
    }
    s.insert(node(n+1,n+1,0));
    for(int i=1;i<=m;i++){
        int op=random()%4+1,l=random()%n+1,r=random()%n+1;
        if(l>r) swap(l,r);
        int x,y;
        if(op==3) x=random()%(r-l+1)+1;
        else x=random()%mx+1;
        if(op==4) y=random()%mx+1;
        if(op==1) add(l,r,x);
        else if(op==2) assign(l,r,x);
        else if(op==3) cout<<rk(l,r,x)<<endl;
        else cout<<sum(l,r,x,y)<<endl;
    }
    return 0;
}