#include<stdio.h>
#include<stdlib.h>
#include<math.h>

typedef long long int int64;

#define MAX(a,b) ((a)>(b)?(a):(b))
#define MIN(a,b) ((a)<(b)?(a):(b))
#define ABS(a) ((a)>(0)?(a):-(a))

#define POS(i,j) ((i)*(n+1)+(j))

typedef struct length{
  int l;
  int index;
} len;

int cmp(const void *a,const void *b){
  return -(((len *)a)->l-((len *)b)->l);
}

void flip(int l,char *s,int n,int *t){
  int a=t[0];
  int b=t[1];
  if(a==t[2]){
    int i;
    for(i=0;i<l;i++){
      int index=POS(a,b+i);
      s[index]=s[index]=='0'?'1':'0';
    }
  } else {
    int i;
    for(i=0;i<l;i++){
      int index=POS(a+i,b);
      s[index]=s[index]=='0'?'1':'0';
    }
  }
  return;
}

void func(int l,char *s,const int n,int *a){
  int max=-l-1;
  int i,j;
  for(j=0;j<n;j++){
    int c=0;
    for(i=0;i<l;i++){
      c+=(s[POS(i,j)]=='0'?-1:1);
    }
    i=0;
    if(c>max){
      max=c;
      a[0]=i;
      a[1]=j;
      a[2]=i+l-1;
      a[3]=j;
    }
    while(i+l<n){
      c-=(s[POS(i,j)]=='0'?-1:1);
      c+=(s[POS(i+l,j)]=='0'?-1:1);
      i++;
      if(c>max){
	max=c;
	a[0]=i;
	a[1]=j;
	a[2]=i+l-1;
	a[3]=j;
      }
    }
  }
  for(i=0;i<n;i++){
    int c=0;
    for(j=0;j<l;j++){
      c+=(s[POS(i,j)]=='0'?-1:1);
    }
    j=0;
    if(c>max){
      max=c;
      a[0]=i;
      a[1]=j;
      a[2]=i;
      a[3]=j+l-1;
    }
    while(j+l<n){
      c-=(s[POS(i,j  )]=='0'?-1:1);
      c+=(s[POS(i,j+l)]=='0'?-1:1);
      j++;
      if(c>max){
	max=c;
	a[0]=i;
	a[1]=j;
	a[2]=i;
	a[3]=j+l-1;
      }
    }
  }
  flip(l,s,n,a);
  return;
}

int score(const char *s,const int n){
  int res=0;
  int i,j;
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      res+=(s[POS(i,j)]=='0'?1:0);
    }
  }
  return res;
}

void algo1(char *s,int n,len *array,int k,int *ans,int *tmp){
  const int parameter=4;
  int i;
  for(i=0;i<k && array[i].l>=parameter;i++){
    int l=array[i].l;
    int index=array[i].index;
    func(l,s,n,ans+4*index);
  }
  int maxScore=score(s,n);
  const int m=i;
  int iter=10;
  while(iter--){
    for(i=0;i<4*k;i++) tmp[i]=ans[i];
    for(i=0;i<m;i++){
      int p=rand()%m;
      flip(array[p].l,s,n,tmp+4*array[p].index);
      func(array[p].l,s,n,tmp+4*array[p].index);
    }
    for(i=0;i<m;i++){
      int l=array[i].l;
      int index=array[i].index;
      flip(l,s,n,tmp+4*index);
      func(l,s,n,tmp+4*index);
    }
    int localScore=score(s,n);
    if(localScore>maxScore){
      maxScore=localScore;
      for(i=0;i<4*k;i++) ans[i]=tmp[i];
    } else {
      for(i=0;i<m;i++){
	flip(array[i].l,s,n,tmp+4*array[i].index);
	flip(array[i].l,s,n,ans+4*array[i].index);
      }
    }
  }
  for(i=m;i<k;i++){
    int l=array[i].l;
    int index=array[i].index;
    func(l,s,n,ans+4*index);
  }
  return;
}

void swap(len *array,int x,int y){
  len t=array[x];
  array[x]=array[y];
  array[y]=t;
  return;
}

void algo2(char *s,int n,len *array,int k,int *ans,int *tmp){
  const int parameter=5;
  int i;
  for(i=0;i<k && array[i].l>=parameter;i++){
    int l=array[i].l;
    int index=array[i].index;
    func(l,s,n,ans+4*index);
  }
  int maxScore=score(s,n);
  const int m=i;
  int iter=20;
  while(iter--){
    for(i=0;i<4*k;i++) tmp[i]=ans[i];
    for(i=0;i<m;i++) swap(array,i,rand()%m);
    for(i=0;i<m;i++) flip(array[i].l,s,n,tmp+4*array[i].index);
    for(i=0;i<m;i++) func(array[i].l,s,n,tmp+4*array[i].index);
    int localScore=score(s,n);
    if(localScore>maxScore){
      maxScore=localScore;
      for(i=0;i<4*k;i++) ans[i]=tmp[i];
    } else {
      for(i=0;i<m;i++){
	flip(array[i].l,s,n,tmp+4*array[i].index);
	flip(array[i].l,s,n,ans+4*array[i].index);
      }
    }
  }
  for(i=m;i<k;i++){
    int l=array[i].l;
    int index=array[i].index;
    func(l,s,n,ans+4*index);
  }
  return;
}

void run(void){
  int n,k;
  scanf("%d%d",&n,&k);
  char *s=(char *)calloc((n+1)*n,sizeof(char));
  len *array=(len *)malloc(sizeof(len)*k);
  int *ans=(int *)malloc(sizeof(int)*4*k);
  int *tmp=(int *)malloc(sizeof(int)*4*k);
  int i;
  for(i=0;i<k;i++){
    int t;
    scanf("%d",&t);
    array[i].l=t;
    array[i].index=i;
  }
  qsort(array,k,sizeof(len),cmp);
  for(i=0;i<n;i++) scanf("%s",s+POS(i,0));
  algo2(s,n,array,k,ans,tmp);
  for(i=0;i<k;i++){
    printf("%d %d %d %d\n",ans[i*4+0]+1,ans[i*4+1]+1,ans[i*4+2]+1,ans[i*4+3]+1);
  }
  return;
}

int main(void){
  run();
  return 0;
}