BZOJ 4140 凸包+二进制分组

思路：

$(x_0-x)^2+(y_0-y)^2<=x^2+y^2$

$y>=(-x_0/y_0)x+(x_0^2+y_0^2)/2y0$

这显然就是凸包了

以一个斜率不断向下（上）走  找到第一个接触到的点

离线可以用cdq分治一发 他左边的判一判是不是都符合在圆内

这不能离线就很难搞了

感觉自己代码能力捉鸡   写棵splay估计要调死

强烈安利二进制分组做法

非常强

在log个凸包上 二分一下斜率 判一判 就好了

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=;
const double pi=acos(-),eps=1e-;
int n,op,flg,tot,num,top,rec,size[];
struct Point{double x,y;Point(double X=,double Y=):x(X),y(Y){}}jy,lst[N],q[N];
Point operator-(Point a,Point b){return Point(a.x-b.x,a.y-b.y);}
double operator*(Point a,Point b){return a.x*b.y-a.y*b.x;}
bool operator<(Point a,Point b){return a.x==b.x?a.y<b.y:a.x<b.x;}
bool cmp(Point a,Point b){return (a-lst[])*(b-lst[])>-eps;}
vector<Point>vec[];vector<double>ang[];
void Tubao(){
    for(int i=;i<=num;i++)if(lst[i]<lst[])swap(lst[],lst[i]);
    sort(lst+,lst++num,cmp),top=;
    for(int i=;i<=num;i++){
        while(top>&&(lst[i]-q[top])*(q[top]-q[top-])>-eps)top--;
        q[++top]=lst[i];
    }
}
double A(double x){return x<=-pi/?x+*pi:x;}
void insert(Point p){
    vec[++tot].push_back(p),size[tot]=,lst[num=]=p;
    while(tot>&&size[tot]==size[tot-]){
        for(int i=;i<vec[tot-].size();i++)lst[++num]=vec[tot-][i];
        size[tot-]+=size[tot],size[tot]=,vec[tot].clear(),ang[tot].clear(),tot--;
    }
    Tubao(),vec[tot].clear(),ang[tot].clear();
    if(top>){
        for(int i=;i<top;i++)vec[tot].push_back(q[i]),ang[tot].push_back(A(atan2(q[i+].y-q[i].y,q[i+].x-q[i].x)));
        vec[tot].push_back(q[top]),ang[tot].push_back(A(atan2(q[].y-q[top].y,q[].x-q[top].x)));
    }
    else vec[tot].push_back(q[]);
}
bool query(){
    double angle;
    angle=jy.y>?A(atan2(-jy.x,jy.y)):A(atan2(jy.x,-jy.y));
    for(int i=;i<=tot;i++){
        op=vec[i].size()==?:lower_bound(ang[i].begin(),ang[i].end(),angle)-ang[i].begin();
        if(*jy.x*vec[i][op].x+*jy.y*vec[i][op].y-jy.x*jy.x-jy.y*jy.y<eps)return ;
    }return ;
}
int main(){
    scanf("%d",&n);
    for(int i=;i<=n;i++){
        scanf("%d%lf%lf",&op,&jy.x,&jy.y),jy.x+=rec,jy.y+=rec;
        if(op==)flg=,insert(jy);
        else flg&&query()?(puts("Yes"),rec++):puts("No");
    }
}