POJ1981 Circle and Points（计算几何）

题意：

给出一系列点，这些点之间的距离都在2附近，求一个单位圆最多能包含几个点。

要点：

极限情况就是2个点在圆的边上，这样就可以求出圆心，再枚举看其他点在不在圆内即可，复杂度是n^3。

#include<cstdio>  
#include<cstdlib>  
#include<cstring>  
#include<cmath>  
#include<algorithm>  
#define eps 1e-8
using namespace std;
struct node
{
	double x, y;
}p[310];

double dis(node a, node b)
{
	return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y));
}
void get_center(node a, node b, node& center)//已知两点和半径求圆心
{
	node mid;
	mid.x = (a.x + b.x) / 2.0;
	mid.y = (a.y + b.y) / 2.0;
	double angle = atan2(a.x - b.x, b.y - a.y);//注意这里
	double d = sqrt(1 - dis(mid, a)*dis(mid, a));
	center.x = mid.x + d*cos(angle);
	center.y = mid.y + d*sin(angle);
	//printf("%lf   %lf\n", center.x, center.y);
}

int main()
{
	int n, i, j, k;
	while (scanf("%d", &n), n)
	{
		for (i = 1; i <= n; i++)
			scanf("%lf%lf", &p[i].x,&p[i].y);
		int ans = 1;
		for(i=1;i<=n;i++)
			for (j = i + 1; j <= n; j++)
			{
				if (dis(p[i], p[j]) > 2.0) continue;
				node center;
				get_center(p[i], p[j],center);
				int cnt = 0;
				for (k = 1; k <= n; k++)
				{
					if (dis(center, p[k]) < 1.0 + eps)
						cnt++;
				}
				ans = max(cnt, ans);
			}
		printf("%d\n", ans);
	}
	return 0;
}