Buggy python submission (verdict: Runtime Error):
from collections import defaultdict, deque, Counter
from heapq import heappush, heappop, heapify
import math
import bisect
import random
from itertools import permutations, accumulate, combinations
import sys
import string
from bisect import bisect_left, bisect_right
from math import factorial, ceil, floor


INF = float('inf')
def LI(): return list(map(int, sys.stdin.readline().split()))
def I(): return int(sys.stdin.readline())
def LS(): return sys.stdin.readline().split()
def S(): return sys.stdin.readline().strip()
def IR(n): return [I() for i in range(n)]
def LIR(n): return [LI() for i in range(n)]
def SR(n): return [S() for i in range(n)]
def LSR(n): return [LS() for i in range(n)]
def SRL(n): return [list(S()) for i in range(n)]
def MSRL(n): return [[int(j) for j in list(S())] for i in range(n)]
mod = 1000000007


h, w = LI()
L = LIR(h)
ans = []
for y in range(h):
    flag = False
    for x in range(w-1):
        if L[y][x] % 2:
            flag ^= True
        if flag:
            ans += [[y+1, x+1, y+1, x+2]]
    L[y][w-1] += flag



for y in range(h-1):
    flag = False
    if L[y][x] % 2:
        flag ^= True
    if flag:
        ans += [[y + 1, h-1, y + 2, h-1]]



print(len(ans))
for i in ans:
    print(*i)

Fixed submission:
from collections import defaultdict, deque, Counter
from heapq import heappush, heappop, heapify
import math
import bisect
import random
from itertools import permutations, accumulate, combinations
import sys
import string
from bisect import bisect_left, bisect_right
from math import factorial, ceil, floor


INF = float('inf')
def LI(): return list(map(int, sys.stdin.readline().split()))
def I(): return int(sys.stdin.readline())
def LS(): return sys.stdin.readline().split()
def S(): return sys.stdin.readline().strip()
def IR(n): return [I() for i in range(n)]
def LIR(n): return [LI() for i in range(n)]
def SR(n): return [S() for i in range(n)]
def LSR(n): return [LS() for i in range(n)]
def SRL(n): return [list(S()) for i in range(n)]
def MSRL(n): return [[int(j) for j in list(S())] for i in range(n)]
mod = 1000000007


h, w = LI()
L = LIR(h)
ans = []
for y in range(h):
    flag = False
    for x in range(w-1):
        if L[y][x] % 2:
            flag ^= True
        if flag:
            ans += [[y+1, x+1, y+1, x+2]]
    L[y][w-1] += flag


flag = False
for y in range(h-1):
    if L[y][w-1] % 2:
        flag ^= True
    if flag:
        ans += [[y+1, w, y+2, w]]



print(len(ans))
for i in ans:
    print(*i)

Buggy python submission (verdict: Runtime Error):
from collections import defaultdict, deque, Counter
from heapq import heappush, heappop, heapify
import math
import bisect
import random
from itertools import permutations, accumulate, combinations
import sys
import string
from bisect import bisect_left, bisect_right
from math import factorial, ceil, floor


INF = float('inf')
def LI(): return list(map(int, sys.stdin.readline().split()))
def I(): return int(sys.stdin.readline())
def LS(): return sys.stdin.readline().split()
def S(): return sys.stdin.readline().strip()
def IR(n): return [I() for i in range(n)]
def LIR(n): return [LI() for i in range(n)]
def SR(n): return [S() for i in range(n)]
def LSR(n): return [LS() for i in range(n)]
def SRL(n): return [list(S()) for i in range(n)]
def MSRL(n): return [[int(j) for j in list(S())] for i in range(n)]
mod = 1000000007


h, w = LI()
L = LIR(h)
for y in range(h):
    flag = False
    for x in range(w-1):
        if L[y][x] % 2:
            flag ^= True
        if flag:
            print(y+1, x+1, y+1, x+2)
    L[y][w-1] += flag



for y in range(h-1):
    flag = False
    if L[y][x] % 2:
        flag ^= True
    if flag:
        print(y + 1, h-1, y + 2, h-1)

Fixed submission:
from collections import defaultdict, deque, Counter
from heapq import heappush, heappop, heapify
import math
import bisect
import random
from itertools import permutations, accumulate, combinations
import sys
import string
from bisect import bisect_left, bisect_right
from math import factorial, ceil, floor


INF = float('inf')
def LI(): return list(map(int, sys.stdin.readline().split()))
def I(): return int(sys.stdin.readline())
def LS(): return sys.stdin.readline().split()
def S(): return sys.stdin.readline().strip()
def IR(n): return [I() for i in range(n)]
def LIR(n): return [LI() for i in range(n)]
def SR(n): return [S() for i in range(n)]
def LSR(n): return [LS() for i in range(n)]
def SRL(n): return [list(S()) for i in range(n)]
def MSRL(n): return [[int(j) for j in list(S())] for i in range(n)]
mod = 1000000007


h, w = LI()
L = LIR(h)
ans = []
for y in range(h):
    flag = False
    for x in range(w-1):
        if L[y][x] % 2:
            flag ^= True
        if flag:
            ans += [[y+1, x+1, y+1, x+2]]
    L[y][w-1] += flag


flag = False
for y in range(h-1):
    if L[y][w-1] % 2:
        flag ^= True
    if flag:
        ans += [[y+1, w, y+2, w]]



print(len(ans))
for i in ans:
    print(*i)