POJ1328《Radar Installation》方法:贪心+快排

题目描述：在x轴上布置雷达，当雷达数最少的时候，使得可以包含全部的岛屿。

方法：首先将岛屿以雷达为半径，转化为x轴上的一个个区间，然后对这些区间进行排序。雷达第一次放在最左边雷达的右端点，然后对下一个雷达的左端点和之前的右端点比较，若该左端点在参照雷达右端点的左端，需要判断对应右端点在雷达右端点的左边还是右边，若在左边，则参照值被赋予该右端点，否则，继续判断下一个端点。若在右边，直接加雷达。

#include <iostream> #include <cstdio> #include <cmath> using namespace std; int part(double p[1000][2], int low, int high) //霍尔分割 { double sentinel[2]; sentinel[0] = p[low][0]; sentinel[1] = p[low][1]; while (low < high) { while (p[high][0] >= sentinel[0] && low < high) --high; p[low][0] = p[high][0]; p[low][1] = p[high][1]; while (p[low][0] <= sentinel[0] && low < high) ++low; p[high][0] = p[low][0]; p[high][1] = p[low][1]; } p[low][0] = sentinel[0]; p[low][1] = sentinel[1]; return low; } void quicksort(double p[1000][2], int left, int right) { if (left < right) { int i = part(p, left, right); quicksort(p, left, i-1); quicksort(p, i+1, right); } } int addlist(double p[1000][2], int length, int r) { int x, y, cnt; double temp; bool isOut = false; for (int i = 0; i < length; ++i) { cin >> x >> y; if (y > r || r < 0) isOut = true; if (isOut == false) { p[i][0] = x - sqrt(r*r-y*y); p[i][1] = x + sqrt(r*r-y*y); } } if (isOut) return -1; quicksort(p, 0, length-1); cnt = 1; temp = p[0][1]; //第一个雷达在默认在最左边界的右端点 for (int i = 1; i < length; ++i) { if (p[i][0] <= temp) { if (p[i][1] < temp) temp = p[i][1]; } else { ++cnt; temp = p[i][1]; } } return cnt; } int main() { int t, r; //t总数, r半径 int cnt = 0; //freopen("temp.txt", "r", stdin); while (true) { double p[1000][2]; //处理后的左端点, 右端点 cin >> t >> r; if (t == 0 && r == 0) break; cout << "Case " << ++cnt << ": " << addlist(p, t, r) << endl; } return 0; }