poj1328 Radar Installation 区间贪心

题目大意：

　　在X轴选择尽量少的点作为圆心，作半径为d的圆。使得这些圆能覆盖所有的点。

思路：

　　把每个点都转化到X轴上。也就是可以覆盖这个点的圆心的位置的范围[a,b]。然后按照每个点对应的a从小到大排序。第一点需要特殊处理，我们赋值r=b₀ 。也就是使得第一个圆的圆心的横坐标尽量大。然后遍历剩下的点。对于i点，如果该点的a_i大于r, 就需要增加一个圆，圆心为b_i ；否则的话，把r更新为r与b_i中小的值。

代码：

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 #include<algorithm>
 6 using namespace std;
 7
 8 const int N = 1010;
 9 struct node
10 {
11     double x,y;
12     double a,b;
13 }p[N];
14 bool cmp(node n1,node n2)
15 {
16     return n1.a<n2.a;
17 }
18 int main()
19 {
20     //freopen("test.txt","r",stdin);
21     int n,d,ca=1;
22     while(scanf("%d%d",&n,&d)!=EOF&&n)
23     {
24         int i;
25         bool flag=0;
26         for(i=0;i<n;i++)
27         {
28             scanf("%lf%lf",&p[i].x,&p[i].y);
29             if(p[i].y>d) flag=1;
30         }
31         printf("Case %d: ",ca++);
32         if(flag) printf("-1\n");
33         else
34         {
35             for(i=0;i<n;i++)
36             {
37                 double x=p[i].x,y=p[i].y;
38                 p[i].a=-sqrt(d*d-y*y)+x;
39                 p[i].b=sqrt(d*d-y*y)+x;
40             }
41             sort(p,p+n,cmp);
42             double r=p[0].b;
43             int ans=1;
44             for(i=1;i<n;i++)
45             {
46                 if(p[i].a>r){
47                     r=p[i].b;
48                     ans++;
49                 }
50                 if(p[i].b<r) r=p[i].b;
51             }
52             printf("%d\n",ans);
53         }
54     }
55     return 0;
56 }