Problem:
Score : 600 points

Problem Statement

Given are N points (x_i, y_i) in a two-dimensional plane.

Find the minimum radius of a circle such that all the points are inside or on it.

Constraints

2 \leq N \leq 50

0 \leq x_i \leq 1000

0 \leq y_i \leq 1000

The given N points are all different.

The values in input are all integers.

Input

Input is given from Standard Input in the following format:

N
x_1 y_1
:
x_N y_N

Output

Print the minimum radius of a circle such that all the N points are inside or on it.

Your output will be considered correct if the absolute or relative error from our answer is at most 10^{-6}.

Sample Input 1

2
0 0
1 0

Sample Output 1

0.500000000000000000

Both points are contained in the circle centered at (0.5,0) with a radius of 0.5.

Sample Input 2

3
0 0
0 1
1 0

Sample Output 2

0.707106781186497524

Sample Input 3

10
10 9
5 9
2 0
0 0
2 7
3 3
2 5
10 0
3 7
1 9

Sample Output 3

6.726812023536805158

If the absolute or relative error from our answer is at most 10^{-6}, the output will be considered correct.

Solution in Kotlin (1.0.0) (judge verdict: Wrong Answer)

import java.lang.Exception

fun main(args: Array<String>) {
    fun calcDist(one: Pair<Double, Double>, two: Pair<Double, Double>) =
            Math.sqrt(Math.pow(one.first - two.first, 2.0) + Math.pow(one.second - two.second, 2.0))

    fun calcCenter(one: Pair<Double, Double>, two: Pair<Double, Double>, three: Pair<Double, Double>): Pair<Double, Double> {
        val a = two.first - one.first
        val b = two.second - one.second
        val c = three.first - one.first
        val d = three.second - one.second

        if ((a * d - b * c) == 0.0) {
            val maxx = Math.max(one.first, Math.max(two.first, three.first))
            val minx = Math.min(one.first, Math.min(two.first, three.first))
            val maxy = Math.max(one.second, Math.max(two.second, three.second))
            val miny = Math.min(one.second, Math.min(two.second, three.second))
            return Pair((maxx + minx) / 2.0, (maxy + miny) / 2.0)
        }

        val ox =
                if (a == 0.0 && c == 0.0) {
                    one.first
                } else {
                    one.first + (d * (a * a + b * b) - b * (c * c + d * d)) / (a * d - b * c) / 2.0
                }
        val oy =
                if (b != 0.0) {
                    (a * (one.first + two.first - ox - ox) + b * (one.second + two.second)) / b / 2.0;
                } else if (d != 0.0) {
                    (c * (one.first + three.first - ox - ox) + d * (one.second + three.second)) / d / 2.0;
                } else {
                    one.second
                }
        return Pair(ox, oy)
    }


    val n = readLine()!!.toInt()
    val xyList = mutableListOf<Pair<Double, Double>>()
    (1..n).forEach {
        val (x, y) = readLine()!!.split(' ').map(String::toInt)
        xyList.add(Pair(x.toDouble(), y.toDouble()))
    }
    if (n == 2) {
        println(calcDist(xyList[0], xyList[1]) / 2.0)
        return
    }

    var minDist = Double.MAX_VALUE
    for (i in 0 until n - 2) {
        for (j in i + 1 until n - 1) {
            for (k in j + 1 until n) {
                try {
                    val center = calcCenter(xyList[i], xyList[j], xyList[k])
                    if (center.first.isNaN() || center.second.isNaN()
                            || center.first.isInfinite() || center.second.isInfinite()) continue
                    var harfDist = calcDist(xyList[i], center)
                    harfDist = Math.max(harfDist, calcDist(xyList[j], center))
                    harfDist = Math.max(harfDist, calcDist(xyList[k], center))
                    if (xyList.all { calcDist(it, center) <= harfDist }) {
                        minDist = Math.min(minDist, harfDist)
                    }
                } catch (e: Exception) {
//                    do nothing
                }
            }
        }
    }
    for (i in 0 until n - 1) {
        for (j in 1 until n) {
            val calcDist = calcDist(xyList[i], xyList[j])
            val center = Pair((xyList[i].second + xyList[j].second) / 2.0, (xyList[i].second + xyList[j].second) / 2.0)
            if (xyList.all { calcDist(it, center) <= calcDist }) {
                minDist = Math.min(minDist, calcDist)
            }
        }
    }
    println(minDist)
}