bzoj3533

题意：
维护一个向量集合，在线支持以下操作：
A x y (|x|，|y| < =10^8)：加入向量(x，y);
Q x y l r (|x|，|y| < =10^8，1 < =L < =R < =T，其中T为已经加入的向量个数)询问第L个到第R个加入的向量与向量(x，y)的点积的最大值。
集合初始时为空。强制在线。
1 < =N < =4×10^5

#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<iostream>
#include<vector>
#define N 410000
#define pb push_back
#define LL long long
using namespace std;
struct node1{LL x,y;}p[N];
struct node{int l,r,lc,rc,mx,mi;bool b;vector<node1> v1,v2;}lt[2*N];
bool online;
int n,tl,tail;
LL ans,inf=1ll<<50;
inline int decode (int x) {return x ^ (ans & 0X7FFFFFFF);}
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    x*=f;
    if(online) x=decode(x);
    return x;
}
void init()
{
    n=read();char c;scanf("%c",&c);while(c==' ') scanf("%c",&c);
    if(c!='E') online=1;
}
void bt(int l,int r)
{
    int now=++tl;
    lt[now].l=l;lt[now].r=r;
    if(l<r)
    {
        int mid=(l+r)/2;
        lt[now].lc=tl+1;bt(l,mid);
        lt[now].rc=tl+1;bt(mid+1,r);
    }
}
double get_k(node1 t1,node1 t2)
{
    double x=t1.x-t2.x,y=t1.y-t2.y;
    return y/x;
}
void merge1(int now,int lc,int rc)
{
    int l1=lt[lc].v1.size(),l2=lt[rc].v1.size(),i1=0,i2=0;
    tail=0;
    while(i1<l1 || i2<l2)
    {
        node1 t;
        if(i1<l1 && i2<l2)
        {
            if(lt[lc].v1[i1].x==lt[rc].v1[i2].x)
            {
                if(lt[lc].v1[i1].y>lt[rc].v1[i2].y) t=lt[lc].v1[i1];
                else t=lt[rc].v1[i2];
                i1++;i2++;
            }
            else if(lt[lc].v1[i1].x<lt[rc].v1[i2].x) {t=lt[lc].v1[i1];i1++;}
            else {t=lt[rc].v1[i2];i2++;}
        }
        else if(i1<l1) {t=lt[lc].v1[i1];i1++;}
        else {t=lt[rc].v1[i2];i2++;}
        while(tail>1)
        {
            double k1=get_k(p[tail],t),k2=get_k(p[tail],p[tail-1]);
            if(k1>=k2) tail--;
            else break;
        }
        p[++tail]=t;
    }
    for(int i=1;i<=tail;i++) lt[now].v1.pb(p[i]);
}
void merge2(int now,int lc,int rc)
{
    int l1=lt[lc].v2.size(),l2=lt[rc].v2.size(),i1=0,i2=0;
    tail=0;
    while(i1<l1 || i2<l2)
    {
        node1 t;
        if(i1<l1 && i2<l2)
        {
            if(lt[lc].v2[i1].x==lt[rc].v2[i2].x)
            {
                if(lt[lc].v2[i1].y<lt[rc].v2[i2].y) t=lt[lc].v2[i1];
                else t=lt[rc].v2[i2];
                i1++;i2++;
            }
            else if(lt[lc].v2[i1].x<lt[rc].v2[i2].x) {t=lt[lc].v2[i1];i1++;}
            else {t=lt[rc].v2[i2];i2++;}
        }
        else if(i1<l1) {t=lt[lc].v2[i1];i1++;}
        else {t=lt[rc].v2[i2];i2++;}
        while(tail>1)
        {
            double k1=get_k(p[tail],t),k2=get_k(p[tail],p[tail-1]);
            if(k1<=k2) tail--;
            else break;
        }
        p[++tail]=t;
    }
    for(int i=1;i<=tail;i++) lt[now].v2.pb(p[i]);
}
void upd(int now)
{
    int lc=lt[now].lc,rc=lt[now].rc;
    lt[now].mx=max(lt[lc].mx,lt[rc].mx);
    lt[now].mi=min(lt[lc].mi,lt[rc].mi);
    if(lt[now].b==0 && (lt[lc].b&lt[rc].b))
    {
        lt[now].b=1;
        merge1(now,lc,rc);
        merge2(now,lc,rc);
    }
}
void change(int now,int k,int a,int b)
{
    int lc=lt[now].lc,rc=lt[now].rc,mid=(lt[now].l+lt[now].r)/2;
    if(lt[now].l==lt[now].r) 
    {lt[now].b=1;lt[now].mi=lt[now].mx=a;node1 t=(node1){a,b};lt[now].v1.pb(t);lt[now].v2.pb(t);return;}
    if(mid>=k) change(lc,k,a,b);
    else change(rc,k,a,b);
    upd(now);
}
void find(int now,int l,int r,int &mx,int &mi)
{
    int lc=lt[now].lc,rc=lt[now].rc,mid=(lt[now].l+lt[now].r)/2;
    if(lt[now].l==l && lt[now].r==r) {mx=lt[now].mx;mi=lt[now].mi;return;}
    if(mid>=r) find(lc,l,r,mx,mi);
    else if(l>mid) find(rc,l,r,mx,mi);
    else {int m1,m2;find(lc,l,mid,mx,mi);find(rc,mid+1,r,m1,m2);mx=max(mx,m1);mi=min(mi,m2);}
}
void make1(int now,int a,int b)
{
    double k=-(double)a/(double)b,kk,x,y;
    int siz=lt[now].v1.size();
    int l=1,r=siz-1,w=siz-1;
    while(l<=r)
    {
        int mid=(l+r)/2;
        x=lt[now].v1[mid-1].x-lt[now].v1[mid].x;
        y=lt[now].v1[mid-1].y-lt[now].v1[mid].y;
        kk=y/x;
        if(kk<=k) w=mid-1,r=mid-1;
        else l=mid+1;
    }
    LL t=(LL)lt[now].v1[w].x*a+lt[now].v1[w].y*b;
    ans=max(ans,t);
}
void make2(int now,int a,int b)
{
    double k=-(double)a/(double)b,kk,x,y;
    int siz=lt[now].v2.size();
    int l=1,r=siz-1,w=siz-1;
    while(l<=r)
    {
        int mid=(l+r)/2;
        x=lt[now].v2[mid-1].x-lt[now].v2[mid].x;
        y=lt[now].v2[mid-1].y-lt[now].v2[mid].y;
        kk=y/x;
        if(kk>=k) w=mid-1,r=mid-1;
        else l=mid+1;
    }
    LL t=(LL)lt[now].v2[w].x*a+lt[now].v2[w].y*b;
    ans=max(ans,t);
}
void query(int now,int l,int r,int a,int b)
{
    int lc=lt[now].lc,rc=lt[now].rc,mid=(lt[now].l+lt[now].r)/2;
    if(lt[now].l==l && lt[now].r==r)
    {
        if(b>0) make1(now,a,b);
        else make2(now,a,b);
        return;
    }
    if(mid>=r) query(lc,l,r,a,b);
    else if(l>mid) query(rc,l,r,a,b);
    else query(lc,l,mid,a,b),query(rc,mid+1,r,a,b);
}
void solve()
{
    bt(1,n);int k=0,l,r,a,b,mx,mi;
    for(int i=1;i<=n;i++)
    {
        char c;scanf("%c",&c);
        while(c!='A' && c!='Q') scanf("%c",&c);
        if(c=='A')
        {
            a=read();b=read();
            change(1,++k,a,b);
        }
        else
        {
            a=read();b=read();l=read();r=read();
            if(l>r) swap(l,r);if(l>k || l==0) continue;
            r=min(r,k);
            ans=-inf;
            if(b==0)
            {
                find(1,l,r,mx,mi);
                ans=max(a*mx,a*mi);
            }
            else query(1,l,r,a,b);
            printf("%lld\n",ans);
        }
        int oo=1;
    }
}
int main()
{
    init();
    solve();
    return 0;
}

