POJ1328Radar Installation（贪心）

对于每一个点，可以找到他在x轴上的可行区域，这样的话就变为了对区间的贪心。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<map>
#include<vector>
#include<set>
#include<stack>
#include<queue>
#include<algorithm>
#include<cmath>
#include<stdlib.h>
using namespace std;
#define MAX(a,b) (a > b ? a : b)
#define MIN(a,b) (a < b ? a : b)
#define MAXN  10000001
#define INF 1000000007
#define mem(a) memset(a,0,sizeof(a))
#define eps 1e-15

struct node{double s,t;}ma[1001];
double R;
int N;


const int island_max=1000;


void get_len(int index,double a,double b)
{
    double temp = sqrt(R*R - b*b);
    ma[index].s=a-temp;
    ma[index].t=a+temp;
}

int cmp(node a,node b)
{
    if(a.t!=b.t)return (a.t < b.t);
    return (b.s > a.s);
}

int main()
{
    int ca = 1;
    while(~scanf("%d%lf",&N,&R) && (N || R))
    {
        int i;
        double a,b;
        int flag = 0;
        for(i=0;i<N;i++)
        {
            scanf("%lf %lf",&a,&b);
            if(b<=R && !flag)
            {
                get_len(i,a,b);
            }
            else flag = 1;
        }
        if(flag){printf("Case %d: -1
",ca++);continue;}
        sort(ma,ma+N,cmp);
        double key = ma[0].t;
        int ans=1;
        for(i=1;i<N;i++)
        {
            if(ma[i].s-key >eps)
            {
                key = ma[i].t;
                ans ++;
            }
        }
         printf("Case %d: %d
",ca++,ans);
    }
    return 0;
}