//https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8 より macro_rules! input { (source = $s:expr, $($r:tt)*) => { let mut iter = $s.split_whitespace(); input_inner!{iter, $($r)*} }; ($($r:tt)*) => { let s = { use std::io::Read; let mut s = String::new(); std::io::stdin().read_to_string(&mut s).unwrap(); s }; let mut iter = s.split_whitespace(); input_inner!{iter, $($r)*} }; } macro_rules! input_inner { ($iter:expr) => {}; ($iter:expr, ) => {}; ($iter:expr, $var:ident : $t:tt $($r:tt)*) => { let $var = read_value!($iter, $t); input_inner!{$iter $($r)*} }; } macro_rules! read_value { ($iter:expr, ( $($t:tt),* )) => { ( $(read_value!($iter, $t)),* ) }; ($iter:expr, [ $t:tt ; $len:expr ]) => { (0..$len).map(|_| read_value!($iter, $t)).collect::>() }; ($iter:expr, chars) => { read_value!($iter, String).chars().collect::>() }; ($iter:expr, usize1) => { read_value!($iter, usize) - 1 }; ($iter:expr, $t:ty) => { $iter.next().unwrap().parse::<$t>().expect("Parse error") }; } // fn run() { input! { n: usize, m: usize, a: [[i64; m]; n], } let mut x = vec![0; m]; for a in a.iter() { for (x, a) in x.iter_mut().zip(a.iter()) { *x += *a; } } let mut y = vec![0; m]; let mut a = a; let mut ans = x.iter().fold(0, |s, a| s + a.pow(2)); for _ in 0..(n / 2) { let mut b = vec![]; for (i, a) in a.iter().enumerate() { let mut sum = 0; for (a, (x, y)) in a.iter().zip(x.iter().zip(y.iter())) { sum += -2 * *a * (*x + *y); } b.push((sum, i)); } b.sort(); let p = b[0].1; let q = b[1].1; ans += b[1].0; for (x, (y, a)) in x.iter_mut().zip(y.iter_mut().zip(a[q].iter())) { *x -= *a; *y += *a; } if p > q { a.remove(p); a.remove(q); } else { a.remove(q); a.remove(p); } } println!("{}", ans); } fn main() { run(); }