#include <cstdio>
using namespace std;

int n,k;
int lc0[26],lc[26],lcs[26],li[26][500];
char a0[64][64],aa[64][64],as[64][64];
int ra[500],rb[500],rc[500],rd[500];
int ra1[500],rb1[500],rc1[500],rd1[500];
int score0,score,score1;
int nazo1;
int nazo2;

void reg(int l,int a,int c,int b, int d){
	int i=li[l][--lc[l]];
	ra[i]=a;
	rb[i]=b;
	rc[i]=c;
	rd[i]=d;
	for(int y=a;y<=c;++y){
		for(int x=b;x<=d;++x){
			aa[y][x]^=1;
		}
	}
}

void solve(int ux,int uy,int vx,int vy){
	{
		int ov=vx+vy<0?n+1:0;
		int ou=ux+uy<0?n+1:0;
		for(int y=1;y<=n;++y){
			for(int x=1;x<=n;++x){
				aa[y][x]=a0[ov+vx*x+vy*y][ou+ux*x+uy*y];
			}
		}
	}
	for(int i=1;i<=25;++i){
		lc[i]=lc0[i];
	}

	for(int j=0;j<2;++j){
		for(int y=1;y<=n;){
			int bl=0,bx=0;
			int bv=435-nazo1;
			for(int x=1;x<=n;++x){
				int v=x*-2-y;
				for(int l=0;l<25&&x+l<=n;++l){
					v+=aa[y][x+l]!=aa[y-1][x+l]?100:-nazo1-200;
					v+=aa[y][x+l]!=aa[y+1][x+l]?100:nazo1;
					if(bv<=v&&lc[l+1]){
						bv=v; bx=x; bl=l+1;
					}
				}
			}
			if(bl){
				reg(bl,y,y,bx,bx+bl-1);
			}else{
				++y;
			}
		}
	}

	bool first=true;
	for(nazo2=183;nazo2<=206;++nazo2){
		if(first){
			for(int i=1;i<=25;++i){
				lcs[i]=lc[i];
			}
			for(int y=1;y<=n;++y){
				for(int x=1;x<=n;++x){
					as[y][x]=aa[y][x];
				}
			}
		}else{
			for(int i=1;i<=25;++i){
				lc[i]=lcs[i];
			}
			for(int y=1;y<=n;++y){
				for(int x=1;x<=n;++x){
					aa[y][x]=as[y][x];
				}
			}
		}
		first=false;
	
		for(int y=1;y<=n;++y){
			for(int x=1;x<=n;++x) {
				if(aa[y][x]&(y==n|aa[y+1][x])){
					int bl=0,bv=0,v=0;
					for(int l=0;l<25&&y+l<=n;++l){
						v+=aa[y+l][x]*6-1;
						int d=aa[y+l][x]-aa[y+l+1][x];
						int v1=v+(d>0?200:d<0?nazo2:0);
						if(bv<=v1&&lc[l+1]){
							bv=v1;
							bl=l+1;
						}
					}
					if(bl){
						reg(bl,y,y+bl-1,x,x);
					}
				}
			}
		}

		for(int bl=25;bl>=1;--bl){
			while(lc[bl]){
				int by=1,bx=1,bv=0;
				for(int y=1;y<=n;++y){
					int v=0;
					for(int x=1;x<=bl;++x){
						v+=aa[y][x];
					}
					for(int x=1;x+bl-1<=n;++x){
						if(bv<v){
							bv=v; by=y; bx=x;
						}
						v+=aa[y][x+bl]-aa[y][x];
					}
				}
				reg(bl,by,by,bx,bx+bl-1);
			}
		}

		score=score0;
		for(int y=1;y<=n;++y){
			for(int x=1;x<=n;++x){
				score-=aa[y][x];
			}
		}
		if(score1<score){
			score1=score;
			for(int i=0;i<k;++i){
				ra1[i]=vy>0?ra[i]:vy<0?n+1-rc[i]:vx>0?rb[i]:vx<0?n+1-rd[i]:0;
				rb1[i]=uy>0?ra[i]:uy<0?n+1-rc[i]:ux>0?rb[i]:ux<0?n+1-rd[i]:0;
				rc1[i]=vy>0?rc[i]:vy<0?n+1-ra[i]:vx>0?rd[i]:vx<0?n+1-rb[i]:0;
				rd1[i]=uy>0?rc[i]:uy<0?n+1-ra[i]:ux>0?rd[i]:ux<0?n+1-rb[i]:0;
			}
		}
	}
}

main(){
	scanf("%d%d",&n,&k);
	for(int i=0;i<k;++i){
		int l;
		scanf("%d ",&l);
		li[l][lc0[l]++]=i;
	}
	for(int y=1;y<=n;++y){
		gets(a0[y]+1);
		for(int x=1;x<=n;++x){
			score0+=a0[y][x]&=1;
		}
	}
	for(nazo1=-55;nazo1<=9;nazo1+=2){
		solve( 1, 0, 0, 1);
		solve(-1, 0, 0, 1);
		solve( 1, 0, 0,-1);
		solve(-1, 0, 0,-1);
		solve( 0, 1, 1, 0);
		solve( 0,-1, 1, 0);
		solve( 0, 1,-1, 0);
		solve( 0,-1,-1, 0);
	}
	for(int i=0;i<k;++i){
		printf("%d %d %d %d\n",ra1[i],rb1[i],rc1[i],rd1[i]);
	}
}