题目连接: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)