#include <bits/stdc++.h>

using namespace std;

int main() {
  int N, K, L[500];
  bool t[60][60];

  cin >> N >> K;
  for(int i = 0; i < K; i++) {
    cin >> L[i];
  }
  int pv = 0;
  for(int i = 0; i < N; i++) {
    string s;
    cin >> s;
    for(int j = 0; j < N; j++) {
      t[i][j] = s[j] == '0';
      pv += t[i][j];
    }
  }

  vector< int > length[26];
  for(int i = 0; i < K; i++) {
    length[L[i]].push_back(i);
  }
  int A[500], B[500], C[500], D[500];

  int ret = 0;

  for(int i = 25; i >= 0; i--) {
    for(int j : length[i]) {
      int best = -55;
      int a, b, c;

      for(int x = 0; x < N; x++) {
        for(int y = 0; y <= N - i; y++) {
          int add = 0;
          for(int k = 0; k < i; k++) {
            add += t[y + k][x] ? -1 : 1;
          }
          if(best < add) {
            best = add;
            a = y, b = x, c = 0;
          }
        }
      }
      for(int y = 0; y < N; y++) {
        for(int x = 0; x <= N - i; x++) {
          int add = 0;
          for(int k = 0; k < i; k++) {
            add += t[y][x + k] ? -1 : 1;
          }
          if(best < add) {
            best = add;
            a = y, b = x, c = 1;
          }
        }
      }

      ret += best;

      if(c == 0) {
        for(int k = 0; k < i; k++) t[a + k][b] ^= 1;
        A[j] = a, B[j] = b, C[j] = a + i - 1, D[j] = b;
      } else {
        for(int k = 0; k < i; k++) t[a][b + k] ^= 1;
        A[j] = a, B[j] = b, C[j] = a, D[j] = b + i - 1;
      }
    }
  }

  cout << ret << endl;

  for(int i = 0; i < K; i++) {
    cout << A[i] + 1 << " " << B[i] + 1 << " " << C[i] + 1 << " " << D[i] + 1 << endl;
  }
}