POJ 1328 Radar Installation

小小的一道贪心算法题，

由于读数据的时候未读完数据居然卡了两个小时。

不过这次代码写的自我感觉不错，结构也清楚，这种小题现在对我来讲毫无压力啊。。。。

#include<iostream>
#include<stdlib.h>
#include<math.h>

using namespace std;

struct ln
{
    double x1;
    double x2;
}section[1010];

int cmp(const void *a,const void *b)
{
    struct ln* c=(ln *)a;
    struct ln* d=(ln *)b;
    return c->x1>d->x1?1:-1;
    //return (*(ln*)a)->x1>(*(ln*)b)->x1?1:-1;//注意解引用操作符。。。
}

int main()
{
    int n,d;
    int cases=0;
    while(cin>>n>>d&&!(n==0&&d==0))
    {
        cases++;
        bool flag=false;
        int solution=1;
        for(int i=0;i<n;i++)
        {
            int x,y;
            cin>>x>>y;
            if(y>d||y<(-1)*d||d<0)
            {
                flag=true;
                solution=-1;
         //     break; //就是这条万恶的语句，导致runtime error。。。。
            }
            section[i].x1=(double)x-sqrt((double)(d*d-y*y));
            section[i].x2=(double)x+sqrt((double)(d*d-y*y));
        }
        if(!flag)
        {
            qsort(section,n,sizeof(section[0]),cmp);
            for(int i=0;i<n-1;i++)
            {
                if(section[i].x2>=section[i+1].x1)
                {
                    if(section[i+1].x2>section[i].x2)
                        section[i+1].x2=section[i].x2;

                }
                else    solution++;
            }
        }
        cout<<"Case "<<cases<<": "<<solution<<endl;
    }
    return 0;
}