const INF: f64 = (1usize << 60) as f64; const ITER: usize = 200; fn calc(d: f64, a: &Vec) -> (f64, f64) { let mut lower = -INF; let mut upper = INF; for _ in 0..ITER { let lmiddle = (upper + 2. * lower) / 3.; let umiddle = (2. * upper + lower) / 3.; let lucost = calc_c(d, lmiddle, a); let ulcost = calc_c(d, umiddle, a); if lucost > ulcost { lower = lmiddle; } else { upper = umiddle; } } (calc_c(d, upper, a), upper) } fn calc_c(d: f64, x: f64, a: &Vec) -> f64 { (0..a.len()).map(|i| { (d * i as f64 + x - a[i]) * (d * i as f64 + x - a[i]) }) .sum::() } fn main() { let mut n = String::new(); std::io::stdin().read_line(&mut n).ok(); let n: usize = n.trim().parse().unwrap(); let mut a = String::new(); std::io::stdin().read_line(&mut a).ok(); let a: Vec = a.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let mut lower_d = -INF; let mut upper_d = INF; let mut result_d = 0.; let mut result_x = 0.; for _ in 0..ITER { let middle_ld = (upper_d + 2. * lower_d) / 3.; let middle_ud = (2. * upper_d + lower_d) / 3.; let (lucost, lux) = calc(middle_ld, &a); let (ulcost, ulx) = calc(middle_ud, &a); if lucost > ulcost { lower_d = middle_ld; result_d = middle_ud; result_x = ulx; } else { upper_d = middle_ud; result_d = middle_ld; result_x = lux; } } println!("{} {}", result_x, result_d); println!("{}", calc_c(result_d, result_x, &a)); }