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 Rust (1.15.1) (judge verdict: Wrong Answer)

fn main() {
    let (r, w) = (std::io::stdin(), std::io::stdout());
    let mut sc = IO::new(r.lock(), w.lock());
    let n: usize = sc.read();
    let points = (0..n)
        .map(|_| (sc.read::<f64>(), sc.read::<f64>()))
        .collect::<Vec<_>>();

    if n == 2 {
        let (x1, y1) = points[0];
        let (x2, y2) = points[1];
        let dx = x1 - x2;
        let dy = y1 - y2;
        let r = (dx * dx + dy * dy).sqrt() / 2.0;
        println!("{}", r);
        return;
    }

    let mut ans = 1e15;
    for i in 0..n {
        for j in 0..i {
            for k in 0..j {
                let (a, b) = points[i];
                let (c, d) = points[j];
                let (e, f) = points[k];

                let aa = a * a;
                let bb = b * b;
                let cc = c * c;
                let dd = d * d;
                let ee = e * e;
                let ff = f * f;

                let py = ((e - a) * (aa + bb - cc - dd) - (c - a) * (aa + bb - ee - ff))
                    / 2.0
                    / ((e - a) * (b - d) - (c - a) * (b - f));
                let px = if c == a {
                    (2.0 * (b - f) * py - aa - bb + ee + ff) / 2.0 / (e - a)
                } else {
                    (2.0 * (b - d) * py - aa - bb + cc + dd) / 2.0 / (c - a)
                };
                let rr = (px - a) * (px - a) + (py - b) * (py - b);

                let ok = points.iter().all(|&(x, y)| {
                    let dx = x - px;
                    let dy = y - py;
                    dx * dx + dy * dy <= rr
                });
                if ok {
                    let r = rr.sqrt();
                    //                    eprintln!("r={} ijk={:?} o={:?}", r, (i, j, k), (px, py));
                    if ans > r {
                        ans = r;
                    }
                }
            }
        }
    }

    println!("{}", ans);
}

pub struct IO<R, W: std::io::Write>(R, std::io::BufWriter<W>);

impl<R: std::io::Read, W: std::io::Write> IO<R, W> {
    pub fn new(r: R, w: W) -> IO<R, W> {
        IO(r, std::io::BufWriter::new(w))
    }
    pub fn write<S: std::ops::Deref<Target = str>>(&mut self, s: S) {
        use std::io::Write;
        self.1.write(s.as_bytes()).unwrap();
    }
    pub fn read<T: std::str::FromStr>(&mut self) -> T {
        use std::io::Read;
        let buf = self
            .0
            .by_ref()
            .bytes()
            .map(|b| b.unwrap())
            .skip_while(|&b| b == b' ' || b == b'\n' || b == b'\r' || b == b'\t')
            .take_while(|&b| b != b' ' && b != b'\n' && b != b'\r' && b != b'\t')
            .collect::<Vec<_>>();
        unsafe { std::str::from_utf8_unchecked(&buf) }
            .parse()
            .ok()
            .expect("Parse error.")
    }
    pub fn vec<T: std::str::FromStr>(&mut self, n: usize) -> Vec<T> {
        (0..n).map(|_| self.read()).collect()
    }
    pub fn chars(&mut self) -> Vec<char> {
        self.read::<String>().chars().collect()
    }
}