求最多散点覆盖圆（求已知R的圆最多能覆盖的散点）

问题描述： 给n个点，给出圆的r，求在给定r下，圆最多能覆盖多少点？

O(n³)： 显然肯定有2个点在最多覆盖圆的边界上。可以枚举每2个点，求出以该两点连线为弦的圆，再枚举每个点是否在圆内。

O(n²logn)： 类似扫描线的做法，以每一个点为圆心化圆，枚举与其相交得圆，保存交点和角度，按角度排序后，扫一遍。
参考：[POJ 1981] Circle and Points(单位圆覆盖最多的点）

例题：POJ 1981 Circle and Points

//O(n2logn)做法
int GetCircleMostCoverNum(const Point *p, int n, double R)//p[0,n)
{
    //只测试了R=1的情况,即单位圆
    struct Node {double ang; int in;} arc[n*n];
    auto cmp=[&](const Node&a,const Node&b)->bool{
        if(dcmp(a.ang-b.ang)==0)
            return a.in>b.in;
        return dcmp(a.ang-b.ang)<0;
    };

    int cnt = 0;
    int ans = 0,maxans = 0;
    for(int i=0; i<n; i++){
        ans = 0; cnt = 0;
        for(int j=0; j<n; j++){
            if(Length(p[i]-p[j])>R*2) continue;
            double ang1 = atan2(p[i].y-p[j].y,p[i].x-p[j].x);
            double ang2 = acos(Length(p[i]-p[j])/(R*2));
            arc[cnt].ang = ang1-ang2;
            arc[cnt].in = 1; 
            cnt++;
            arc[cnt].ang = ang1+ang2;
            arc[cnt].in = -1; 
            cnt++;
        }
        sort(arc,arc+cnt,cmp);

        for(int j=0; j<cnt; j++) {
            ans+=arc[j].in;
            if(ans>maxans){
                maxans=ans;
            }
        }
    }
    return maxans;
}