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):
    if turn == 0:
        tmp = []
        for i in range(N):
            if used[i]:
                continue
            s = 0
            for j, a in enumerate(A[i]):
                s += a * a + 2 * X[j] * a
            tmp.append((s, i))
        tmp.sort()
        i = tmp[-1][1]
        used[i] = 1
        for j, a in enumerate(A[i]):
            X[j] += a
    else:
        tmp = []
        for i in range(N):
            if used[i]:
                continue
            s = 0
            for j, a in enumerate(A[i]):
                s += a * a + 2 * Y[j] * a
            tmp.append((s, i))
        tmp.sort()
        i = tmp[-1][1]
        used[i] = 1
        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)