POJ 1328 区间点选择 贪心

题目链接

题目大意：选取最少的放雷达的点，使得每一个点都可以被探测到。

思路：

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;

}