标签：Willem itl int CF896C ll odt Seniorious ret split

写一种数据结构，支持：

• \(1\) \(l\) \(r\) \(x\) ：将\([l,r]\) 区间所有数加上\(x\)
• \(2\) \(l\) \(r\) \(x\) ：将\([l,r]\) 区间所有数改成\(x\)
• \(3\) \(l\) \(r\) \(x\) ：输出将\([l,r]\) 区间从小到大排序后的第\(x\) 个数是的多少(即区间第\(x\) 小，数字大小相同算多次，保证 \(1\leq\) \(x\) \(\leq\) \(r-l+1\) )
• \(4\) \(l\) \(r\) \(x\) \(y\) ：输出\([l,r]\) 区间每个数字的\(x\) 次方的和模\(y\) 的值(即(\(\sum^r_{i=l}a_i^x\) ) \(\mod y\) )

其中数据在本机上根据程序生成（生成器是随机数据）， \(n,m\leq 10^5,x,y\leq 10^9\) 。

珂朵莉树的元老题。首先声明珂朵莉树只有在随机数据下保证复杂度正确，大多数情况下只能用来骗分。然后用这一题来写一下ODT的基本流程。

首先珂朵莉树并不是树，只是借用了STL::set来解决区间问题，其中set中存储的结构体 \((l,r,v)\) 表示 \([l,r]\) 区间存储的值为 \(v\) 。同时支持合并 \(assign\) 和分割 \(split\) 操作，分别表示将 \([l,r]\) 区间的值全部附为 \(v\) ，将 \([l,r]\) 区间改为 \([l,x),[x,r]\) 两段，并返回后一段的指针。然后其他操作直接在STL上暴力 \(++it\) 即可。

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

namespace QuickMul{
	ll Qmul(ll a,ll b,ll p){
		ll ret=1;
		while(b){
			if(b&1)
				ret=(ret*a)%p;
			b>>=1;
			a=((a%p)*(a%p))%p;
		}
		return ret;
	}
} 

namespace ChthollyTree{
	struct node_c{
		ll l,r;
		mutable ll v;
		node_c(ll il,ll ir,ll iv){
			l=il,r=ir,v=iv;
		}
		inline bool operator <(const node_c o)const{
			return l<o.l;
		}
	};
	struct rank{
		ll num,cnt;
		bool operator < (const rank a)const{
			return num<a.num;
		}
		rank(ll inum,ll icnt){
			num=inum;
			cnt=icnt;
		}
		rank(){
			num=0;
			cnt=0;
		}
	};
	#define iter set<node_c>::iterator
	set<node_c> odt;
	iter split(ll x){
		iter it=odt.lower_bound(node_c(x,0,0));
		if(it!=odt.end() && it->l == x){
			return it;
		}
		it--;
		if(it->r < x)
			return odt.end();
		ll l=it->l;
		ll r=it->r;
		ll v=it->v;
		odt.erase(it);
		odt.insert(node_c(l,x-1,v));
		return odt.insert(node_c(x,r,v)).first; 
	}
	void assign(ll l,ll r,ll v){
		iter itr=split(r+1),itl=split(l);
		odt.erase(itl,itr);
		odt.insert(node_c(l,r,v));
	}
	void add(int l,int r,int x){
		iter itr=split(r+1),itl=split(l);
		for(;itl!=itr;itl++){
			itl->v+=x;
		}
	}
	ll getkth(int l,int r,int k){
		iter itr=split(r+1),itl=split(l);
		vector<rank> v;
		for(;itl!=itr;itl++){
			v.push_back(rank(itl->v,(itl->r)-(itl->l)+1));
		}
		sort(v.begin(),v.end());
		for(int i=0;i<(int)v.size();++i){
			if(v[i].cnt>=k)
				return v[i].num;
			else
				k-=v[i].cnt;
		}
	}
	ll getpow(ll l,ll r,ll x,ll y){
		iter itr=split(r+1),itl=split(l);
		ll ret=0;
		for(;itl!=itr;++itl){
			ret=(ret+((((itl->r)-(itl->l)+1)%y)*(QuickMul::Qmul(itl->v,x,y))%y))%y;
		}
		return ret;
	}
}
using namespace ChthollyTree;
ll n,m,seed,vmax,a[100005];

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

int main(){
	cin>>n>>m>>seed>>vmax;
	for(int i=1;i<=n;++i)
		a[i]=(rnd()%vmax)+1;
	for(int i=1;i<=n;++i){
		odt.insert(node_c(i,i,a[i]));
	}
	for(int i=1;i<=m;++i){
		ll op=(rnd()%4)+1;
		ll l=(rnd()%n)+1;
		ll r=(rnd()%n)+1;
		ll x,y;
		if(l>r)
			swap(l,r);
		if(op==3){
			x=(rnd()%(r-l+1))+1;
			printf("%lld\n",getkth(l,r,x));
		}
		else{
			x=(rnd()%vmax)+1;
		}
		if(op==4){
			y=(rnd()%vmax)+1;
			printf("%lld\n",getpow(l,r,x,y));
		}
		if(op==1)
			add(l,r,x);
		if(op==2)
			assign(l,r,x); 
	}
	
}

标签：Willem,itl,int,CF896C,ll,odt,Seniorious,ret,split
来源： https://www.cnblogs.com/zhouzizhe/p/16639839.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。