玄学数据结构——珂朵莉树

关于珂朵莉

こんなにも、たくさんの幸せをあの人に分けてもらった
だから、きっと
今の、私は
誰が何と言おうと

世界一、 幸せな女の子だ

————クトリ

美图欣赏：

(图片来源：萌娘百科)

中国珂学院

好了好了说正事

关于这奇怪的名字

本来是起源于CF896C Willem, Chtholly and Seniorious(一个CF在国内的镜像网站 )，当时CF上有一名叫做Old Driver的用户(就是lxl)在提交了线段树的正解之后又补充了这么一份奇怪的数据结构，得以命名为珂朵莉树(又称ODT(老司机树))

在某谷上的评级是一道黑题，但是日常被爆切

前置知识

std::set的简单操作，一点区间操作的知识，没了

正篇

珂朵莉树可以用来解决含有区间覆盖(或者是区间统一赋值)的操作，本质上是基于(暴力)std::set的数据结构，对于随机数据来说效率很高甚至有时候比正解还高

注意这种数据结构只在数据随机的情况下表现良好

珂朵莉树的定义非常简单，就是把一段有着相同值的区间变成set里面的一个元素，以左区间位置为关键字维护的一个set。所以我们就可以写出珂朵莉树的节点定义和初始化：

struct node {
	int l, r;//左右区间
	mutable ll val;//值，这里的mutable意味着在集合的任何位置上都是课修改的
	node(int l_ = -1, int r_ = -1, ll val_ = 0) {//构造函数
		l = l_, r = r_, val = val_;
	}
	bool operator < (const node &a) const {//运算符重载
		return l < a.l;
	}
};
set<node> st;//搞一个集合

split操作

split操作是珂朵莉树里面最重要的操作，其实很简单。对于某一个位置pos,我们将其所处的区间分为[l,pos),[pos,r]并且返回后面那个区间的迭代器，具体的看代码吧：

IT split(int pos) {
	IT it = st.lower_bound(node(pos));//查找后继
	if (it != st.end() && it->l == pos) return it;//如果此时这个区间的左边界就是pos的话就直接返回
	--it;//否则就是前面的区间
	node tmp = *it;//保存这个区间的信息
	st.erase(it);//删除这个区间
	st.insert(node(tmp.l, pos - 1, tmp.val));//插入两个新区间
	return st.insert(node(pos, tmp.r, tmp.val)).first;
}

assign操作

assign操作就是区间赋值/区间推平(像远野志贵那样)，有了split操作之后我们就可以很简单地写出assign的代码了。

珂朵莉树的复杂度全部来源于assign操作，如果我们只有split操作的化复杂度就螺旋上天

由于数据随机，所以有1/4的操作都是assign, assign操作会合并具有相同值的区间，使set的大小下降，最后趋近于logn, 所以时间复杂度是mlogn

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

这里要注意一点，就是但凡我们要split的话肯定是优先split(r+1)，这里应该是为了保证不让在split(l)的时候把r的位置也给搞掉了。

接下来就是一些非常暴力的运算方式了，这里也给一下代码：

区间加：

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

这里要注意，由于我们set存储的是区间，所以在解决区间求和之类的问题的时候不要忘记考虑区间长度

区间求和：

ll querySum(int l, int r) {
	IT itr = split(r + 1), itl = split(l);
	ll res = 0;
	for (IT it = itl; it != itr; ++it) {
		res += (it->r - it->l + 1) * it->val;
	}
	return res;
}

区间第k小：

ll queryKth(int l, int r, int k) {
	vector< pair<ll, int> > vec(0);
	IT itr = split(r + 1), itl = split(l);
	for (IT it = itl; it != itr; ++it) {
		vec.push_back(make_pair(it->val, it->r - it->l + 1));
	}
	sort(vec.begin(), vec.end());
	for (vector< pair<ll, int> >::iterator it = vec.begin(); it != vec.end(); ++it) {
		if ((k -= it->second) <= 0) return it->first;
	}
	return -1;
}

这里我们其实就是直接把整个区间都push到一个vector里面，然后再暴力sort接着输出。。。注意我们要把值放在pair的第一个关键字哪里，因为我们是要给值排序的嘛

这里给出一个板子：

//Chtholly tree
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
struct node {
	int l, r;
	mutable ll val;
	node(int l_ = -1, int r_ = -1, ll val_ = 0) {
		l = l_, r = r_, val = val_;
	}
	bool operator < (const node &a) const {
		return l < a.l;
	}
};
typedef set<node>::iterator IT;
set<node> st;
IT split(int pos) {
	IT it = st.lower_bound(node(pos));
	if (it != st.end() && it->l == pos) return it;
	--it;
	node tmp = *it;
	st.erase(it);
	st.insert(node(tmp.l, pos - 1, tmp.val));
	return st.insert(node(pos, tmp.r, tmp.val)).first;
}

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

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

ll querySum(int l, int r) {
	IT itr = split(r + 1), itl = split(l);
	ll res = 0;
	for (IT it = itl; it != itr; ++it) {
		res += (it->r - it->l + 1) * it->val;
	}
	return res;
}

ll queryKth(int l, int r, int k) {
	vector< pair<ll, int> > vec(0);
	IT itr = split(r + 1), itl = split(l);
	for (IT it = itl; it != itr; ++it) {
		vec.push_back(make_pair(it->val, it->r - it->l + 1));
	}
	sort(vec.begin(), vec.end());
	for (vector< pair<ll, int> >::iterator it = vec.begin(); it != vec.end(); ++it) {
		if ((k -= it->second) <= 0) return it->first;
	}
	return -1;
}

int main(void) {
	
}

复杂度证明

虽然前面我们估算(瞎猜)了一下含有assign操作复杂度，但是这里我找了一份相对比较完全的复杂度证明

复杂度证明
(在骗分数据结构底下寻找复杂度证明是不是搞错了什么)

后记

还有很多别的操作，其本质就是对分离出的那一段区间进行暴力运算，这里就不再赘述了

题目的话这里有一道：P4344 [SHOI2015]脑洞治疗仪
2020/03/09的时候还没有卡珂朵莉树，但是如果要用珂朵莉树的话需要吸氧(开O2)才能过

好多题目珂朵莉树都没法AC(可见还是好好打线段树才是王道)

总的来说，珂朵莉树就是用来骗分的(一题骗到84分还是很可观的)，学会打了之后可以应对某些题目的时候拿稍微多一点的分数

可能有些地方讲的有点问题，欢迎dalao提出