题目12：喷水装置（二）

题目连接：http://acm.nyist.net/JudgeOnline/problem.php?pid=12
描述
有一块草坪，横向长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

算法思想：这道题目可以抽象成一个区间覆盖问题，可以使用贪心算法实现。

首先将每个装置能覆盖的线段的左边界left和右边界right算出，并对输入的数据进行预处理，将不能将上下边界完全覆盖的删除，然后对数据的按左边界进行排序预处理；
循环遍历，设当前已经覆盖的位置为current，查找left<=current（为了保证能连续覆盖）且右边界right最大的装置（为了保证使用最小装置）。然后将current更新为right，并将计数值加1。
退出循环条件：
1. 遍历到已经没有任何装置可用，退出循环；
2. left>current（因为已经按left排序预处理过，故后面的装置能覆盖到的最左边肯定大于current，不能实现连续覆盖），退出循环。
最后，判断current>left?
如果是，则输出计数值；
否则，输出0（代表不存在方案）。

源代码：

#include <iostream>
#include <vector>
#include <algorithm>
#include <math.h>
using namespace std;
struct  cycle
{
    double a;
    double r;
    double left;
    double right;
};
bool less_cmp(const cycle &x, const cycle &y)
{
    return x.left < y.left;
}

int main()
{
    int N, n, w, h, result;
    vector<cycle> vec;
    cycle cyc;
    double left, right, current;
    cin >> N;
    while (N--)
    {
        cin >> n >> w >> h;
        result = 1;
        vec.clear();
        while (n--)
        {
            cin >> cyc.a >> cyc.r;
            double t = cyc.r * cyc.r - h * h / 4.0;
            if (t <= 0)
            {
                continue;
            }
            left = cyc.a - sqrt(t);
            right = cyc.a + sqrt(t);
            cyc.left = left;
            cyc.right = right;
            vec.push_back(cyc);
        }
        //对圆的最左边进行排序
        sort(vec.begin(),vec.end(),less_cmp);
        current = 0;
        double temp = 0;
        int flag = 0;
        vector<cycle>::iterator it;
        for (it = vec.begin(); it != vec.end() && temp < w; it++)
        {
            if ((*it).left <= current )
            {
                temp = max(temp,(*it).right);
            }
            else if ((*it).left > current && (*it).left <= temp)
            {
                result++;
                current = temp;
                temp = (*it).right;
            }
            else
            {
                break;
            }
        }
        //cout << temp << " @" << w << endl;
        if (temp < w)
            cout <<  0 << endl;
        else
            cout << result << endl;
    }
    return 0;
}

算法时间复杂度：该算法只需循环一遍，故时间复杂度O(n)