poj 1328

这是一道经典的贪心题

题意：用一个平面坐标系来表示一个地面，其中x轴是海岸线，y<0的半平面表示地面，而y>0的半平面表示海洋，海洋上有n个岛屿，我们需要在海岸线（x轴）上建立一些雷达站（半径全部为d），使得将所有岛屿都包含在雷达站的范围内。

这题比较容易看出是贪心，对于一个给定的岛屿（x,y），容易知道只有在坐标范围是 l = （x-sqrt（d*d-y*y），0）到

r = （x+sqrt（d*d-y*y））建立雷达站时，才能将岛屿包含进来，所以对于每个给定岛屿，我们记录它的这两个值，然后先根据 r 进行排序，如果 r 值相同，则根据 l 值进行排序，我们假设前k个岛屿已经被修好的雷达站给包含进去了，且第k+1个岛屿还没被包含，于是我们要在第k+1个岛屿的 l 到 r 之间修一个雷达站，但这个雷达站修在哪里好呢？思考一下我们知道，这跟后面岛屿的坐标有关，假设第k+2个岛屿的值为 l1 和 r1 ，如果 l1 小于 r ，说明两个岛屿之间是有交界的，如果 r1<r 那么雷达站修在 l1 到 r1 之间就能将两个雷达站包含进来，如果 r1>r ，那么就将雷达站修在 l1 到 r 之间，接下来再考虑第k+3个岛屿。如此下去，知道所有岛屿都被包含进来

代码如下：

#include<iostream> #include<algorithm> #include<cmath> using namespace std; struct point { double l,r; }p[1010]; bool operator <(const point a,const point b) { if(a.l==b.l) return a.r<b.r; return a.l<b.l; } int main() { int n,m,i,ans,Case=0; double x,y,rec,flag; while(scanf("%d %d",&n,&m),n||m) { Case++; flag=ans=1; for(i=1;i<=n;i++) { scanf("%lf %lf",&x,&y); if(fabs(y)>m) flag=0; p[i].l=x-sqrt(m*m-y*y); p[i].r=x+sqrt(m*m-y*y); } sort(p+1,p+1+n); rec=p[1].r; for(i=2;i<=n;i++) { if(p[i].l>rec) { ans++; rec=p[i].r; } else { if(p[i].r<rec) rec=p[i].r; } } if(flag) printf("Case %d: %d/n",Case,ans); else printf("Case %d: -1/n",Case); } return 0; }