#include <iostream>
#ifdef MBP
#include "stdc++.h"
#else
#include <bits/stdc++.h>
#endif

#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

#define REP(i, n) for (long long int i = 0, i##_len = (n); i < i##_len; ++i)
#define REP_R(i, n) for (long long int i = (n - 1); i >= 0; --i)
#define ALL(x) (x).begin(), (x).end()
#define SZ(x) ((long long int)(x).size())
#define BIT(n) (1LL << (n))
#define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end());

#define INF INT_MAX
#define INFLL LLONG_MAX
#define PB push_back
#define MOD 1000000007LL
#define LSB(x) __builtin_ctzll(x)
#define MSB(x) 64LL - __builtin_clzll(x) - 1LL
#define MSBINT(x) 32LL - __builtin_clz(x)
#define BITCNT(x) __builtin_popcountll(x)

using namespace std;
using uif = uint_fast64_t;
using ll = long long int;
using vec = vector<ll>;
using mat = vector<vector<ll>>;
using pa = pair<ll, ll>;
using tTree = __gnu_pbds::tree<ll, __gnu_pbds::null_type, less<ll>, __gnu_pbds::rb_tree_tag, __gnu_pbds::tree_order_statistics_node_update>;

template <class T>
using PQ = priority_queue<T>;
template <class T>
using PQG = priority_queue<T, vector<T>, greater<T>>;

const ll dx[4] = {1, 0, -1, 0};
const ll dy[4] = {0, 1, 0, -1};
const ll dx8[8] = {1, 1, 0, -1, -1, -1, 0, 1};
const ll dy8[8] = {0, 1, 1, 1, 0, -1, -1, -1};

template <class T>
bool chmax(T &, const T &);
template <class T>
bool chmin(T &, const T &);

ll qp(ll a, ll b);
ll qp(ll a, ll b, ll mo);

ll mod_pow(ll, ll);
ll mod_add(ll, ll);
ll mod_mul(ll, ll);
ll mod_div(ll, ll);
void make_fact(ll);
void make_invfact(ll);
void init_fact();

ll combination(ll, ll);

#define FACSIZE 200002LL

ll invfac[FACSIZE];
ll fac[FACSIZE];

ll h, w;

vector<string> f(pa p1, pa p2)

{
  vector<string> ans;
  auto cur = p1;

  while (cur != p2)
  {

    auto dir = 1 - (cur.first % 2) * 2;
    auto nex = cur;
    if (dir == 1 && cur.second == w - 1)
    {
      nex.first++;
    }
    else if (dir == -1 && cur.second == 0)
    {
      nex.first++;
    }
    else
    {
      nex.second += dir;
    }
    //cout << cur.first + 1 << " " << cur.second + 1 << " " << nex.first + 1 << " " << nex.second + 1 << endl;
    string s;
    vec v{cur.first + 1, cur.second + 1, nex.first + 1, nex.second + 1};
    /*
    for (auto x : v)
    {
      s += to_string(x);
      s += " ";
    }
    */
    s += to_string(v[0]) + ' ' + to_string(v[1]) + ' ' + to_string(v[2]) + ' ' + to_string(v[3]);
    ans.PB(s);

    cur = nex;
  }
  return ans;
}

int main(void)
{

  cin >> h >> w;
  mat A(h, vec(w, 0));
  ll even = 0;
  ll odd = 0;
  queue<pa> pos;

  deque<pa> PosDeq[h];

  REP(i, h)
  {
    REP(j, w)
    {
      ll a;
      cin >> a;
      if (a % 2)
      {
        odd++;
        //pos.push(make_pair(i, j));
        if (i % 2)
        {
          PosDeq[i].push_front(make_pair(i, j));
        }
        else
        {
          PosDeq[i].push_back(make_pair(i, j));
        }
      }
      else
      {
        even++;
      }
    }
  }
  //cout << even + odd / 2 << endl;

  REP(i, h)
  {
    //while (!empty(PosDeq[i]))
    while (!PosDeq[i].empty())
    {
      auto p = PosDeq[i].front();
      PosDeq[i].pop_front();
      pos.push(p);
    }
  }
  vector<vector<string>> ans;

  while (pos.size() > 1)
  {
    auto p1 = pos.front();
    pos.pop();
    auto p2 = pos.front();
    pos.pop();
    //cout << p1.first << "," << p1.second << endl;
    //cout << p2.first << "," << p2.second << endl;
    auto v = f(p1, p2);
    ans.PB(v);
  }
  ll cnt = 0;
  for (auto v : ans)
    for (auto s : v)
      cnt++;
  cout << cnt << endl;
  for (auto v : ans)
    for (auto s : v)
      cout << s << endl;

  return 0;
}

template <class T>
bool chmax(T &a, const T &b)
{
  if (a < b)
  {
    a = b;
    return 1;
  }
  return 0;
}
template <class T>
bool chmin(T &a, const T &b)
{
  if (b < a)
  {
    a = b;
    return 1;
  }
  return 0;
}

ll qp(ll a, ll b)
{
  ll ans = 1LL;
  do
  {
    if (b & 1LL)
      ans = 1LL * mod_mul(ans, a) % MOD;
    a = 1LL * mod_mul(a, a) % MOD;
  } while (b >>= 1LL);
  return ans;
}
ll qp(ll a, ll b, ll mo)
{
  ll ans = 1LL;
  do
  {
    if (b & 1LL)
      ans = 1LL * (ans % mo) * (a % mo);
    a = 1LL * (a % mo) * (a % mo);
  } while (b >>= 1LL);
  return ans;
}

ll mod_pow(ll x, ll r)
{
  if (r == 0)
    return 1;
  else if (r == 1)
    return x % MOD;
  else if (r % 2 == 0)
  {
    ll t = mod_pow(x, r / 2) % MOD;
    return t * t % MOD;
  }
  else
  {
    ll t = mod_pow(x, r / 2) % MOD;
    ll k = t * t % MOD;
    return k * x % MOD;
  }
}
ll mod_add(ll a, ll b)
{
  return (a + b) % MOD;
}

ll mod_mul(ll a, ll b)
{
  return (a * b) % MOD;
}

ll mod_div(ll a, ll b)
{
  return mod_mul(a, mod_pow(b, MOD - 2));
}

void make_fact(ll n)
{
  fac[0] = 1;
  REP(i, n)
  {
    fac[i + 1] = mod_mul(fac[i], i + 1);
  }
}

void make_invfact(ll n)
{
  invfac[n] = mod_pow(fac[n], MOD - 2);
  REP_R(i, n)
  {
    invfac[i] = mod_mul(invfac[i + 1], i + 1);
  }
}

void init_fact()
{
  make_fact(FACSIZE - 1);
  make_invfact(FACSIZE - 1);
}

ll combination(ll n, ll c)
{

  return mod_mul(fac[n], mod_mul(invfac[c], invfac[n - c]));
}
///////////////fact
/*
  init_fact();
  cout << fac[5] << endl;
  cout << invfac[5] << endl;
  cout << combination(8, 4) << endl;
  cout << combination(31, 8) << endl;
  //*/
///////////////////

Measured memory_bytes: 17324