BZOJ 3533 Sdoi2014 向量集

题目大意：给定一个平面，维护下列操作：

1.插入一个向量

2.询问第l到r个插入的向量中与某个向量点积的最大值

强制在线

Sol:

向量点积????欺负我没学向量喵(虽然的确没学，看来要去找数学老师了QAQ)

题目给了点积公式(x,y)^(a,b)=xa+yb

恩，，，推一波式子

按给定向量的时间顺序建线段树，每个节点维护凸包

二分的话要二分斜率，直接乘的话害怕爆long long又不想用double，就写三分好了反正复杂度差不多

但是。。。怎么动态维护凸包呢？没做一次合并是O(size)的，如果插入一次所影响的log个节点都重建一遍的话单次就O(nlogn)了

如果离线的话可以一次全部建好，复杂度是O(nlogn)的

在线的话，我们发现，因为是顺序插入，如果有一个节点没有被填满的话就不会被用到

这提示我们可以当一个点被填满后再建凸包，复杂度...看起来不会比离线多吧

然后对于一个查询，划分成了log个区间，我们需要把这些区间合并成一个凸包吗？显然不，每个区间的凸包上三分一下最后去Max就行啦！

Code:

#include<bits/stdc++.h>
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define decode(x) ((x)^(ans&0x7fffffff))
typedef long long ll;
using namespace std;
const int maxn = 400009;
struct Point
{
	ll x,y;
	Point(ll _x,ll _y):x(_x),y(_y){}
	Point operator + (const Point &X){return Point(x+X.x,y+X.y);}
	Point operator - (const Point &X){return Point(x-X.x,y-X.y);}
	long long operator ^ (const Point &X){return x*X.x+y*X.y;}
	long long operator * (const Point &X){return x*X.y-y*X.x;}
	bool operator < (const Point &X)
	{
		if(x==X.x) return y<X.y;
		return x<X.x;
	}
};
typedef vector<Point> Convex;
struct sg_tree{Convex a,b;}node[maxn<<2];
int n,tot,flag;
ll ans;

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

ll calc(Convex &Hull,Point x)
{
	int l=0,r=Hull.size()-1;
	ll res=-1e17;
	while(r-l>2)
	{
		int lmid=(l+l+r)/3,rmid=(l+r+r)/3;
		if((Hull[lmid]^x)>(Hull[rmid]^x)) r=rmid;
		else l=lmid;
	}
	for(int i=l;i<=r;i++) res=max(res,(Hull[i]^x));
	return res;
}

void Cover(Convex &tar,Convex &l,Convex &r)
{
	int t1=l.size(),t2=r.size(),p1=0,p2=0,top=0;
	while(p1<t1||p2<t2)
	{
		Point now(0,0);
		if(p1==t1) now=r[p2++];else if(p2==t2) now=l[p1++];
		else if(l[p1]<r[p2]) now=l[p1++];else now=r[p2++];
		while(top>=2&&(tar[top-1]-tar[top-2])*(now-tar[top-1])>=0) tar.pop_back(),top--;
		top++;tar.push_back(now);
	}
}
void Below(Convex &tar,Convex &l,Convex &r)
{
	int t1=l.size(),t2=r.size(),p1=0,p2=0,top=0;
	while(p1<t1||p2<t2)
	{
		Point now(0,0);
		if(p1==t1) now=r[p2++];else if(p2==t2) now=l[p1++];
		else if(l[p1]<r[p2]) now=l[p1++];else now=r[p2++];
		while(top>=2&&(tar[top-1]-tar[top-2])*(now-tar[top-1])<=0) tar.pop_back(),top--;
		top++;tar.push_back(now);
	}
}

void Build(int rt)
{
	Cover(node[rt].a,node[rt<<1].a,node[rt<<1|1].a);
	Below(node[rt].b,node[rt<<1].b,node[rt<<1|1].b);
}
void insert(int l,int r,int rt,int pos,Point x)
{
	if(l==r)
	{
		node[rt].a.push_back(x);
		node[rt].b.push_back(x);
		return ;
	}int mid=(l+r)>>1;
	if(pos<=mid) insert(lson,pos,x);
	else insert(rson,pos,x);
	if(pos==r) Build(rt);
}
ll query(int l,int r,int rt,int left,int right,Point x)
{
	if(l==left&&right==r)
	{
		if(x.y<0) return calc(node[rt].b,x);
		else return calc(node[rt].a,x);
	}int mid=(l+r)>>1;
	if(right<=mid) return query(lson,left,right,x);
	else if(left>mid) return query(rson,left,right,x);
	else return max(query(lson,left,mid,x),query(rson,mid+1,right,x));
}

int main()
{
	char opt[5];
	n=read();scanf("%s",opt+1);
	if(opt[1]!='E') flag=1;
	for(int i=1;i<=n;i++)
	{
		scanf("%s",opt+1);
		int x=read(),y=read();
		if(flag) x=decode(x),y=decode(y);
		if(opt[1]=='A')
		{
			tot++;
			insert(1,n,1,tot,Point(x,y));
		}
		else
		{
			int l=read(),r=read();
			if(flag) l=decode(l),r=decode(r);
			printf("%lld\n",ans=query(1,n,1,l,r,Point(x,y)));
		}
	}
	return 0;
}