题意:
维护一个向量集合,在线支持以下操作:
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<[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<by2−y1x2−x1
当b>0
y2−y1x2−x1>−ab
当b<0
y2−y1x2−x1<−ab
当b=0时找出y的区间最大最小值比较一下就好。。
观察
y2−y1x2−x1
这是两点间的斜率,就是说两点谁比较优是与他们间的斜率有关的
按照套路,我们维护一个斜率单调的凸壳就可以二分了。。
显然,b>0时要维护斜率单调递减,b<0时要维护斜率单调递增。
注意只会在尾部加入,没有中间的修改,让我想起了以前做(乱搞)过的一道题。我是O(n),正解是高逼格的二进制分组线段树。
那我们也用二进制分组线段树,线段树节点维护这个区间的两个凸壳。每个点只会往父亲合并一次,合并是O(n log n)的。合并时遇到两个x坐标相同,即没有斜率的点对时手动讨论一发就知道谁比较优了。而查询则要在log个凸壳上二分,是O(
nlog2n
)的。
二进制分组就是这个线段树节点包含区间已经塞满后再向父亲合并。