注意两题输入格式不同。。。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;
}