#include <bits/stdc++.h>
using namespace std;
using lint = long long int;
struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_;
#define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i<i##_end_;i++)
#define REP(i, n) FOR(i,0,n)
template<typename T> istream &operator>>(istream &is, vector<T> &vec){ for (auto &v : vec) is >> v; return is; }
#define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl;

int main()
{
    int N, M;
    cin >> N >> M;
    vector<vector<lint>> A(N, vector<lint>(M));
    cin >> A;
    vector<lint> tot(M);
    REP(i, N) REP(j, M) tot[j] += A[i][j];

    vector<pair<lint, vector<lint>>> vp;
    REP(i, N)
    {
        lint val = 0;
        REP(j, M) val += tot[j] * A[i][j];
        vp.emplace_back(-val, A[i]);
    }
    sort(vp.begin(), vp.end());
    vector<vector<lint>> d(2, vector<lint>(M));
    bool b = false;
    for (auto x : vp)
    {
        REP(i, M) d[b][i] += x.second[i];
        b = !b;
    }
    lint ret = 0;
    REP(i, M) ret += d[0][i] * d[0][i] - d[1][i] * d[1][i];
    cout << ret << '\n';
}