#include <bits/stdc++.h>
using namespace std;
const long long INF = 1e12;
const int MAXN = 1e5 + 1, MODM = 1e9 + 7, MULT1 = 30, MULT2 = 239, LUCKY = 13,
          ELEVEN = 11, MAXM = 101;
struct cl {
  long long val, x, y;
  cl() { val = -1; }
};
long long mas[2001][2001];
long long le[2001][2001], rr[2][5001];
long long ans[2001][2001];
cl f_mx[2][2], s_mx[2][2];
int main() {
  ios_base::sync_with_stdio(0);
  int n;
  cin >> n;
  for (int i = 0; i < n; i++)
    for (int j = 0; j < n; j++) {
      cin >> mas[i][j];
      le[i - min(i, j)][j - min(i, j)] += mas[i][j];
      rr[0][j + i] += mas[i][j];
    }
  for (int i = 0; i < n; i++)
    for (int j = 0; j < n; j++) {
      ans[i][j] = -mas[i][j] + le[i - min(i, j)][j - min(i, j)] + rr[0][j + i];
      int i1 = i % 2, j1 = j % 2;
      if (f_mx[i1][j1].val < ans[i][j]) {
        f_mx[i1][j1].val = ans[i][j];
        f_mx[i1][j1].x = i;
        f_mx[i1][j1].y = j;
      }
    }
  long long sum = -1, x1 = 0, y1 = 0, x2 = 0, y2 = 0;
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 2; j++) {
      if (sum < f_mx[i][j].val + f_mx[!i][j].val) {
        sum = f_mx[i][j].val + f_mx[!i][j].val;
        x1 = f_mx[i][j].x;
        y1 = f_mx[i][j].y;
        x2 = f_mx[!i][j].x;
        y2 = f_mx[!i][j].y;
      }
      if (sum < f_mx[i][j].val + f_mx[i][!j].val) {
        sum = f_mx[i][j].val + f_mx[i][!j].val;
        x1 = f_mx[i][j].x;
        y1 = f_mx[i][j].y;
        x2 = f_mx[i][!j].x;
        y2 = f_mx[i][!j].y;
      }
    }
  cout << sum << "\n"
       << x1 + 1 << " " << y1 + 1 << " " << x2 + 1 << " " << y2 + 1;
  return 0;
}