题解：
先推一发点积的式子(虽然有斜率什么高逼格的方法但我不太会><)
询问向量(a,b)。两个向量(x1,y1),(x2,y2) $x1<x2$，什么时候选x2更优？
$x1a+y1b<x2a+y2b$
$-a<b\frac{y2-y1}{x2-x1}$
当b>0
$\frac{y2-y1}{x2-x1}>-\frac{a}{b}$
当b<0
$\frac{y2-y1}{x2-x1}<-\frac{a}{b}$
当b=0时找出y的区间最大最小值比较一下就好。。
观察$\frac{y2-y1}{x2-x1}$这是两点间的斜率，就是说两点谁比较优是与他们间的斜率有关的
按照套路，我们维护一个斜率单调的凸壳就可以二分了。。
显然，b>0时要维护斜率单调递减，b<0时要维护斜率单调递增。
注意只会在尾部加入，没有中间的修改，让我想起了以前做(乱搞)过的一道题。我是O(n)，正解是高逼格的二进制分组线段树。
那我们也用二进制分组线段树，线段树节点维护这个区间的两个凸壳。每个点只会往父亲合并一次，合并是O(n log n)的。合并时遇到两个x坐标相同，即没有斜率的点对时手动讨论一发就知道谁比较优了。而查询则要在log个凸壳上二分，是O($nlog^2n$)的。
二进制分组就是这个线段树节点包含区间已经塞满后再向父亲合并。