题目大意:选取最少的放雷达的点,使得每一个点都可以被探测到。
思路:
1.错误的贪心算法
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <iostream>
using namespace std;
const int N = 1000+10;
int n,d;
struct node {
double l;
double r;
};
node ans[N];
double input[N][2];
int cnt = 0;
bool flag = true;
void calculate() {
for (int i=0 ; i<n; i++ ) {
if ( d<0 || d<input[i][1] ) {
flag = false;
break;
} else {
double dis = sqrt( d*d-input[i][1]*input[i][1] ) ;
ans[i].l = input[i][0]-dis;
ans[i].r = input[i][0]+dis;
}
}
}
bool cmp( node a,node b ) {
if ( a.l!=b.l ) return a.l<b.l;
else return a.r>b.r;
}
int main()
{
int cases = 0;
while ( cin>>n>>d ) {
++cases;
if (n==0&&d==0) break;
memset(ans,0,sizeof(ans));
flag = true;
for (int i=0 ;i<n; i++) {
double x,y;
cin>>x>>y;
input[i][0] = x;
input[i][1] = y;
}
calculate();
if ( !flag || d<0 ) {
cout<<"Case "<<cases<<": "<<-1<<endl;
} else {
sort( ans,ans+n,cmp );
int all = 0;
double e;
int index = 0 ;
while ( index<n ) {
e = ans[index].r;
++index;
while ( index<n && ans[index].l<=e )
++index;
all++;
}
printf("Case %d: %d\n",cases,all);
}
}
return 0;
}
在这种情况下,我依据左端点进行排序,相同情况下按照右端点递减顺序。然后以第一个区间的末尾来选择后面的区间,如果后面的区间的起点小于该末尾,
那我认为该两个区间有交集,但是我忽略了如下的情况。
,假设我现在的区间是A,我的算法选择的是A之后的每一个区间与A是否有交集,而不是A之后(包含A)的每一个区间都有交集 的情况。
正确做法
最小的点数覆盖的充分条件是每一个点都包含于尽可能多的区间。
以右端点从小到大排序(右端点相等,左端点的排序方式无所谓),来筛选左端点,这种情况可以保证每个区间都有交集。我们将第一区间的末尾作为标准点,因为右端点是增长的,
所以我们只需要和左端点进行比较即可。
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <iostream>
using namespace std;
const int N = 1000+10;
int n,d;
struct node {
double l;
double r;
};
node ans[N];
double input[N][2];
int cnt = 0;
bool flag = true;
void calculate() {
for (int i=0 ; i<n; i++ ) {
if ( d<0 || d<input[i][1] ) {
flag = false;
break;
} else {
double dis = sqrt( d*d-input[i][1]*input[i][1] ) ;
ans[i].l = input[i][0]-dis;
ans[i].r = input[i][0]+dis;
}
}
}
bool cmp( node a,node b ) {
if ( a.r!=b.r ) return a.r<b.r;
else return a.l>b.l;
}
int main()
{
int cases = 0;
while ( cin>>n>>d ) {
++cases;
if (n==0&&d==0) break;
memset(ans,0,sizeof(ans));
flag = true;
for (int i=0 ;i<n; i++) {
double x,y;
cin>>x>>y;
input[i][0] = x;
input[i][1] = y;
}
calculate();
if ( !flag || d<0 ) {
cout<<"Case "<<cases<<": "<<-1<<endl;
} else {
sort( ans,ans+n,cmp );
int all = 0;
double e;
int index = 0 ;
while ( index<n ) {
e = ans[index].r;
++index;
while ( index<n && ans[index].l<=e )
++index;
all++;
}
printf("Case %d: %d\n",cases,all);
}
}
return 0;
}