poj 1328 贪心算法和快排（快排细节）

把海岸线看做直线，海上有小岛，岸上有雷达站
雷达站只能覆盖半径为d的区域
把地图看做直角坐标系，海岸线为X轴，陆地在下面
给定小岛坐标和雷达站的覆盖半径。
任务是建立最小的雷达站以覆盖所有的小岛
多测试
n d
n line
换行开始下个test
-1表示无解

思路：把每一个小岛当做半径为d的圆的圆心，算出圆与x轴的交点。记录左右两个坐标，进行快排后，区间重叠的个数就是雷达站个数。思路很简单，但是我在这题奉献了不少WA，原因是没弄清qsort的返回值，qsort的返回值是整形，所以在对double型数据排序的时候要注意返回值的问题，所以qosrt得cmp就得改变为：

int cmp(const void *a,const void *b)
{
return ((double *)a)[0]>((double *)b)[0] ?1:-1;
}

也就是加了个 ?1:-1;这样就能限定返回值是整数了。这样就不会WA了

贴上用数组实现和用结构实现的AC代码和测试数据，同时我也会更新代码集中的代码。

#include<stdio.h>
#include<string.h>
#include <iostream>
#include<queue>
#include<math.h>
#include<stdlib.h>
using namespace std;
#define maxn 1005
int cmp(const void *a,const void *b)
{
return ((double *)a)[0]>((double *)b)[0] ?1:-1;
}
double dir[maxn][2];
int main()
{
int n,ans;
double d;
int ca=1;
while(scanf("%d %lf",&n,&d)!=EOF)
{
if(n==0&&d==0)
return 0;
int jg=0;
for(int i=0; i<n; i++)
{
double tmp1,tmp2;
scanf("%lf %lf",&tmp1,&tmp2);
if(tmp2>d)
jg=1;
double tmp3=sqrt(d*d-tmp2*tmp2);
dir[i][0]=tmp1-tmp3;
dir[i][1]=tmp1+tmp3;
}
if(!jg)
{
qsort(dir,n,sizeof(dir[0]),cmp);
/*for(int i=0;i<n;i++)
printf("%lf\n",dir[i][0]);*/
double l=dir[0][0];
double r=dir[0][1];
ans=1;
for(int i=1; i<n; i++)
{
if(dir[i][0]>r)
{
ans++;
l=dir[i][0];
r=dir[i][1];
}
else if(dir[i][1]<r)
{
l=dir[i][0];
r=dir[i][1];
}
}
}
else
ans=-1;
printf("Case %d: %d\n",ca,ans);
ca++;
}
}

#include<stdio.h>
#include<string.h>
#include <algorithm>
#include <iostream>
#include<queue>
#include<math.h>
#include<stdlib.h>
using namespace std;
#define maxn 1005
typedef struct Node{double xr;double xl;} ode;
ode node[maxn];
int cmp(const void *a,const void *b)
{
return ((ode *)a)->xl>((ode *)b)->xl? 1:-1;
}
int main()
{
int n,ans;
double d;
int ca=1;
while(scanf("%d %lf",&n,&d)!=EOF)
{
if(n==0&&d==0)
return 0;
int jg=0;
for(int i=0; i<n; i++)
{
double tmp1,tmp2;
scanf("%lf %lf",&tmp1,&tmp2);
if(fabs(tmp2)>d)
jg=1;
node[i].xl=tmp1*1.0-sqrt(d*d-tmp2*tmp2);
node[i].xr=tmp1*1.0+sqrt(d*d-tmp2*tmp2);
}
if(!jg)
{
qsort(node,n,sizeof(node[0]),cmp);
//for(int i=0;i<n;i++)
//printf("%lf\n",node[i].xl);
ode tmp=node[0];
ans=1;
for(int i=1; i<n; i++)
{
if(node[i].xl>tmp.xr)
{
ans++;
tmp=node[i];
}
else if(node[i].xr<tmp.xr)
{
tmp=node[i];
}
}
}
else
ans=-1;
printf("Case %d: %d\n",ca,ans);
ca++;
}
}

大神的测试数据：

测试数据，来源于poj讨论组

2 5-3 4-6 34 5-5 3-3 52 33 320 8-20 7-18 6-5 8-21 8-15 7-17 5-1 5-2 3-9 61 22 33 44 55 66 77 88 79 610 50 02 30 22 32 30 21 33 31 2-3 22 48 52 4-4 4-3 3-3 1-3 0-1 00 56 03 01 2-3 12 13 21 2-3
12 11 20 22 30 22 34 -54 34 32 36 -93 -31 2-3 22 16 21 21 21 2-3 12 10 01 20 22 30 21 33 101 102 34 53 51 102 34 54 71 102 34 50 03 91 102 34 52 50 38 30 0

运行结果：

Case 1: 1

Case 2: 2

Case 3: 4

Case 4: 1

Case 5: 1

Case 6: -1

Case 7: 3

Case 8: -1

Case 9: 2

Case 10: 1

Case 11: 1

Case 12: -1

Case 13: -1

Case 14: 2

Case 15: 1

Case 16: 1

Case 17: 1

Case 18: -1

Case 19: -1

Case 20: -1

Case 21: 1