【SDOI2014】【BZOJ3533】向量集

最新推荐文章于 2019-03-06 19:29:00 发布

CreationAugust

最新推荐文章于 2019-03-06 19:29:00 发布

阅读量1.6k

点赞数

分类专栏：随便搞搞文章标签：线段树凸包

本文链接：https://blog.csdn.net/CreationAugust/article/details/50487792

版权

随便搞搞专栏收录该内容

403 篇文章 0 订阅

订阅专栏

Description

维护一个向量集合，在线支持以下操作：
“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)的点积的最大值。
集合初始时为空。

Input

输入的第一行包含整数N和字符s，分别表示操作数和数据类别；
接下来N行，每行一个操作，格式如上所述。
请注意s≠'E'时，输入中的所有整数都经过了加密。你可以使用以下程序

得到原始输入：
inline int decode (int x long long lastans) {
return x ^ (lastans & Ox7fffffff);
}
function decode
begin
其中x为程序读入的数，lastans为之前最后一次询问的答案。在第一次询问之前，lastans=0。

注：向量(x，y)和(z，W)的点积定义为xz+yw。

Output

对每个Q操作，输出一个整数表示答案。

Sample Input

6 A

A 3 2

Q 1 5 1 1

A 15 14

A 12 9

Q 12 8 12 15

Q 21 18 19 18

Sample Output

解释：解密之后的输入为

6 E

HINT

1 < =N < =4×10^5

新加数据一组..2015.315

Source

Round 1 Day 2

线段树维护一下凸包然后三分
也可以二进制分组
才知道BZOJ上WA有时候也会显示RE…

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 400010
#define GET (ch>='0'&&ch<='9')
#define LL long long
#define MAXINT 0x7fffffff
#define MAXLL 1ll<<60
#define lchild rt<<1,l,mid
#define rchild rt<<1|1,mid+1,r
#define ln rt<<1
#define rn rt<<1|1
#define MAXSIZE 50
using namespace std;
int n,cnt,top,top1,top2;
LL lastans;
char Flag,ch;
struct Point
{
    int x,y;
    friend Point operator -(Point a,Point b)    {return (Point){a.x-b.x,a.y-b.y};}
    friend LL operator *(Point a,Point b)   {return (LL)a.x*b.x+(LL)a.y*b.y;}
    friend LL operator ^(Point a,Point b)   {return (LL)a.x*b.y-(LL)a.y*b.x;}
}s[MAXN],tmp[MAXN],p,sta1[MAXN<<3],sta2[MAXN<<3];
bool cmp1(Point a,Point b)  {return a.x==b.x?a.y>b.y:a.x<b.x;}
bool cmp2(Point a,Point b)  {return a.x==b.x?a.y<b.y:a.x<b.x;}
struct seg  {int l,r,l1,r1,l2,r2;}tree[MAXN<<2];
void in(int &x)
{
    char ch=getchar();x=0;int flag=1;
    while (!GET)    flag=ch=='-'?-1:1,ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();x*=flag;
}
void build(int rt=1,int l=1,int r=n)
{
    tree[rt].l=l;tree[rt].r=r;
    if (l==r)   return;
    int mid=(l+r)>>1;build(lchild);build(rchild);
}
LL Dir(Point a,Point b,Point c) {Point A=b-a,B=c-b;return B^A;}//有向面积 
void insert(int rt)
{
    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;
    if (R==cnt&&R-L>=MAXSIZE)
    {
        top=0;
        for (int i=L;i<=R;i++)  tmp[++top]=s[i];
        sort(tmp+1,tmp+top+1,cmp1);sta1[++top1]=tmp[1];tree[rt].l1=top1;
        for (int i=2;i<=top;i++)
            if (tmp[i].x!=tmp[i-1].x)
            {
                while (top1>tree[rt].l1&&Dir(sta1[top1-1],sta1[top1],tmp[i])<=0)    top1--;
                sta1[++top1]=tmp[i];
            }
        tree[rt].r1=top1;sort(tmp+1,tmp+top+1,cmp2);sta2[++top2]=tmp[1];tree[rt].l2=top2;
        for (int i=2;i<=top;i++)
            if (tmp[i].x!=tmp[i-1].x)
            {
                while (top2>tree[rt].l2&&Dir(sta2[top2-1],sta2[top2],tmp[i])>=0)    top2--;
                sta2[++top2]=tmp[i];
            }
        tree[rt].r2=top2;
    }
    if (L==R)   return;
    if (cnt<=mid)   insert(ln);else insert(rn);
}
LL ask1(int l,int r)
{
    LL ret=-MAXLL;
    while (l<=r)
    {
        int len=(r-l)/3,mid1=l+len,mid2=r-len;LL dot1=p*sta1[mid1],dot2=p*sta1[mid2];
        if (dot1>dot2)  r=mid2-1,ret=max(ret,dot1);
        else    l=mid1+1,ret=max(ret,dot2);
    }
    return ret;
}
LL ask2(int l,int r)
{
    LL ret=-MAXLL;
    while (l<=r)
    {
        int len=(r-l)/3,mid1=l+len,mid2=r-len;LL dot1=p*sta2[mid1],dot2=p*sta2[mid2];
        if (dot1>dot2)  r=mid2-1,ret=max(ret,dot1);
        else    l=mid1+1,ret=max(ret,dot2);
    }
    return ret;
}
LL query(int rt,int l,int r)
{
    int L=tree[rt].l,R=tree[rt].r,mid=(L+R)>>1;LL ret=-MAXLL;
    if (l<=L&&r>=R)
    {
        if (R-L<MAXSIZE)    for (int i=L;i<=R;i++)  ret=max(ret,p*s[i]);
        else    ret=max(ret,p.y>0?ask1(tree[rt].l1,tree[rt].r1):ask2(tree[rt].l2,tree[rt].r2));
        return ret;
    }
    if (r<=mid) return query(ln,l,r); 
    if (l>mid)  return query(rn,l,r);
    return max(query(ln,l,mid),query(rn,mid+1,r));
}
void Getch(char &x) {x=getchar();while (x==' '||x=='\n'||x=='\r')   x=getchar();}
int calc(int x) {if (Flag!='E') x^=(lastans&MAXINT);return x;}
int main()
{
    in(n);Getch(Flag);build();int l,r;
    while (n--)
    {
        Getch(ch);
        if (ch=='A')    in(s[++cnt].x),in(s[cnt].y),s[cnt].x=calc(s[cnt].x),s[cnt].y=calc(s[cnt].y),insert(1);
        else    in(p.x),in(p.y),in(l),in(r),p.x=calc(p.x),p.y=calc(p.y),l=calc(l),r=calc(r),lastans=query(1,l,r),printf("%lld\n",lastans);
    }
}