利用线段树维护凸壳（李超树）

Text

之前听lyd_729神犇讲过关于李超树的思想，大概是标记永久化一类的

要求支持动态插入直线（线段），动态查询上/下凸壳与x=k的交点，k为整数
今天做模拟赛的时候刚好碰到一题，然而忘得差不多了
考场上自己YY了一个：

线段树每一个区间都保存一条直线，表示在这个区间中这条直线有可能有一部分为凸壳中的边

当我们需要插入一条直线（线段），与当前线段树区间保存的直线进行比较，看哪一个斜率大，哪一个更优的部分占整个区间更多（就是比较区间中点和直线交点）

将占整个区间更多的保存为当前区间的直线，另外一条则递归到相应的儿子去

如果一条直线在这个区间内完全都比另一条优，那就可以直接保存返回了

比完赛看了一下网上关于李超树的讲解，觉得自己真是个SB
其实可以直接比较这两条直线哪个斜率大，比较这两条在区间中点时的值然后直接判断即可，就不需要计算交点了

模板

相当的丑
大佬们将就着看吧。。

LL calc(LL k,LL h,LL w)
{
    return k*w+h;
}
int nwp(int k)
{
    if(k) return k;
    tr[++n1].l=tr[n1].r=tr[n1].k=tr[n1].h=0;
    return n1;
}
LL get(int k,int l,int r,int t)
{
    if(l==r&&l==t) return calc(tr[k].k,tr[k].h,t);
    int mid=(l+r)/2;
    LL v=0;
    tr[k].l=nwp(tr[k].l),tr[k].r=nwp(tr[k].r);
    if(t<=mid) v=get(tr[k].l,l,mid,t);
    else v=get(tr[k].r,mid+1,r,t);
    return max(calc(tr[k].k,tr[k].h,t),v);
}
double dot(LL k1,LL h1,LL k2,LL h2)
{
    return (h2-h1)*1.0/(k1-k2);
}
void ins(int k,int l,int r,LL dx,LL hi)
{
    if(l==r)
    {
        if(calc(dx,hi,l)>calc(tr[k].k,tr[k].h,l)) tr[k].k=dx,tr[k].h=hi;
        return;
    }
    if(calc(dx,hi,l)>=calc(tr[k].k,tr[k].h,l)&&calc(dx,hi,r)>=calc(tr[k].k,tr[k].h,r)) 
    {
        tr[k].k=dx,tr[k].h=hi;
        return;
    }
    if(calc(dx,hi,l)<calc(tr[k].k,tr[k].h,l)&&calc(dx,hi,r)<calc(tr[k].k,tr[k].h,r)) return;
    int mid=(l+r)/2;
    tr[k].l=nwp(tr[k].l),tr[k].r=nwp(tr[k].r);
    double w=dot(dx,hi,tr[k].k,tr[k].h);
    bool pd=0;
    if(dx>tr[k].k)
    {
        if(w>mid) pd=1;
        else swap(tr[k].k,dx),swap(tr[k].h,hi);
    }
    else
    {
        if(w>=mid+1) swap(tr[k].k,dx),swap(tr[k].h,hi),pd=1; 
    }
    if(!pd) ins(tr[k].l,l,mid,dx,hi);
    else ins(tr[k].r,mid+1,r,dx,hi);
}