2019网易笔记--最小覆盖圆

2019网易春招第三题，用到了类似最小覆盖圆的原理
题目大致和这个题相同：
【题目描述】
H族的领地是一片富饶之地，有丰富的资源，无数的黄金。尽管异族侵略者被压制几百多年，他们并没有放弃对这片土地的渴望。最近，他们研制了一枚重磅炸弹。试图用它摧毁H族的防御工事。H族的防御工事由多个基地构成，侵略者试图用炸弹摧毁尽量多的防御基地，现在，他们急需你的帮助。如果你能为他们解决难题，会有巨大的奖励。（该炸弹的攻击范围是一个半径为r单位的圆）

【输入】

第一行为样例的个数T。
每个样例中有若干行，第一行为炸弹的攻击半径r，r有可能不是整数(1<=r<=100)
第二行是一个整数n，表示H族防御基地的个数(1<=n<=100)
接下来n行，每行有两个数，代表H族某个防御基地的坐标
【输出】

只有一个数，代表炸弹能最多摧毁的H族防御基地数量。
【输入样例】

1
1
3
1 0
0 1
4 0
【输出样例】

2
思路：
暴力法，O（nnn）
前面n*n确定假想圆的圆心，然后再逐个点判断这个点是不是在假想圆里，同时记录最大数量

#define eps 1e-8
#include <iostream>
#include <math.h>

using namespace std;

struct point
{
	int x;
	int y;
	point(){}
	point(double zx, double zy) { x = zx; y = zy; }
}p[100];

double distance(point a, point b)
{
	return sqrt(pow((a.x - b.x), 2) + pow((a.y - b.y), 2));
}

point getCircleCenter(point a, point b,int r)
{
	point mid=point((a.x + b.x) / 2, (a.y + b.y) / 2);
	double angle = atan2(a.x - b.x, a.y - b.y);
	double d = sqrt(r*r - distance(mid, a)*distance(mid, a));
	return point(mid.x + d*cos(angle), mid.y + d*sin(angle));
}

int main()
{
	int r, num,res=INT_MIN;
	cin >> r >> num;
	for (int i = 0; i < num; ++i)
		cin >> p[i].x >> p[i].y;
	for (int i = 0; i < num; ++i)
	{
		for (int j = i + 1; j < num; ++j)
		{
			if (distance(p[i], p[j])>2.0*r+eps)continue;
			point center=getCircleCenter(p[i], p[j],r);
			int count = 0;
			for (int k = 0; k < num; ++k)
				if (distance(center, p[k]) < 1.0*r + eps)
					count++;
			res = (res > count) ? res : count;
		}
	}
	cout << res << endl;
	system("pause");
}