/** * author: otera **/ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define int long long typedef long long ll; typedef unsigned long long ul; typedef unsigned int ui; typedef long double ld; const int inf=1e9+7; const ll INF=1LL<<60 ; const ll mod=1e9+7 ; #define rep(i,n) for(int i=0;i=0;i--) #define Rep(i,sta,n) for(int i=sta;i=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef complex Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair P; typedef pair LDP; typedef pair LP; #define fr first #define sc second #define all(c) c.begin(),c.end() #define pb push_back #define debug(x) cerr << #x << " = " << (x) << endl; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } void solve() { int n, m; cin >> n >> m; vector> a(n, vector(m + 1)); rep(i, n) { int sum = 0; rep(j, m) { cin >> a[i][j]; sum += a[i][j]; } a[i][m] = sum; } sort(all(a), [&](vector x, vector y){ return x[m] > y[m];}); vector xcnt(m, 0), ycnt(m, 0); rep(i, n) { if(i % 2 == 0) { rep(j, m) { xcnt[j] += a[i][j]; } } else { rep(j, m) { ycnt[j] += a[i][j]; } } } int ans = 0; rep(i, m) { ans += xcnt[i] * xcnt[i]; } rep(i, m) { ans -= ycnt[i] * ycnt[i]; } cout << ans << endl; } signed main() { ios::sync_with_stdio(false); cin.tie(0); //cout << fixed << setprecision(10); //int t; cin >> t; rep(i, t)solve(); solve(); return 0; }