算法实验 通信覆盖

一、 实验题目
• 题目描述：
假设海岸线是一条无限延伸的直线，陆地在海岸线的一侧，海洋在另一侧。每一个小岛相当于海洋侧的一个点，坐落在海岸线的基站只能覆盖半径为d的范围。应用直角坐标系，将海岸线作为x轴，设海洋在x轴上方。给定海洋中各小岛的位置，以及通信基站的覆盖半径，要求用最少的基站覆盖所有小岛，使得每个小岛都能和陆地通过某基站通信。

• 输入:
小岛的数目，基站的覆盖半径，各小岛的坐标（第一行为小岛数与通信半径，后面各行为小岛坐标）
• 输出：
覆盖小岛的最少基站数目及基站坐标，如果无解，输出-1
• 实验要求：

1. 用贪心算法解决该问题；
2. 至少运行3组不同的数据。
   二、 实验实现
   1． 编程语言:
   实验用C完成算法设计和程序设计并调试通过。
   2．解题方法：贪心算法
   将小岛的x坐标按从小到大排序x1<x2<…<xn，以左边第一个小岛作为第一个基站的边界作半径为d的圆，圆心为c1=x1+sqrt（dd-y1y1），从左到右检查每个小岛，如果小岛在该圆的覆盖范围内，则检查下一个；如果小岛不在覆盖范围内，则以该点为边界点做下一个圆c2。直至所有点都被覆盖。
3. 证明：
   先假设小岛的坐标集合S={x1,x2,…,xn}已经按x轴从小到大排序。
   第一次选取横坐标最小的点作为第一个圆周上的点，圆心为c1（c1=x1+sqrt（dd-y1y1），0），S’={xi，…，xn}为去除第一个圆所包含的点的集合。此时选取x1作为左侧圆周上的点，重复上述操作，直到所有的点都在圆中。因为每次选取子问题 的点可以得到最优解O={a1，a2，…，ai}。
   a=1c1，当a2！=c2，将去除集合中的a2换成c2，这样既能覆盖点x2，又能覆盖右边最多的点，因此选取c2覆盖的点不会更少，则解不会更差，因此贪心选择它。
4. 结果截图
   输入1： 输入2： 输入3：

结果1： 结果2： 结果3：

5. 时间复杂度：
   算法运行过程中，只需检索每个小岛，小岛检索完程序也就运行结束，所以时间复杂度为O（n）。
   四、源程序

#include<iostream>
#include<algorithm>
#include<cmath>

using namespace std;

int n, num;
double d, len;

struct Node {
	double x;
	double y;
};

struct Node a[100];

double Center[100];
int F[100];


bool squence(Node a, Node b) {
	if (a.x < b.x) return 1;
	return 0;
}

void Greedy() {
	//点无法覆盖
	for (int i = 0; i < n; i++) {
		if (a[i].y > d) printf_s("-1");
	}
	len = a[0].x + sqrt(double(d * d) - a[0].y * a[0].y);
	Center[0] = len;
	F[0] = 1;
	for (int i = 1; i < n; i++) {
		if ((a[i].x - len) * (a[i].x - len) + a[i].y * a[i].y > d* d) {
			len = a[i].x + sqrt(double(d * d) - a[i].y * a[i].y);
			Center[i] = len;
			F[i] = 1;
			num++;
		}
		else if ((a[i].x - len) * (a[i].x - len) + a[i].y * a[i].y <= d * d) {
			Center[i] = len;
			continue;
		}
	}
}

int main() {
	memset(F, 0, sizeof(F));
	cin >> n >> d;
	for (int i = 0; i < n; i++) {
		cin >> a[i].x >> a[i].y;
	}

	sort(a, a + n, squence);

	num = 1;
	Greedy();
	printf_s("%d\n", num);
	for (int i = 0; i < n; i++) {
		if (F[i] == 1) {
			printf_s("%f 0\n", Center[i]);
		}
	}
	return 0;
}