时间限制: 2000ms 内存限制: 256MB
描述
在一条河的一侧,分布着 N 个村庄。这些村庄平日里需要一些贸易往来,然而商人们来回走遍每一座村庄是非常辛苦的,于是他们决定每个月都在河边举行一次集会,大家都来集会上购买需要的物品。然而在集会地点的选择上,大家却有分歧,因为谁都不愿意集会的地点离自己村庄非常远。经过一番激烈的讨论之后,大家决定要将集会地点挑选在这样一个位置:它离最远的村庄的距离要尽可能的近。
我们把河看做一条足够长的直线,河岸就是平面坐标系上 y = 0 的这条线,y < 0 的区域是河水,而所有村庄都在 y ≥ 0 的区域里。现在给出所有村庄的平面坐标,你要在河岸上找到这样一个位置,使得它到所有村庄的最远距离最小。
输入
输入文件包含多组测试数据。
第一行,给出一个整数 T,为数据组数。接下来依次给出每组测试数据。
每组数据的第一行是一个整数 N,表示村庄的数量。接下来 N 行,每行有两个实数 xi 和 yi,表示每一个村庄的坐标。
输出
对于每组测试数据,输出一行"Case #X: Y",其中 X 表示测试数据编号,Y 表示集会地点的 x 坐标值,要求与正确答案的绝对误差在10-6以内。所有数据按读入顺序从 1 开始编号。
数据范围
小数据:T ≤ 100, 0 < N ≤ 50, 0 ≤ |xi|, yi ≤ 10000
大数据:T ≤ 10, 0 < N ≤ 50000, 0 ≤ |xi|, yi ≤ 10000
-
样例输入
-
1 5 0 8 1 6 4 4 -5 7 -6 1
样例输出
-
Case #1: -1.000000
import java.util.ArrayList;
import java.util.List;
import java.util.Scanner;
public class Main {
public static List<Main.Point> points = new ArrayList<Main.Point>();
public Main(List<Main.Point> list){
points = list;
}
static class Point{
public Point(double x, double y){
this.x = x;
this.y = y;
}
double x;
double y;
}
private double getMaxDis(double curx){
double maxDistance = 0.0;
for(Main.Point p : points){
double dis = Math.pow(curx-p.x, 2) + Math.pow(p.y, 2);
if(dis > maxDistance){
maxDistance = dis;
}
}
return maxDistance;
}
public double getPos(double leftx, double rightx){
if(rightx - leftx <= Math.pow(10.0, -6)){
return (leftx + rightx)/2;
}else{
double maxDisB = getMaxDis((leftx + rightx)/2.0);
double maxDisA = getMaxDis((3*leftx + rightx)/4.0);
double maxDisC = getMaxDis((leftx + 3*rightx)/4.0);
if(maxDisA < maxDisB && maxDisB < maxDisC){
return getPos(leftx, (leftx + 3*rightx)/4.0);
}else if(maxDisA > maxDisB && maxDisB > maxDisC){
return getPos((3*leftx + rightx)/4.0, rightx);
}else{
return getPos((3*leftx + rightx)/4.0, (leftx + 3*rightx)/4.0);
}
}
}
/**
* @param args
*/
public static void main(String[] args) {
Scanner sc = new Scanner(System.in);
int groupNum = sc.nextInt();
double[] result = new double[groupNum];
for(int i = 0; i < groupNum; i++){
int pointNum = sc.nextInt();
List<Main.Point> points = new ArrayList<Main.Point>(pointNum);
double minx = Double.MAX_VALUE;
double maxx = Double.MIN_VALUE;
for(int j = 0; j < pointNum; j++){
double pointx = sc.nextDouble();
double pointy = sc.nextDouble();
if(pointx > maxx){
maxx = pointx;
}else if(pointx < minx){
minx = pointx;
}
points.add(new Main.Point(pointx, pointy));
}
Main main = new Main(points);
double pos = main.getPos(minx, maxx);
result[i] = pos;
}
for(int i = 0; i < result.length; i++){
StringBuilder s = new StringBuilder();
s.append("Case #");
s.append(i+1);
s.append(": ");
s.append(result[i]);
System.out.println(s.toString());
}
}
}
import java.util.ArrayList;
import java.util.List;
import java.util.Scanner;
public class Main {
public static List<Main.Point> points = new ArrayList<Main.Point>();
public Main(List<Main.Point> list){
points = list;
}
static class Point{
public Point(double x, double y){
this.x = x;
this.y = y;
}
double x;
double y;
}
private double getMaxDis(double curx){
double maxDistance = 0.0;
for(Main.Point p : points){
double dis = Math.pow(curx-p.x, 2) + Math.pow(p.y, 2);
if(dis > maxDistance){
maxDistance = dis;
}
}
return maxDistance;
}
public double getPos(double leftx, double rightx){
if(rightx - leftx <= Math.pow(10.0, -6)){
return (leftx + rightx)/2;
}else{
double maxDisB = getMaxDis((leftx + rightx)/2.0);
double maxDisA = getMaxDis((3*leftx + rightx)/4.0);
double maxDisC = getMaxDis((leftx + 3*rightx)/4.0);
if(maxDisA < maxDisB && maxDisB < maxDisC){
return getPos(leftx, (leftx + 3*rightx)/4.0);
}else if(maxDisA > maxDisB && maxDisB > maxDisC){
return getPos((3*leftx + rightx)/4.0, rightx);
}else{
return getPos((3*leftx + rightx)/4.0, (leftx + 3*rightx)/4.0);
}
}
}
/**
* @param args
*/
public static void main(String[] args) {
Scanner sc = new Scanner(System.in);
int groupNum = sc.nextInt();
double[] result = new double[groupNum];
for(int i = 0; i < groupNum; i++){
int pointNum = sc.nextInt();
List<Main.Point> points = new ArrayList<Main.Point>(pointNum);
double minx = Double.MAX_VALUE;
double maxx = Double.MIN_VALUE;
for(int j = 0; j < pointNum; j++){
double pointx = sc.nextDouble();
double pointy = sc.nextDouble();
if(pointx > maxx){
maxx = pointx;
}else if(pointx < minx){
minx = pointx;
}
points.add(new Main.Point(pointx, pointy));
}
Main main = new Main(points);
double pos = main.getPos(minx, maxx);
result[i] = pos;
}
for(int i = 0; i < result.length; i++){
StringBuilder s = new StringBuilder();
s.append("Case #");
s.append(i+1);
s.append(": ");
s.append(result[i]);
System.out.println(s.toString());
}
}
}