[二进制分组维护凸包]BZOJ 4140—— 共点圆加强版

题目描述

在平面直角坐标系中，Wayne需要你完成n次操作，操作只有两种：

1.0 x y。表示在坐标系中加入一个以(x, y)为圆心且过原点的圆。

2.1 x y。表示询问点(x, y)是否在所有已加入的圆的内部（含圆周），且至少在一个圆内部（含圆周）。

为了减少你的工作量，题目保证圆心严格在x轴上方（纵坐标为正），且横坐标非零。

解题思路

如果不强制在线CDQ分治也是可以的。

现在强制在线了，自然想到二进制分组。

题意即满足$(x_i-x)^2+(y_i-y)^2<=x_i^2+y_i^2$

化简成$y_i>=(x^2+y^2)/(2*y)-(x/y)*x_i$就可以凸包了

即每个点都在该直线的上方，维护一个下凸壳就可以了。

查询答案时用二分斜率逼近就好了。

#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;
}