一个小贪心 区间覆盖问题

一个小贪心，东西很简单，但毕竟是自己做的：

喷水装置（二）

时间限制：3000 ms  |  内存限制：65535 KB

难度：4

描述

　　有一块草坪，横向长w,纵向长为h,在它的橫向中心线上不同位置处装有n(n<=10000)个点状的喷水装置，每个喷水装置i喷水的效果是让以它为中心半径为Ri的圆都被润湿。请在给出的喷水装置中选择尽量少的喷水装置，把整个草坪全部润湿。

输入

　　第一行输入一个正整数N表示共有n次测试数据。

　　每一组测试数据的第一行有三个整数n,w,h，n表示共有n个喷水装置，w表示草坪的横向长度，h表示草坪的纵向长度。

　　随后的n行，都有两个整数xi和ri,xi表示第i个喷水装置的的横坐标（最左边为0），ri表示该喷水装置能覆盖的圆的半径。

输出

　　每组测试数据输出一个正整数，表示共需要多少个喷水装置，每个输出单独占一行。

　　如果不存在一种能够把整个草坪湿润的方案，请输出0。

样例输入

2 2 8 6 1 1 4 5 2 10 6 4 5 6 5

样例输出

1 2

C代码如下

#include "stdio.h"
#include "math.h"
typedef struct _Equip//区间的结构体
{
    float L;
    float R;
}Equip;

void Q_sort(Equip **a,int low,int heigh)//对区间按开始端进行升序排序
{
    if(low < heigh)
    {
        int i = low,j = heigh;
        Equip * temp;
        temp = a[i];
        while(i < j)
        {
            while(i < j && a[j]->L > temp->L)
            {
                j --;
            }
            if(i < j)
            {
                a[i] = a[j];
                i++;
            }
            while(i < j && a[i]->L <= temp->L)
            {
                i ++;
            }
            if(i < j)
            {
                a[j] = a[i];
                j --;
            }
        }
        a[i] = temp;
        Q_sort(a,low,i - 1);
        Q_sort(a,i+1,heigh);
    }
}
Equip buf[10000];
Equip *Si[10000];

//寻找从begin开始到end-1之间符合条件的一个区间
//即开始端 <= Rc 的所有小区间中结束端最大的那个小区间 
int findmaxR(int begin,int end)
{
    int k = begain;
    for(int i = begain;i < end;i++)
    {
        if(Si[i]->R >= Si[k]->R)
        {
            k = i;
        }
    }
    return k;
}
int main()
{
    int i,j,m,n,count,s;
    float w,h,Rc,Lc,x,r;
    while(scanf("%d",&m) != EOF)
    {
        for(i = 0;i < m;i++)
        {
            scanf("%d%f%f",&n,&w,&h);
            s = 0;
            for(j = 0;j < n;j++)
            {
                scanf("%f%f",&x,&r);
                if(r > h/2.0)
                {
                    buf[j].L =x - sqrt(r * r - h * h / 4.0);
                    if(buf[j].L < 0.0)
                    {
                        buf[j].L = 0.0;
                    }
                    buf[j].R =x + sqrt(r * r - h * h / 4.0);
                    if(buf[j].R > w)
                    {
                        buf[j].R = w;
                    }
                    Si[s++] = &buf[j];
                }
            }
            Q_sort(Si,0,s-1);//从这里到前边都是准备工作，后边开始贪心了
            Rc = 0.0;
            for(j = 0;j < s && Rc < w;)
            {
                if(Si[j]->L <= Rc)
                {
                    int k;
                    for(k = j;k < s;k++)
                    {
                        if(Si[k]->L > Rc)
                        {
                            break;
                        }
                    }
                    int maxR = findmaxR(j,k);
                    count ++;
                    Rc = Si[maxR]->R;
                    j = maxR +1;
                }
                else
                {
                    break;
                }
            }
            if(Rc < w)
                printf("0\n");
            else
                printf("%d\n",count);
        }
    }
}

作者：风波

mail : fengbohello@qq.com