bzoj 4311: 向量 线段树分治+凸包+三分

题意

你要维护一个向量集合，支持以下操作：
1.插入一个向量(x,y)
2.删除插入的第i个向量
3.查询当前集合与(x,y)点积的最大值是多少。如果当前是空集输出0
n<=200000 1<=x,y<=10^6

分析

如果直接用类似二维数点的方法来做显然不行。
考虑点积最大值的几何意义，也就是一条与该向量垂直的直线从无穷远处飞过来，碰到的第一个点就是答案。不难发现答案一定在凸包上面。
如果不带删除操作的话，我们可以用cdq分治。但如果有删除操作呢？那就要用线段树分治了。
具体来说就是对所有询问建线段树，然后把向量插入到对应的询问区间，最后每个询问在叶节点到根的每个点中查询一次就好了。
剩下的凸包+三分就是套路了。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;

typedef long long LL;

const int N=200005;

int n,tim[N],top;
LL ans[N];
struct pts{int x,y;}p[N],q[N],stack[N];
vector<pts> vec[N*4];

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();}
    return x*f;
}

bool cmp(pts a,pts b)
{
    return a.x<b.x||a.x==b.x&&a.y<b.y;
}

LL point(pts a,pts b)
{
    return (LL)a.x*b.x+(LL)a.y*b.y;
}

LL cross(pts c,pts a,pts b)
{
    return (LL)(a.x-c.x)*(b.y-c.y)-(LL)(a.y-c.y)*(b.x-c.x);
}

void ins(int d,int l,int r,int x,int y,pts p)
{
    if (l==x&&r==y) {vec[d].push_back(p);return;}
    int mid=(l+r)/2;
    if (y<=mid) ins(d*2,l,mid,x,y,p);
    else if (x>mid) ins(d*2+1,mid+1,r,x,y,p);
    else ins(d*2,l,mid,x,mid,p),ins(d*2+1,mid+1,r,mid+1,y,p);
}

void query(int id)
{
    pts x=q[id];
    int l=1,r=top;
    while (r-l>=3)
    {
        int mid1=l+(r-l)/3,mid2=r-(r-l)/3;
        if (point(x,stack[mid1])<=point(x,stack[mid2])) l=mid1;
        else r=mid2;
    }
    for (int i=l;i<=r;i++) ans[id]=max(ans[id],point(x,stack[i]));
}

void solve(int d,int l,int r)
{
    sort(vec[d].begin(),vec[d].end(),cmp);
    top=0;
    for (int i=0;i<vec[d].size();i++)
    {
        pts x=vec[d][i];
        while (top>1&&cross(stack[top-1],stack[top],x)>=0) top--;
        stack[++top]=x;
    }
    for (int i=l;i<=r;i++) if (q[i].x) query(i);
    if (l==r) return;
    int mid=(l+r)/2;
    solve(d*2,l,mid);solve(d*2+1,mid+1,r);
}

int main()
{
    n=read();int tot=0;
    for (int i=1;i<=n;i++)
    {
        int op=read(),x;
        if (op==1) p[++tot].x=read(),p[tot].y=read(),tim[tot]=i;
        else if (op==2) x=read(),ins(1,1,n,tim[x],i,p[x]),tim[x]=0;
        else q[i].x=read(),q[i].y=read();
    }
    for (int i=1;i<=tot;i++) if (tim[i]) ins(1,1,n,tim[i],n,p[i]);
    solve(1,1,n);
    for (int i=1;i<=n;i++) if (q[i].x) printf("%lld\n",ans[i]);
    return 0;
}