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;
}