[HEOI2013] Segment

传送门

这道题其实是一个李超线段树的模板题
李超线段树（$lichaosegmenttree$），用来解决多个函数在某坐标位置上的最值的问题，其实李超线段树其实是线段树的一类题目的思想，而且比较好理解，但是极难写对

---

直接做显然不太好做，我们考虑标记永久化，让每个点储存的是当前区间的中点的优势最大线段，然后我们查询的时候查询所有包含$x=k$的区间，然后再找优势最大的线段就可以了

然后考虑对于一个新加入的线段如何更新

首先只有一个区间完整的被一条线段覆盖才能被更新
然后我们对于各种情况分类讨论

1.如果当前区间还没有线段，直接加入，然后结束，而不用下放，因为我们

2.如果当前线段在区间内全部的部分都在当前优势最大的线段的下方

直接$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♝ 的博客