这道题其实是一个李超线段树的模板题
李超线段树(
l
i
c
h
a
o
s
e
g
m
e
n
t
t
r
e
e
lichao\ segment\ tree
lichao segment tree),用来解决多个函数在某坐标位置上的最值的问题,其实李超线段树其实是线段树的一类题目的思想,而且比较好理解,但是极难写对
直接做显然不太好做,我们考虑标记永久化,让每个点储存的是当前区间的中点的优势最大线段,然后我们查询的时候查询所有包含 x = k x=k x=k的区间,然后再找优势最大的线段就可以了
然后考虑对于一个新加入的线段如何更新
首先只有一个区间完整的被一条线段覆盖才能被更新
然后我们对于各种情况分类讨论
1.如果当前区间还没有线段,直接加入,然后结束,而不用下放,因为我们
2.如果当前线段在区间内全部的部分都在当前优势最大的线段的下方
直接
r
e
t
u
r
n
return
return掉,不管他
3.如果当前线段的在区间的全部部分全部低于当前优势最大线段
用当前的线段覆盖原来的优势最大线段
4.如果当前线段和原来优势线段在这一段区间内有交点,分四种情况进行讨论
(1)新线段的左端点高于原优势最大线段,并且交点位于区间中点左侧
他不能成为这个区间的优势最大线段
所以尝试拿这个线段更新左子区间
(2)新线段的左端点高于原优势最大线段,并且交点位于区间中点右侧
那么他在这个区间就是最强的了,拿他覆盖原来的线段,但是注意这里,因为原来的优势最大的线段在交点右边还有可能是最好的,所以尝试拿原来的做优线段去更新右子树
(3)新线段的右端点高于原优势最大线段,并且交点位于区间中点右侧
同理,利用新线段更新右子区间
(4)新鲜度胺的右端点高于原优势最大线段,并且交点位于区间中点左侧
和(2)同理,覆盖当前区间,拿优势最大线段尝试更新左子区间
然后就是代码了,代码细节还是挺多的,还是挺恶心的(
#include <bits/stdc++.h>
using namespace std;
# define Rep(i,a,b) for(int i=a;i<=b;i++)
# define _Rep(i,a,b) for(int i=a;i>=b;i--)
# define RepG(i,u) for(int i=head[u];~i;i=e[i].next)
typedef long long ll;
const int N=1e5+5;
template<typename T> void read(T &x){
x=0;int f=1;
char c=getchar();
for(;!isdigit(c);c=getchar())if(c=='-')f=-1;
for(;isdigit(c);c=getchar())x=(x<<1)+(x<<3)+c-'0';
x*=f;
}
int n,segt;
int lastans;
struct seg{
double k,b;
double func(int x){
return k*x+b;
}
}a[N];
struct node{
int l,r;
int id;
}lic[N<<2];
# define lc (u<<1)
# define rc (u<<1|1)
int lichaomax(int x,int y,int p){
if(!x)return y;
if(!y)return x;
return a[x].func(p)>a[y].func(p)?x:y;
}
void build(int u,int l,int r){
lic[u].l=l,lic[u].r=r;
if(l==r)return;
int mid=l+r>>1;
build(lc,l,mid);
build(rc,mid+1,r);
}
void update(int u,int l,int r,int id){
int mid=lic[u].l+lic[u].r>>1;
if(l<=lic[u].l&&r>=lic[u].r){
if(!lic[u].id){lic[u].id=id;return;}
double yl=a[lic[u].id].func(lic[u].l),yr=a[lic[u].id].func(lic[u].r);
double xl=a[id].func(lic[u].l),xr=a[id].func(lic[u].r);
if(xl<=yl&&xr<=yr)return;
if(xl>=yl&&xr>=yr){lic[u].id=id;return;}
double cross=(a[id].b-a[lic[u].id].b)/(a[lic[u].id].k-a[id].k);
if(xl>=yl){
if(cross<=mid)update(lc,l,r,id);
else update(rc,l,r,lic[u].id),lic[u].id=id;
}
else{
if(cross> mid)update(rc,l,r,id);
else update(lc,l,r,lic[u].id),lic[u].id=id;
}
return;
}
if(l<=mid)update(lc,l,r,id);
if(r>mid)update(rc,l,r,id);
}
int query(int u,int x){
if(lic[u].l==lic[u].r)return lic[u].id;
int mid=lic[u].l+lic[u].r>>1;
int res=lic[u].id;
if(x<=mid)res=lichaomax(res,query(lc,x),x);
else res=lichaomax(res,query(rc,x),x);
return res;
}
int main()
{
read(n);
build(1,1,40000);
Rep(i,1,n){
int opt;
read(opt);
if(!opt){
int k;
read(k);
k=(k+lastans-1)%39989+1;
lastans=query(1,k);
printf("%d\n",lastans);
}
else{
int x1,y1,x2,y2;
read(x1),read(y1),read(x2),read(y2);
x1=(x1+lastans-1)%39989+1;
x2=(x2+lastans-1)%39989+1;
y1=(y1+lastans-1)%1000000000+1;
y2=(y2+lastans-1)%1000000000+1;
segt++;
if(x1==x2){
a[segt]=(seg){0,max(y1,y2)};
update(1,x1,x1,segt);
}
else{
a[segt].k=1.*(y2-y1)/(x2-x1);
a[segt].b=1.*(y1-1.*x1*a[segt].k);
update(1,min(x1,x2),max(x1,x2),segt);
}
}
}
return 0;
}
图片转载自♗Wendigo♝ 的博客