珂朵莉树(Chtholly Tree)学习笔记

---

珂朵莉树原理

其原理在于运用一颗树（set，treap，splay......）其中要求所有元素有序，并且支持基本的操作（删除，添加，查找......）来实现区间压缩。

那么区间压缩的意义在于区间推平这是珂朵莉树的核心（如果没有这个操作实际上不一定需要这种算法）

ps：若保证有连续相等甚至递增的区间，也可以的（吧？)。

可想而知它的操作在于对区间的分裂和合并操作

（为什么？因为这样可以方便而快捷的区间推平）

珂朵莉树的实现

在众多树中因为set这个c++自带，所以决定选择它。

结构

一个node结构体——把它叫区间(已typedef long long LL)

struct node

{

	int l,r;LL v;//这里官方写法加了一个mutable，看个人写法

	node(int L,int R,LL V):l(L),r(R),v(V){}

	bool operator < (const node& o) const

	{

		return l<o.l;

	}

};

声明集合

set<node>s;

初始化

for(int i=1;i<=n;++i)

    s.insert(i,i,val);

s.insert(n+1,n+1,0);//见下面中it!=s.end()的前提

分裂

#define it_ set<node>::iterator

it_ split(int pos)

{

	it_ it=s.lower_bound(node(pos));

	if(it!=s.end()&& it->l==pos) return it;

	--it;

	int ll=it->l,rr=it->r;

	LL vv=it->v;

	s.erase(it);

	s.insert(node(ll,pos-1,vv));

	return s.insert(node(pos,rr,vv)).first;

}

区间推平

void assign(int l,int r,LL val=0)

{

	it_ itl=split(l),itr=split(r+1);

	s.erase(itl,itr);

	s.insert(node(l,r,val));

}

如上。

---

栗例子

见cf896c

仅仅给出代码

#include <bits/stdc++.h>

#define it_ set<node>::iterator

// #define swap(X,Y) {int (tmp)=(X),(X)=(Y),(Y)=(tmp);}

using namespace std;

typedef long long LL;

const int M7=1e9+7;

const int maxn=1e5+5;

int n,m;

LL seed,vmax;

LL a[maxn];

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 pos)

{

	it_ it=s.lower_bound(node(pos));

	if(it!=s.end()&& it->l==pos) return it;

	--it;

	int ll=it->l,rr=it->r;

	LL vv=it->v;

	s.erase(it);

	s.insert(node(ll,pos-1,vv));

	return s.insert(node(pos,rr,vv)).first;

}

void assign(int l,int r,LL val=0)

{

	it_ itl=split(l),itr=split(r+1);

	s.erase(itl,itr);

	s.insert(node(l,r,val));

}

void add(int l,int r,LL val=0)

{

	it_ itl=split(l),itr=split(r+1);

	for(it_ i=itl;i!=itr;++i)

	{

		i->v+=val;

	}

}

LL so(int l,int r,int k)

{

	it_ itl=split(l),itr=split(r+1);

	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());

	for(vector<pair<LL,int>>::iterator it=v.begin();it!=v.end();++it)

	{

		k-=it->second;

		if(k<=0)return it->first;

	}

}

LL pow(LL a,LL b,LL mod)

{

	LL ret=1;a%=mod;

	while(b)

	{

		if(b&1) ret=ret*a%mod;

		a=a*a%mod;

		b>>=1;

	}

	return ret%mod;

}

LL ok(int l,int r,LL k,LL mod)

{

	it_ itl=split(l),itr=split(r+1);

	LL ans=0;

	for(;itl!=itr;++itl)

	{

		ans+=(LL)(itl->r-itl->l+1)*pow(itl->v,k,mod)%mod;

		ans%=mod;

	}

	return ans%mod;

}

LL rnd()

{

	LL ret=seed;

	seed=(seed*7+13)%M7;

	return ret;

}

int main(int argc, char const *argv[])

{

	ios::sync_with_stdio(0);

	cin>>n>>m>>seed>>vmax;

	for(int i=1;i<=n;++i)

	{

		a[i]=(rnd()%vmax)+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=(rnd()%4)+1,

		l=(rnd()%n)+1,

		r=(rnd()%n)+1;LL x,y;

		if(l>r)

			swap(l,r);

		if(op==3)

			x=(rnd()%(LL)(r-l+1))+1;

		else

			x=(rnd()%vmax)+1;

		if(op==4)

			y=(rnd()%vmax)+1;

		if(op==1) add(l,r,x);

		else if(op==2) assign(l,r,x);

		else if(op==3) cout<<so(l,r,x)<<endl;

		else if(op==4) cout<<ok(l,r,x,y)<<endl;

	}

	return 0;

}