贪心 POJ 1328 Radar Installation

 

题目地址：http://poj.org/problem?id=1328

 1 /*  2  贪心  3  (转载)题意：有一条海岸线，在海岸线上方是大海，海中有一些岛屿，  4  这些岛的位置已知，海岸线上有雷达，雷达的覆盖半径知道，  5  问最少需要多少个雷达覆盖所有的岛屿。  6  (错误)思路：我开始是想从最左边的点雷达能探测的到的最右的位置出发，判断右边其余的点是否与该点距离小于d  7  是，岛屿数-1；不是，雷达数+1，继续。。。  8  (正确)思路：每个岛屿的座标已知，以雷达半径为半径画圆，与x轴有两个交点。  9  也就是说，若要覆盖该岛，雷达的位置范围是这两个交点。因此转化为覆盖区间的问题。  10  参考代码：http://www.cnblogs.com/kuangbin/archive/2011/07/30/2121838.html  11 */  12 #include <cstdio>  13 #include <iostream>  14 #include <algorithm>  15 #include <cstring>  16 #include <string>  17 #include <cmath>  18 using namespace std;  19  20 const int MAXN = 1e3 + 10;  21 const int INF = 0x3f3f3f3f;  22 struct NODE  23 {  24 int x, y;  25 double l, r;  26 }node[MAXN];  27  28 bool cmp(NODE a, NODE b)  29 {  30 return a.l < b.l;  31 }  32  33 bool ok(int n, int d)  34 {  35 if (d < 0) return false;  36 for (int i=1; i<=n; ++i)  37  {  38 if (node[i].y > d) return false;  39  }  40  41 return true;  42 }  43  44 void work(int n, int d)  45 {  46 int cnt = 1;  47 double now = node[1].r;  48 for (int i=2; i<=n; ++i)  49  {  50 if (now > node[i].r) now = node[i].r;  51 if (now < node[i].l)  52  {  53 now = node[i].r;  54 cnt++;  55  }  56  }  57 printf ("%d\n", cnt);  58 }  59  60 int main(void) //POJ 1328 Radar Installation  61 {  62 //freopen ("I.in", "r", stdin);  63  64 int n, d;  65 int cnt = 0;  66 while (~scanf ("%d%d", &n, &d) && n && d)  67  {  68 for (int i=1; i<=n; ++i)  69  {  70 scanf ("%d%d", &node[i].x, &node[i].y);  71 node[i].l = (double)node[i].x - sqrt ((double)d * d - node[i].y * node[i].y);  72 node[i].r = (double)node[i].x + sqrt ((double)d * d - node[i].y * node[i].y);  73  }  74 sort (node+1, node+1+n, cmp);  75 printf ("Case %d: ", ++cnt);  76 if (!ok (n, d))  77  {  78 printf ("%d\n", -1); continue;  79  }  80  work (n, d);  81  }  82  83 return 0;  84 }  85  86 /*  87 void work(int n, int d)  88 {  89  int i = 1;  90  int j = 1;  91  int next = 1;  92  int num = 0;  93  int cnt = 0;  94  while (num < n)  95  {  96  double res = node[i].x + sqrt (d * d - node[i].y * node[i].y);  97  cnt++; num++;  98  int flag = 0;  99  fors (j=i+1; j<=n; ++j) 100  { 101  if (pow (node[j].x - res, 2) + pow (node[j].y, 2) <= d * d) 102  { 103  num++; next = j; flag = 1; 104  } 105  } 106  if (flag) i = next + 1; 107  else i++; 108  } 109  printf ("%d\n", cnt); 110 } 111 */