#include <bits/stdc++.h>
using namespace std;
int n, m;
int t[1000][1000];
long long sol;
long long sqr(long long x) { return x * x; }
long long a[1000];
int opt1D(int N) {
  int ret(0);
  long long retv(1000000000000000000ll);
  for (int i = 0; i <= N; i++) {
    long long cv(0ll);
    int p0(i * 4);
    for (int j = 0; j < N; j++) {
      int p1((j * 4) + 2);
      cv += sqr(p1 - p0) * a[j];
    }
    if (cv < retv) {
      retv = cv;
      ret = i;
    }
  }
  sol += retv;
  return ret;
}
int main() {
  scanf("%d%d", &n, &m);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < m; j++) {
      scanf("%d", &t[i][j]);
      a[i] += t[i][j];
    }
  }
  int x = opt1D(n);
  for (int i = 0; i < m; i++) {
    a[i] = 0;
    for (int j = 0; j < n; j++) a[i] += t[j][i];
  }
  int y = opt1D(m);
  printf("%I64d\n%d %d\n", sol, x, y);
  return 0;
}
