特定半径包围点数

链接：https://ac.nowcoder.com/acm/contest/18428/D
来源：牛客网

天快黑了，蚊子正在攻击渔夫弗洛伊德。弗洛伊德决定最后一次投掷他的圆网，并希望充分利用最后一次投掷。

给定 Floyd 网的大小（半径）和一组鱼的位置 (x,y)，Floyd 一次投掷最多可以钓到多少鱼？也就是说，找到您可以在网络（圆圈）中拥有的最大点数。如果点（鱼）在圆边界的 10-6 范围内，则考虑圆中的点。

输入描述：
第一个输入行提供圆的半径。第二个输入行包含一个整数 n (1 ≤ n ≤ 100)，表示鱼的数量。接下来的 n 个输入行中的每一行都提供一条鱼的位置 (x,y)。

假设所有这些点都是不同的。还假设输入中的所有 x 和 y 值（以及圆半径）都是 1 到 1000 之间的整数，包括 1 和 1000。
输出描述：
打印 Floyd 一次可以钓到的最大鱼。
示例1
输入
20
4
6 7
5 5
190 100
4 4
输出
3
示例2
输入
3
2
1 1
95 4
输出
1

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const double esp=1e-6;
ll r;
struct point1
{
	double x,y;
	point1(){}
	point1(double tx,double ty){x=tx,y=ty;}
}point[1009];
double dis(point1 a,point1 b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
point1 center(point1 a,point1 b)
{
	point1 mid;
	mid.x=(a.x+b.x)/2;
	mid.y=(a.y+b.y)/2;
	double angle=atan2(a.x-b.x,a.y-b.y);
	double d=sqrt(r*r-dis(a,mid)*dis(a,mid));
	return point1(mid.x+d*cos(angle),mid.y+d*sin(angle));
}
int main()
{
	cin>>r;
	ll n;
	cin>>n;
	ll ans=1,cnt;
	for(int i=1;i<=n;i++) cin>>point[i].x>>point[i].y;
	for(int i=1;i<=n;i++)
	{
		for(int j=i+1;j<=n;j++)
		{
			if(dis(point[i],point[j])>2.0*r) continue;
			point1 c=center(point[i],point[j]);
			 cnt=0;
			for(int k=1;k<=n;k++)
			{
				if(dis(c,point[k])<1.0*r+esp) cnt++;
			}
			ans=max(ans,cnt);
		}
	}
	cout<<ans<<endl;
	return 0;
}