问题描述
题解
我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;
}