HDU 1883 Phone Cell

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <iostream>

#include <algorithm>
using namespace std ;
const double eps = 1e-8;
const double PI = acos(-1.0) ;
int n ;
double r ;
struct point
{
    double x,y ;

} p[2010] ;
struct node
{
    double angle ;
    int flag ;
} q[5000] ;
int dcmp(double d)
{
    return d < -eps ? -1 : d > eps ;
}
bool cmp(const node &a,const node &b)///角度区间排序,区间上限永远在上;
{
    if(dcmp(a.angle-b.angle) == 0 ) return a.flag > b.flag ;

    return a.angle < b.angle ;
}
double Cal(double x)
{
    return x*x ;
}
double dist(const point &a,const point &b)
{
    return sqrt(Cal(a.x-b.x)+Cal(a.y-b.y)) ;
}

int main()
{
    while(~scanf("%d%lf",&n,&r))
    {
        if(n == 0) break ;
        for(int i = 0 ; i < n ; i++)
            scanf("%lf%lf",&p[i].x,&p[i].y) ;

        int ans = 0,m = 0;
        for(int i = 0 ; i < n ; i++)
        {
            m = 0 ;
            for(int j = 0 ; j < n ; j++)
            {
                if(i == j) continue ;
                double d = dist(p[i],p[j]) ;
                if(d > 2*r+0.001) continue ;
                double s = atan2(p[j].y-p[i].y,p[j].x-p[i].x) ;
                ///atan2 和 atan 的区别;
                ///1：参数的填写方式不同；
                ///2：atan2 的优点在于 如果 x2-x1等于0 依然可以计算，但是atan函数就会导致程序出错；
                if(s < 0) s += 2*PI ;///角度区间修正
                double ph = acos(d/2.0/r) ;///圆心角转区间
                q[m++].angle = s - ph + 2*PI ;///图中b-a,避免负值+2*PI;
                q[m-1].flag = 1 ;///标记为上界

                q[m++].angle = s + ph + 2*PI ;///图中a+b;
                q[m-1].flag = -1 ;///标记下界
            }
            sort(q,q+m,cmp) ;
            int sum = 0 ;
            for(int j = 0 ; j < m ; j++)///这种记录区间最大交集次数的方式很巧妙,当时我对这里还是挺困惑的,读者若不懂应稍加思考.
                ans = max(ans,sum += q[j].flag) ;
        }
        /*for(int j = 0; j < m; j++)
        {
            cout<<"angle = "<<q[j].angle<<endl;
            cout<<"flag = "<<q[j].flag<<endl;
        }*/
        printf("It is possible to cover %d points.\n",ans+1) ;
    }
    return 0 ;
}