Description
Assume the coasting is an infinite straight line. Land is in one side of coasting, sea in the other. Each small island is a point locating in the sea side. And any radar installation, locating on the coasting, can only cover d distance, so an island in the sea can be covered by a radius installation, if the distance between them is at most d.
We use Cartesian coordinate system, defining the coasting is the x-axis. The sea side is above x-axis, and the land side below. Given the position of each island in the sea, and given the distance of the coverage of the radar installation, your task is to write a program to find the minimal number of radar installations to cover all the islands. Note that the position of an island is represented by its x-y coordinates.
Figure A Sample Input of Radar Installations
We use Cartesian coordinate system, defining the coasting is the x-axis. The sea side is above x-axis, and the land side below. Given the position of each island in the sea, and given the distance of the coverage of the radar installation, your task is to write a program to find the minimal number of radar installations to cover all the islands. Note that the position of an island is represented by its x-y coordinates.
Figure A Sample Input of Radar Installations
Input
The input consists of several test cases. The first line of each case contains two integers n (1<=n<=1000) and d, where n is the number of islands in the sea and d is the distance of coverage of the radar installation. This is followed by n lines each containing two integers representing the coordinate of the position of each island. Then a blank line follows to separate the cases.
The input is terminated by a line containing pair of zeros
The input is terminated by a line containing pair of zeros
Output
For each test case output one line consisting of the test case number followed by the minimal number of radar installations needed. "-1" installation means no solution for that case.
Sample Input
3 2 1 2 -3 1 2 1 1 2 0 2 0 0
Sample Output
Case 1: 2 Case 2: 1
#include <cmath> #include <algorithm> #include <cstdio> using namespace std; struct node { double le, r; } is[1200]; bool cmp(node a, node b) { return a.le < b.le; } int main() { int n, i, num, t; double x, y, temp, d; int flag = 1; while (~scanf("%d %lf", &n, &d) && !(n==0&&d==0)) { t = 0; for (i = 0; i < n; i++) { scanf("%lf %lf", &x, &y); if (y > d || d < 0) { t = 1; } is[i].le = x - sqrt(d * d - y * y); //重点&&难点:将面积转换为长度 is[i].r = x + sqrt(d * d - y * y); } if (t == 1) { printf("Case %d: -1\n", flag); } else //主要算法 { num = 1; sort(is, is + n, cmp); temp = is[0].r; for(i=1;i<n;i++) { if(is[i].r<=temp) { temp=is[i].r; } else if(is[i].le > temp) { num++; temp=is[i].r; } } printf("Case %d: %d\n", flag, num); } // end 算法 flag++; } return 0; }