#include<iostream>
#include<map>
#include<set>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<limits>
#include<sstream> 
#include<unordered_map>      

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef int itn;
typedef vector<int> vi;
typedef vector<ll> vl;

const int inf = numeric_limits<int>::max() >> 2;
const ll linf = numeric_limits<ll>::max() >> 2;
const ull ulinf = numeric_limits<ull>::max() >> 2;
const double pi = acos(-1);
const int dx[4]={0,-1,1,0};
const int dy[4]={-1,0,0,1};
const int dx8[8]={-1,0,1,-1,1,-1,0,1};
const int dy8[8]={-1,-1,-1,0,0,1,1,1};

#define p_queue(i) priority_queue<i> 
#define rp_queue(i) priority_queue<i, vector<i>, greater<i>> 
#define umap(i,j) unordered_map<i,j>
#define P(p) cout<<(p)<<endl
#define PS(p) cout<<(p)<< " "
#define IN cin >> 
#define rep(i,m,n) for(int i = (m); i < (int)(n); i++)
#define rrep(i,m,n) for(int i = (m-1); i >= (int)(n); i--)
#define inrep(n,a) for(int i = 0; i < (int)(n); i++)cin >> a[i];
#define mod(i) ((i)%(ll)(1e9 + 7))
#define divm(a,b) (mod(a * modpow((ll)b,(ll)(1e9+5))))
#define rsort(a,b,c) sort(a, b, greater<c>())
#define vsort(v) sort(v.begin(), v.end());
#define rvsort(v,c) sort(v.begin(), v.end(), greater<c>());
#define ft first
#define sd second
#define pb push_back
#define it insert
#define sz(x) ((int)(x).size())
#define lb(a,n,k) (lower_bound(a,a+n,k) - a) 
#define vlb(a,k)  (lower_bound(a.begin(),a.end(),k) - a.begin())
#define ub(a,n,k) (upper_bound(a,a+n,k) - a) 
#define vub(a,k)  (upper_bound(a.begin(),a.end(),k) - a.begin())
#define YES cout<<"YES"<< endl
#define NO cout<<"NO"<<endl
#define Yes cout<<"Yes"<<endl
#define No cout<<"No"<<endl  
#define yes cout<<"yes"<<endl
#define no cout<<"no"<<endl
#define ret return

ll modpow(ll i,ll j){ ll tmp=1; while(j){ if(j%2)tmp=mod(tmp*i);i=mod(i*i);j/=2;}return tmp;}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll lcm(ll a, ll b) { return a / gcd(a, b) * b; }

vector<string> split(const string &str, char sep){
    vector<string> v;
    stringstream ss(str);
    string buffer;
    while( getline(ss, buffer, sep) ) {
        v.push_back(buffer);
    }
    return v;
}

//ll ncr[100][100];
//ll nCr(int n, int r){if(n==r) ret ncr[n][r] = 1; if(r==0) ret ncr[n][r] = 1; if(r==1) ret ncr[n][r] = n;if(ncr[n][r]) ret ncr[n][r]; ret ncr[n][r] = nCr(n-1,r) + nCr(n-1,r-1);}

//ll npr[100][100]={};
//ll nPr(int n,int r){if(npr[n][r])ret npr[n][r];if(r==0)ret npr[n][r] = 1;if(r==1)ret npr[n][r] = n;ret npr[n][r] = n * nPr(n-1,r-1);}

//ll nHr(int n,int r){ret nCr(n+r-1,r);}

///////////////////////////////////////////////////////////////////////////
int h,w;
int a[505][505] = {};
vector<int> x,xx,y,yy;

int main(){
    cin >> h >> w;
    rep(i,0,h)rep(j,0,w)cin >> a[i][j];


    rep(i,0,h){
        rep(j,0,w){
            if(a[i][j] > 0 && a[i][j]&1){
                int f = 0;
                rep(k,2,4){
                    int nx = j+dx[k],ny = i+dy[k];
                    if(0 <= nx && nx < w && 0 <= ny && ny < h ){
                        
                        if(a[ny][nx]&1){
                            a[ny][nx]++;
                            a[i][j] --;

                            x.pb(j);
                            y.pb(i);
                            xx.pb(nx);
                            yy.pb(ny);
                            f = 1;
                            break;
                        }
                    }
                }
                if(f == 0){
                    if(i < h-1){
                        a[i][j]--;
                        a[i+1][j]++;

                        x.pb(j);
                        y.pb(i);
                        xx.pb(j);
                        yy.pb(i+1);
                    }else if(j < w-1){
                        a[i][j]--;
                        a[i][j+1]++;

                        x.pb(j);
                        y.pb(i);
                        xx.pb(j+1);
                        yy.pb(i);
                    }else{

                    }
                }
            }
        }
    }

    P(sz(x));
    rep(i,0,sz(x))cout << y[i]+1 << " " << x[i]+1 << " " << yy[i]+1 << " " << xx[i]+1 <<endl;
    
   /* rep(i,0,h){
        rep(j,0,w)PS(a[i][j]);
        cout << endl;
    }*/
	return  0;
}

Measured memory_bytes: 5220