HDU 1883 Phone Cell （圆覆盖最多点）

题目链接

题意 ： 给你很多点和一个半径r，这个半径为r的圆能覆盖的最多的点是多少。

思路 ： 对每个点做半径为 r 的圆, 求交集，交集最多的区域的被覆盖次数就是能覆盖的最多的点。贴两个链接，分析的挺好，代码写得挺好

#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[4030] ;

inline 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 Sqrt(double x)
{
    return x*x ;
}
double dist(const point &a,const point &b)
{
    return sqrt(Sqrt(a.x-b.x)+Sqrt(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 ;
        for(int i = 0 ; i < n ; i++)
        {
            int 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) ;
                if(s < 0) s += 2*PI ;//角度区间修正
                double ph = acos(d/2.0/r) ;//圆心角转区间
                q[m++].angle = s - ph + 2*PI ;q[m-1].flag = 1 ;
                q[m++].angle = s + ph + 2*PI ;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) ;
        }
        printf("It is possible to cover %d points.
",ans+1) ;
    }
return 0 ;
}