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);
}