#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>
#include <cassert>
#include <vector>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <bitset>
#include <string>
#include <algorithm>
#include <utility>
#define llint long long
#define inf 1e18
#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)
#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))

using namespace std;
typedef pair<llint, llint> P;

llint n, m;
llint a[305][305], c[305];
llint b[305];
vector<P> vec;
llint x[305];

int main(void)
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	
	cin >> n >> m;
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= m; j++){
			cin >> a[i][j];
			c[j] += a[i][j];
		}
	}
	
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= m; j++){
			b[i] += c[j]*a[i][j];
		}
		vec.push_back(P(b[i], i));
	}
	sort(vec.rbegin(), vec.rend());
	
	for(int i = 0; i < vec.size(); i+=2){
		llint id = vec[i].second;
		for(int j = 1; j <= m; j++){
			x[j] += a[id][j];
		}
	}
	
	llint ans = 0;
	for(int i = 1; i <= m; i++){
		ans += 2*c[i]*x[i] - c[i]*c[i];
	}
	cout << ans << endl;
	
	return 0;
}