[BZOJ2961]共点圆-[凸包+cdq分治]

Description

传送门

Solution

考虑对于每一个点：

设圆的坐标为(x,y)，点的坐标为(x0,y0)。依题意得，当一个点在圆里，需要满足(x-x0)²+(y-y0)²<=x²+y²。

化简得x0²+y0²<=2x0*x+2y0*y。

当y0>0，x*(-x0/y0)+0.5y0+x0²/(2*y0)<=y，这是一个半平面的式子；当y0<0时同理，但是要变号。

所以对于某个点(x0,y0)，我们构造出在它前面所有圆心的凸包。凸包应分为上下。

通过以上式子我们可以得出，当y0>0时应在下凸包上找点(x,y)【该点为直线y=-x0/y0与下凸包的切点，即若此点满足要求，其他任何点都会满足要求。通过这个条件，我们也可以理解把该点理解为2x0*x+2y0*y最小的点】，反之则应该在上凸包上找。

好的让我们假设目前的y0>0，由于凸包上的点(x,y)是按极角排序，x*x0+y*y0是单峰的(这个式子是向量的点乘，其几何意义为：向量a点乘向量b=向量a的长度*向量b的长度*cos(向量a，b的夹角)。所以根据这个定义，证明，就画图吧。qaq正儿八经公式太麻烦了。)上凸包也是一样的。

Code

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
const double eps=1e-;
int n;
bool ans[];
struct W{
    int tp;double x,y;
    friend double operator *(W a,W b){return a.x*b.y-a.y*b.x;}
    friend double operator ^(W a,W b){return a.x*b.x+a.y*b.y;}
    friend W operator -(W a,W b){return W{,a.x-b.x,a.y-b.y};}
    friend bool operator <(W a,W b){return (fabs(a.x-b.x)<eps)?a.y<b.y:a.x<b.x;}
}w[],t[],st[][];//st[0]-上凸包，st[1]-下凸包
double ask(W a,W b){return a.x*b.x+a.y*b.y;}
bool _ok[];
int top0,top1;
bool query(int id)
{
    int l,r,mid1,mid2;double minn=1e12;
    if (w[id].y<)
    {
        l=;r=top0;
        while (r-l>)
        {
            mid1=(l*+r)/;mid2=(r*+l)/;
            if ((w[id]^st[][mid1])<(w[id]^st[][mid2])) r=mid2;
            else l=mid1;
        }
        for (int i=l;i<=r;i++) minn=min(minn,w[id]^st[][i]);

    } else
    {
        l=;r=top1;
        while (r-l>)
        {
            mid1=(l*+r)/;mid2=(r*+l)/;
            if ((w[id]^st[][mid1])<(w[id]^st[][mid2])) r=mid2;
            else l=mid1;
        }
        for (int i=l;i<=r;i++) minn=min(minn,w[id]^st[][i]);
    }
    if (*minn-(w[id].x*w[id].x+w[id].y*w[id].y)<eps) ans[id]=;

}
void solve(int l,int r)
{
    if (l==r) return;
    int mid=(l+r)/,tot=;
    solve(l,mid);
    solve(mid+,r);
    for (int i=l;i<=mid;i++) if (!w[i].tp) t[++tot]=w[i];
    sort(t+,t+tot+);
    top1=,top0=;
    for (int i=;i<=tot;i++)
    {
        while (top0>&&(st[][top0]-st[][top0-])*(t[i]-st[][top0])>-eps) top0--;
        st[][++top0]=t[i];
    }
    st[][top0+].x=st[][top0].x;st[][top0+].y=-1e12;
    for (int i=tot;i;i--)
    {
        while (top1>&&(st[][top1]-st[][top1-])*(t[i]-st[][top1])>-eps) top1--;
        st[][++top1]=t[i];
    }
    st[][top1+].x=st[][top1].x;st[][top1+].y=1e12;
    for (int i=mid+;i<=r;i++) if (w[i].tp&&ans[i]) query(i);
}
bool _is=;
int main()
{
    scanf("%d",&n);
    for (int i=;i<=n;i++)
    {
        scanf("%d%lf%lf",&w[i].tp,&w[i].x,&w[i].y);
        if (!w[i].tp) _is=;else ans[i]=_is;
    }
    solve(,n);
    for (int i=;i<=n;i++) if (w[i].tp)
     if (ans[i]) printf("Yes\n");else printf("No\n");
}