题目描述
在平面直角坐标系中,Wayne需要你完成n次操作,操作只有两种:
1.0 x y。表示在坐标系中加入一个以(x, y)为圆心且过原点的圆。
2.1 x y。表示询问点(x, y)是否在所有已加入的圆的内部(含圆周),且至少在一个圆内部(含圆周)。
为了减少你的工作量,题目保证圆心严格在x轴上方(纵坐标为正),且横坐标非零。
解题思路
如果不强制在线CDQ分治也是可以的。
现在强制在线了,自然想到二进制分组。
题意即满足 (xi−x)2+(yi−y)2<=x2i+y2i
化简成 yi>=(x2+y2)/(2∗y)−(x/y)∗xi 就可以凸包了
即每个点都在该直线的上方,维护一个下凸壳就可以了。
查询答案时用二分斜率逼近就好了。
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=500005;
struct jz{
double x,y;
jz(double x=0,double y=0):x(x),y(y){}
jz operator-(const jz &b)const{return jz(x-b.x,y-b.y);}
bool operator<(const jz &b)const{
if (x==b.x) return y<b.y;
return x<b.x;
}
}a[maxn],q[maxn];
double cross(jz x,jz y){return x.x*y.y-x.y*y.x;}
int n,ans,L[maxn],R[maxn],top,num;
void build(double x,double y){
a[++num]=jz(x,y);L[++top]=num;
while(top>1&&num-L[top]+1==L[top]-L[top-1]) top--;
sort(a+L[top],a+1+num);
int r=L[top];q[r]=a[r];
for (int i=L[top]+1;i<=num;i++){
while(r>L[top]&&cross(q[r]-q[r-1],a[i]-q[r-1])<=0) r--;
q[++r]=a[i];
}
R[top]=r;
}
int find(int L,int R,double k){
int mid,ans,l=L;
while(L<=R){
int mid=L+(R-L>>1);
if (mid==l||k>(q[mid].y-q[mid-1].y)/(q[mid].x-q[mid-1].x)) L=(ans=mid)+1;else R=mid-1;
}
return ans;
}
bool query(double x,double y){
if (!top) return 0;
double k=-x/y,qt=(x*x+y*y)/(y*2);
for (int i=1;i<=top;i++){
int w=find(L[i],R[i],k);
if (q[w].x*k+qt>q[w].y) return 0;
}
return 1;
}
int main(){
freopen("exam.in","r",stdin);
freopen("exam.out","w",stdout);
scanf("%d",&n);
for (int i=1;i<=n;i++){
int h;double x,y;
scanf("%d%lf%lf",&h,&x,&y);
x+=ans,y+=ans;
if (h==0) build(x,y);
else if (query(x,y)) ans++,printf("Yes\n");else printf("No\n");
}
return 0;
}