B - Radar Installation poj 1328【贪心】

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<algorithm>
using namespace std;
#define N 1100
struct node {
    double x1,x2;
}c[N];

double cmp(struct node a,struct node b)
{
    if(a.x1 > b.x1 )
        return a.x1 < b.x1 ;
}
int main()
{
    int i,j,n,R,ans,t=0;
    double x,y,coords;
    while(scanf("%d%d",&n,&R),n!=0&&R!=0)//不能用(n+R)作为循环条件，因为3 -3这种情况不满足 
    {
        ans = 1;
        if(R <= 0)//先判断半径是否满足条件 
            ans = -1;

        for(i = 1; i <= n; i ++)
        {
            scanf("%lf%lf",&x,&y);//坐标是实数 
            if(fabs(y) > R)//纵坐标大于半径时  不满足条件 
            {
                ans = -1;
            }
            else
            {
                coords = sqrt(R*R-y*y);
                c[i].x1 = x - coords;//计算与x轴相交的左右端点 
                c[i].x2 = x + coords;
            }
        }
        if( ans == -1)
        {
            printf("Case %d: %d\n",++t,ans);
            continue;
        }
        sort(c+1,c+1+n,cmp);//结构体排序 
        coords = c[1].x2 ;//初始化为最左坐标的右端点 
        
        for(i = 2; i <= n; i ++)
        {
            if(c[i].x1 > coords)//如果下一个坐标的左端点大于上一个坐标右端点 说明两者没有公共区间 
            {
                ans ++;
                coords = c[i].x2 ;
            }
            else if(c[i].x2 < coords)//如果下一个坐标的右端点小于上一个坐标的右端点，说明存在公共区间 
                coords = c[i].x2 ;
        }
        printf("Case %d: %d\n",++t,ans);
    }
    return 0;
}