hdu 1077 圆覆盖的最大点数

题目大意：给你一个半径为1的圆，再给你n个点，

求这个圆最多覆盖n个点中的几个。

 

思路：枚举两个在圆上的点确定一个圆的圆心，

然后枚举出一共有多少点能被这个圆覆盖，求出其中最多的数目。

注意：要注意会用在圆上的两个点以及圆的半径求出圆心！

#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define INF 99999999
#define uf 1e-8 
#define maxn 400
#define LL long long
using namespace std;
typedef struct node
{
	double x,y;
}pointy;
double dis(pointy a,pointy b)//求两点距离的平方 
{
	double ans;
	ans=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
	return ans;
}
pointy circlepoint(pointy a,pointy b)//求圆心 
{
	pointy m,z;
	double xx=a.x-b.x;
	double yy=a.y-b.y;
	m.x=(a.x+b.x)/2.0;
	m.y=(a.y+b.y)/2.0;
	double an=sqrt(1.0-dis(a,m));
	if(fabs(xx)<uf)//如果两点两线与x轴接近平行，此时tan a 为无穷大 
	{
		z.x=m.x;
		z.y=m.y+an;
	}
	else
	{
		double ang=atan(-xx/yy);
		z.x=m.x+an*cos(ang);
		z.y=m.y+an*sin(ang);
	}
	return z;
}
int main()
{
	int T;
	cin>>T;
	while(T--)
	{
		int n;
		cin>>n;
		pointy a[maxn];
		for(int i=0;i<n;i++)
		{
			cin>>a[i].x>>a[i].y;
		}
		int ans=1;
		int cnt;
		for(int i=0;i<n;i++)
		{
			for(int j=i+1;j<n;j++)
			{
				if(dis(a[i],a[j])<=4.0)//如果两点可以在圆上 
				{
					pointy z=circlepoint(a[i],a[j]);
					cnt=0;
					for(int k=0;k<n;k++)
					{
						if(sqrt(dis(z,a[k]))<=1.00000001)
							cnt++;
					}
					if(ans<cnt)
						ans=cnt;
				}
			}
		}
		cout << ans << endl;
	}
	return 0;	
}