Sdoi2014 向量集

题目描述

题解：

码力太差重构之后才$A……$

首先求向量点积最大很容易想到凸包，

设已知$(x_0,y_0)$，求$(x,y)$满足$(x,y)*(x_0,y_0)>=(x',y')*(x_0,y_0)$

设$(x,y)*(x_0,y_0)=c$

那么$x*x_0+y*y_0=c$，$y=\frac(-x_0,y_0)*x+\frac(c,x_0)$

所以$x_0>0$时，$b$取最大，维护上凸包；

$x_0<0$时，$b$取最小，维护下凸包。

其实$x_0=0$不需要单独维护，放到任意一堆里二分都行。

由于询问与时间序有关，我们可以按时间建线段树，

每个节点挂一个凸包。

每次询问时从对应凸包二分查最大点积。

插入时并不需要每插一个点就把这个凸包拆开重构，

只需要在子树插满，即$r=x$时建立凸包即可。

代码：

#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = ;
typedef long long ll;
const ll inf = 0x3f3f3f3f3f3f3f3fll;
const double eps = 1e16;
int n,tot;
ll ans;
char s0[],s[];
struct Point
{
    ll x,y;
    Point(){}
    Point(ll x,ll y):x(x),y(y){}
    ll operator * (const Point&a)const{return x*a.x+y*a.y;}
    ll operator ^ (const Point&a)const{return y*a.x-x*a.y;}
    Point operator + (const Point&a)const{return Point(x+a.x,y+a.y);}
    Point operator - (const Point&a)const{return Point(x-a.x,y-a.y);}
    bool operator < (const Point&a)const{return x==a.x?y<a.y:x<a.x;}
};
vector<Point>tmp[N<<],c1[N<<],c2[N<<];
struct segtree
{
    void insert(int l,int r,int u,ll x,ll y,int z)
    {
        tmp[u].push_back(Point(x,y));
        if(z==r)
        {
            sort(tmp[u].begin(),tmp[u].end());
            for(int i=;i<tmp[u].size();i++)
            {
                while(c1[u].size()>&&((c1[u][c1[u].size()-]-c1[u][c1[u].size()-])^(tmp[u][i]-c1[u][c1[u].size()-]))<=)c1[u].pop_back();
                c1[u].push_back(tmp[u][i]);
                while(c2[u].size()>&&((c2[u][c2[u].size()-]-c2[u][c2[u].size()-])^(tmp[u][i]-c2[u][c2[u].size()-]))>=)c2[u].pop_back();
                c2[u].push_back(tmp[u][i]);
            }
        }
        if(l==r)return ;
        int mid = (l+r)>>;
        if(z<=mid)insert(l,mid,u<<,x,y,z);
        else insert(mid+,r,u<<|,x,y,z);
    }
    void query(int l,int r,int u,Point p,int ql,int qr)
    {
        if(l==ql&&r==qr)
        {
            if(p.y>=)
            {
                int L = ,R = c1[u].size()-,as = ;
                while(L<=R)
                {
                    int mid = (L+R)>>;
                    if(p*c1[u][mid]>p*c1[u][mid-])as=mid,L=mid+;
                    else R=mid-;
                }
                ans=max(ans,p*c1[u][as]);
            }else
            {
                int L = ,R = c2[u].size()-,as = ;
                while(L<=R)
                {
                    int mid = (L+R)>>;
                    if(p*c2[u][mid]>p*c2[u][mid-])as=mid,L=mid+;
                    else R=mid-;
                }
                ans=max(ans,p*c2[u][as]);
            }
            return ;
        }
        int mid = (l+r)>>;
        if(qr<=mid)query(l,mid,u<<,p,ql,qr);
        else if(ql>mid)query(mid+,r,u<<|,p,ql,qr);
        else query(l,mid,u<<,p,ql,mid),query(mid+,r,u<<|,p,mid+,qr);
    }
}tr;
int main()
{
    scanf("%d%s",&n,s0);
    for(int x,y,l,r,i=;i<=n;i++)
    {
        scanf("%s%d%d",s,&x,&y);
        if(s0[]!='E')x^=ans,y^=ans;
        if(s[]=='A')
        {
            tot++;
            tr.insert(,n,,x,y,tot);
        }else
        {
            scanf("%d%d",&l,&r);
            if(s0[]!='E')l^=ans,r^=ans;
            ans=-inf;
            tr.query(,n,,Point(x,y),l,r);
            printf("%lld\n",ans);
            ans&=0x7fffffff;
        }
    }
    return ;
}