李超线段树学习笔记

III.李超线段树

李超线段树是一种可以维护动态凸包的线段树。更准确地说，其可以支持的常规操作有两种：

1. 在平面直角坐标系中插入一条线段。

2. 询问在坐标系中的一个点 $(x,+\infty)$ 向下看，能看到的点的坐标。换句话说，是一条自无穷高处引下的垂线与所有线段的交点中最高的那个点。

其具体实现如下：

我们以 $x$ 坐标为下标，建一棵线段树（可能需要离散化）。然后，在线段树的每个区间上，我们维护所有完整地包含该区间的线段中，在该区间的中点处最高的那条。具体而言，当我们插入一条线段时，首先要把它在线段树上拆作众多小区间。当递归到某个完整的区间时，考虑将其与该区间上本来记录的中点处最高线段做比较。明显，共有三种可能：

1. 新插入的线段严格不低于原本线段。此时直接替换就行。
2. 新插入的线段严格不高于原本线段。此时直接不替换就行。
3. 新插入的线段与原本线段有交。此时，先在当前区间中记录中点更高的那条，然后明显，在其一个子区间内，更高的一定是当前区间中记录的线段，可以不管；而在其另一个子区间中，更高的就不一定是它了，此时要递归入该子区间继续修改。

下面考虑分析复杂度。显然，如果插入的全部是直线，复杂度为 $n\log n$，因为两个子区间中最多访问其中一个；而如果我们插入的全是线段，因为要先拆到各个子区间上，所以复杂度就为 $n\log^2n$。

然后询问时就仿照标记永久化，直接在一路上的区间处用大区间上记录的线段更新即可。明显是 $O(n\log n)$ 的。

III.I.[JSOI2008]Blue Mary开公司

是板题，甚至插入的都不是线段而是直线，就更好处理了。

真的很好写。

代码：

#include<bits/stdc++.h>
using namespace std;
const int lim=50000;
#define lson x<<1
#define rson x<<1|1
#define mid ((l+r)>>1)
struct SegTree{double lv,rv;}seg[200100];
double calc(double lv,double rv,int l,int r,int P){return lv+(rv-lv)/(r-l)*(P-l);}
void modify(int x,int l,int r,double LV,double RV){
	if(LV+RV>seg[x].lv+seg[x].rv)swap(seg[x].lv,LV),swap(seg[x].rv,RV);
	if(l==r)return;
	if(LV>seg[x].lv)modify(lson,l,mid,LV,calc(LV,RV,l,r,mid));
	if(RV>seg[x].rv)modify(rson,mid+1,r,calc(LV,RV,l,r,mid+1),RV);
}
double query(int x,int l,int r,int P){
	if(l>P||r<P)return 0;
	double ret=calc(seg[x].lv,seg[x].rv,l,r,P);
	if(l!=r)ret=max(ret,query(lson,l,mid,P)),ret=max(ret,query(rson,mid+1,r,P));
	return ret;
}
char s[10];
int T;
int main(){
	scanf("%d",&T);
	for(int i=1;i<=T;i++){
		scanf("%s",s);
		if(s[0]=='Q'){int x;scanf("%d",&x);printf("%d\n",(int)(query(1,1,lim,x)/100));}
		else{double l,r;scanf("%lf%lf",&l,&r),r=l+(lim-1)*r;modify(1,1,lim,l,r);}
	}
	return 0;
}