#include<cstdio>
#include<bitset>
using namespace std;
const int o=2010,O=o*o*2;
int n,m,cnt[o],I[O],J[O],typ[O],K,ch[O][2],tot,h[o],Cnt,st[o][o],tp[o],v[o],t;char s[o][o];bool flg;bitset<o> b[o][2];
struct Edge{int v,p;}e[O*2];
inline void ad(int U,int V){e[++Cnt].v=V;e[Cnt].p=h[U];h[U]=Cnt;}
bool Dfs(int nw,int val){
	if(v[nw]) return v[nw]==val;
	v[st[t][++tp[t]]=nw]=val;
	for(int i=h[nw],j;i;i=e[i].p){
		j=e[i].v;
		if(nw==I[j]&&val-1-((typ[j]-1)&1)) if(!Dfs(J[j],((typ[j]-1)>>1)+1)) return false;
		if(nw==J[j]&&((typ[j]-1)>>1)+1) if(!Dfs(I[j],val-1-((typ[j]-1)&1))) return false;
	}
	return true;
}
void dfs(int nw,int u){
	if(flg) return;
	if(nw>1&&!u){flg=1;printf("-1");return;}
	if(nw>m) return;
	if(cnt[nw]==n){dfs(nw+1,ch[u][1]);return;}
	if(cnt[nw]==0){dfs(nw+1,ch[u][0]);return;}
	if(!v[nw]){
		t=nw;
		if(Dfs(nw,1)) dfs(nw+1,ch[u][0]);
		for(;tp[nw];v[st[nw][tp[nw]--]]=0);
		if(Dfs(nw,2)) dfs(nw+1,ch[u][1]);
		for(;tp[nw];v[st[nw][tp[nw]--]]=0);
	}
	else dfs(nw+1,ch[u][v[nw]-1]);
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;++i) scanf("%s",s[i]+1);
	for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) cnt[j]+=(s[i][j]=='1');
	for(int i=1;i<=m;++i)
		if(cnt[i]==n) I[++K]=i,J[K]=i,typ[K]=4;
		else if(!cnt[i]) I[++K]=i,J[K]=i,typ[K]=1;
	for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) b[j][s[i][j]-'0'].set(i);
	for(int i=1;i<=m;++i) for(int j=i+1;j<=m;++j) if(cnt[i]&&cnt[i]<n&&cnt[j]&&cnt[j]<n){
		if(!(b[i][1]&b[j][1]).count()) I[++K]=i,J[K]=j,typ[K]=1,ad(i,K),ad(j,K);
		if(!(b[i][0]&b[j][1]).count()) I[++K]=i,J[K]=j,typ[K]=2,ad(i,K),ad(j,K);
		if(!(b[i][1]&b[j][0]).count()) I[++K]=i,J[K]=j,typ[K]=3,ad(i,K),ad(j,K);
		if(!(b[i][0]&b[j][0]).count()) I[++K]=i,J[K]=j,typ[K]=4,ad(i,K),ad(j,K);
	}
	for(int i=1;i<=n;++i) for(int u=0,j=1;j<=m;u=ch[u][s[i][j++]-48])
		if(!ch[u][s[i][j]-48]) ch[u][s[i][j]-48]=++tot;
	dfs(1,0);
	if(flg) return 0;
	printf("%d\n",K);
	for(int i=1;i<=K;++i) printf("%d %d %d\n",I[i]-1,J[i]-1,typ[i]);
	return 0;
}