#include <bits/stdc++.h>

#define REP(i,a,b) for(int i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

using namespace std;

int INF = 1<<29;

int a[61][61];
int l[26];
int lsum[66][66];
int N, K;

void calc_xlsum(int i) {
  for (int j = N - 1; j >= 0; j--) {
    lsum[i][j] = (a[i][j] ? 1 : -1) + (j < N - 1 ? lsum[i][j + 1] : 0);
  }
}

void solve() {
//  int raw = 0;
//  rep(i, N) rep(j, N) raw += !a[i][j];
  vector<tuple<int, int, int, int>> vec;
//  int ret = raw;
  rep(k, K) {
    int mx = -INF;
    int idx_i = -1, idx_j = -1;
    rep(i, N) rep(j, N - l[k] + 1) {
      int cand = lsum[i][j] - ( j + l[k] < N ? lsum[i][j + l[k]] : 0 );
      if (mx < cand) {
        mx = cand;
        idx_i = i, idx_j = j;
      }
    }
//    ret += mx;
    vec.emplace_back(idx_i + 1, idx_j + 1, idx_i + 1, idx_j + l[k]);
    REP(j, idx_j, idx_j + l[k] + 1) {
      a[idx_i][j] ^= 1;
    }
    /*
    printf("-------------------\n");
    rep(i, N) {
      rep(j, N) {
        printf("%3d", a[i][j]);
      }
      printf("\n");
    }
    printf("===================\n");
    */
    calc_xlsum(idx_i);
  }
  for (auto e: vec) {
    int x, y, z, w; tie(x, y, z, w) = e;
    printf("%d %d %d %d\n", x, y, z, w);
  }
//  printf("ret = %d\n", ret);
//  printf("diff = %d\n", ret - raw);
}

int main() {
  cin >> N >> K;
  rep(i, K) {
    cin >> l[i];
  }
  sort(l, l + K, greater<int>());

  rep(i, N) {
    string s; cin >> s;
    rep(j, N) a[i][j] = s[j] - '0';
  }

  rep(i, N)
    calc_xlsum(i);

  solve();
}