2013 编程之美挑战赛 集会

时间限制: 2000ms 内存限制: 256MB

描述

在一条河的一侧，分布着 N 个村庄。这些村庄平日里需要一些贸易往来，然而商人们来回走遍每一座村庄是非常辛苦的，于是他们决定每个月都在河边举行一次集会，大家都来集会上购买需要的物品。然而在集会地点的选择上，大家却有分歧，因为谁都不愿意集会的地点离自己村庄非常远。经过一番激烈的讨论之后，大家决定要将集会地点挑选在这样一个位置：它离最远的村庄的距离要尽可能的近。

我们把河看做一条足够长的直线，河岸就是平面坐标系上 y = 0 的这条线，y < 0 的区域是河水，而所有村庄都在 y ≥ 0 的区域里。现在给出所有村庄的平面坐标，你要在河岸上找到这样一个位置，使得它到所有村庄的最远距离最小。

输入

输入文件包含多组测试数据。

第一行，给出一个整数 T，为数据组数。接下来依次给出每组测试数据。

每组数据的第一行是一个整数 N，表示村庄的数量。接下来 N 行，每行有两个实数 x_i 和 y_i，表示每一个村庄的坐标。

输出

对于每组测试数据，输出一行"Case #X: Y"，其中 X 表示测试数据编号，Y 表示集会地点的 x 坐标值，要求与正确答案的绝对误差在10^-6以内。所有数据按读入顺序从 1 开始编号。

数据范围

小数据：T ≤ 100, 0 < N ≤ 50, 0 ≤ |x_i|, y_i ≤ 10000

大数据：T ≤ 10, 0 < N ≤ 50000, 0 ≤ |x_i|, y_i ≤ 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());
		}
	}

}