import sys
input = sys.stdin.buffer.readline
sys.setrecursionlimit(10 ** 7)

N, M = map(int, input().split())
X = [0]*M
Y = [0]*M

used = [0]*N
A = tuple(tuple(map(int, input().split())) for _ in range(N))

turn = 0
for _ in range(N):
    tmp = []
    for i in range(N):
        if used[i]:
            continue
        s = 0
        t = 0
        for j, a in enumerate(A[i]):
            s += a * (X[j] + Y[j])
            t += a * a + X[j] * a * 2
        tmp.append((s, t, i))
    tmp.sort()
    i = tmp[-1][2]
    used[i] = 1
    if turn == 0:
        for j, a in enumerate(A[i]):
            X[j] += a
    else:
        for j, a in enumerate(A[i]):
            Y[j] += a
    turn = 1 - turn

ans = sum(x * x for x in X) - sum(y * y for y in Y)
print(ans)