计算几何 之 hdu 1077 poj 1981 O(n^3)

注意两题输入格式不同。。。O(n^3)，确实有更优的解法 O(n^2logn)   详解+代码，请点击这里

/*
 
题目： 用半径为 1.0001 的圆，尽可能多的包含所给点。
 
第一次做这样的计算几何题，，，感觉挺麻烦的，，特别精度的把握上。。。
 
思路：
1）任意两点间的距离大于直径（2.0002），则答案为 1 ；
2）若不满足1),则在圆上必存在两个点，而这两个点也正好确定了这个圆，
	故遍历所有点，看有多少个点在圆中，满足条件点的个数即答案。
 
关键：
判断点是否在圆中： 可用此点到原点的距离与半径作比较判断
故关键在于如何求原点？
	两个点确定中垂线 + 半径长 => 原点 （具体实现看代码，一些细节东西需要注意）
 
另外：
由于两个点可以确定两个圆心，但取一个就好。
（
严格地数学证明，没找到；但是，想一下三个点的情况：
假设三个点分别为a、b、c,若三个点在同一个半径为1.0001的圆内，
即使a、b没包含c，a、c没包含b，但是b、c一定可以包含a,
前提：代码中求任意两个点所对应的圆心时，或者都 '+' ，或者都 '-' （具体看函数 void Get_Center(int i, int j)）
）
 
*/

// hdu 1077 
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
const int MAX = 305;
const double eps = 1e-6;

double arr_p[MAX][2];
double center[2];

double Distance(double p1[], double p2[])
{
	return sqrt((p2[0] - p1[0])*(p2[0] - p1[0]) + (p2[1] - p1[1])*(p2[1] - p1[1]));
}

void Get_Center(int i, int j){
	double p_mid[2];
	p_mid[0] = (arr_p[i][0] + arr_p[j][0]) / 2.0;
	p_mid[1] = (arr_p[i][1] + arr_p[j][1]) / 2.0;

	double dis_x, dis_y;
	dis_x = arr_p[i][0] - arr_p[j][0];
	dis_y = arr_p[i][1] - arr_p[j][1];

	double dis_mid_center = sqrt(1 - ((p_mid[0] - arr_p[i][0])*(p_mid[0] - arr_p[i][0]) + (p_mid[1] - arr_p[i][1])*(p_mid[1] - arr_p[i][1])));

	if (fabs(dis_y) < eps) { // 中垂线垂直于 x 轴
		center[0] = p_mid[0];
		center[1] = p_mid[1] + dis_mid_center;
		/*
		center[0] = p_mid[0];
		center[1] = p_mid[1] - dis_mid_center;
		*/
	}
	else {
		double k = atan(-(dis_x / dis_y));
		double tep1 = dis_mid_center*cos(k);
		double tep2 = dis_mid_center*sin(k);
		center[0] = p_mid[0] + tep1;  // 皆采取 '+' 的方式
		center[1] = p_mid[1] + tep2;
		/*
		center[0] = p_mid[0] - tep1;  // 皆采取 '-' 的方式
		center[1] = p_mid[1] - tep2;
		*/
	}
}

int main()
{
	//freopen("input.txt", "r", stdin);
	int T;
	scanf("%d", &T);
	while (T--) {
		int N, tep;
		scanf("%d", &N);
		for (int i = 0; i < N; i++) {
			scanf("%lf %lf", &arr_p[i][0], &arr_p[i][1]);
		}
		int ans = 1;
		for (int i = 0; i < N; i++) {
			for (int j = i + 1; j < N; j++) {
				if (Distance(arr_p[i], arr_p[j]) > 2.0) {
					continue;
				}

				Get_Center(i, j);

				tep = 2;
				for (int k = 0; k < N; k++) {
					if (k == i || k == j) {
						continue;
					}
					if (Distance(arr_p[k], center) <= 1.0001) {
						tep++;
					}
				}
				ans = max(ans, tep);
			}
		}
		printf("%d\n", ans);
	}
	return 0;
}

// poj 1981
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
const int MAX = 305;
const double eps = 1e-6;

double arr_p[MAX][2];
double center[2];

double Distance(double p1[], double p2[])
{
	return sqrt((p2[0] - p1[0])*(p2[0] - p1[0]) + (p2[1] - p1[1])*(p2[1] - p1[1]));
}

void Get_Center(int i, int j){
	double p_mid[2];
	p_mid[0] = (arr_p[i][0] + arr_p[j][0]) / 2.0;
	p_mid[1] = (arr_p[i][1] + arr_p[j][1]) / 2.0;

	double dis_x, dis_y;
	dis_x = arr_p[i][0] - arr_p[j][0];
	dis_y = arr_p[i][1] - arr_p[j][1];

	double dis_mid_center = sqrt(1 - ((p_mid[0] - arr_p[i][0])*(p_mid[0] - arr_p[i][0]) + (p_mid[1] - arr_p[i][1])*(p_mid[1] - arr_p[i][1])));

	if (fabs(dis_y) < eps) {
		center[0] = p_mid[0];
		center[1] = p_mid[1] + dis_mid_center;
	}
	else {
		double k = atan(-(dis_x / dis_y));
		double tep1 = dis_mid_center*cos(k);
		double tep2 = dis_mid_center*sin(k);
		center[0] = p_mid[0] + tep1;
		center[1] = p_mid[1] + tep2;
	}
}

int main()
{
	int N;
	while (scanf("%d", &N) && N) {
		int tep;

		for (int i = 0; i < N; i++) {
			scanf("%lf %lf", &arr_p[i][0], &arr_p[i][1]);
		}
		int ans = 1;
		for (int i = 0; i < N; i++) {
			for (int j = i + 1; j < N; j++) {
				if (Distance(arr_p[i], arr_p[j]) > 2.0002) {
					continue;
				}

				Get_Center(i, j);

				tep = 2;
				for (int k = 0; k < N; k++) {
					if (k == i || k == j) {
						continue;
					}
					if (Distance(arr_p[k], center) <= 1.0001) {
						tep++;
					}
				}
				ans = max(ans, tep);
			}
		}
		printf("%d\n", ans);
	}
	return 0;
}