Problem: This is the hard version of the problem. The difference between the versions is in the number of possible operations that can be made. You can make hacks if and only if you solved both versions of the problem.

You are given a binary table of size n × m. This table consists of symbols 0 and 1.

You can make such operation: select 3 different cells that belong to one 2 × 2 square and change the symbols in these cells (change 0 to 1 and 1 to 0).

Your task is to make all symbols in the table equal to 0. You are allowed to make at most nm operations. You don't need to minimize the number of operations.

It can be proved, that it is always possible.

Input

The first line contains a single integer t (1 ≤ t ≤ 5000) — the number of test cases. The next lines contain descriptions of test cases.

The first line of the description of each test case contains two integers n, m (2 ≤ n, m ≤ 100).

Each of the next n lines contains a binary string of length m, describing the symbols of the next row of the table.

It is guaranteed, that the sum of nm for all test cases does not exceed 20000.

Output

For each test case print the integer k (0 ≤ k ≤ nm) — the number of operations.

In the each of the next k lines print 6 integers x_1, y_1, x_2, y_2, x_3, y_3 (1 ≤ x_1, x_2, x_3 ≤ n, 1 ≤ y_1, y_2, y_3 ≤ m) describing the next operation. This operation will be made with three cells (x_1, y_1), (x_2, y_2), (x_3, y_3). These three cells should be different. These three cells should belong to some 2 × 2 square.

Example

Input


5
2 2
10
11
3 3
011
101
110
4 4
1111
0110
0110
1111
5 5
01011
11001
00010
11011
10000
2 3
011
101


Output


1
1 1 2 1 2 2
2 
2 1 3 1 3 2
1 2 1 3 2 3
4
1 1 1 2 2 2 
1 3 1 4 2 3
3 2 4 1 4 2
3 3 4 3 4 4
4
1 2 2 1 2 2 
1 4 1 5 2 5 
4 1 4 2 5 1
4 4 4 5 3 4
2
1 3 2 2 2 3
1 2 2 1 2 2

Note

In the first test case, it is possible to make only one operation with cells (1, 1), (2, 1), (2, 2). After that, all symbols will be equal to 0.

In the second test case:

  * operation with cells (2, 1), (3, 1), (3, 2). After it the table will be: 
    
          
    011  
    001  
    000  
    

  * operation with cells (1, 2), (1, 3), (2, 3). After it the table will be: 
    
          
    000  
    000  
    000  
    




In the fifth test case:

  * operation with cells (1, 3), (2, 2), (2, 3). After it the table will be: 
    
          
    010  
    110  
    

  * operation with cells (1, 2), (2, 1), (2, 2). After it the table will be: 
    
          
    000  
    000  
    

Solution (in PYTHON):

Solution:

# Author : raj1307 - Raj Singh
# Date   : 18.11.2020

from __future__ import division, print_function

import os,sys
from io import BytesIO, IOBase

if sys.version_info[0] < 3:
    from __builtin__ import xrange as range
    from future_builtins import ascii, filter, hex, map, oct, zip


def ii():  return int(input())
def si():  return input()
def mi():  return map(int,input().strip().split(" "))
def msi(): return map(str,input().strip().split(" "))
def li():  return list(mi())

def dmain():
    sys.setrecursionlimit(1000000)
    threading.stack_size(1024000)
    thread = threading.Thread(target=main)
    thread.start()
    
#from collections import deque, Counter, OrderedDict,defaultdict
#from heapq import nsmallest, nlargest, heapify,heappop ,heappush, heapreplace
#from math import log,sqrt,factorial,cos,tan,sin,radians
#from bisect import bisect,bisect_left,bisect_right,insort,insort_left,insort_right
#from decimal import *
#import threading
#from itertools import permutations
#Copy 2D list  m = [x[:] for x in mark] .. Avoid Using Deepcopy

abc='abcdefghijklmnopqrstuvwxyz'
abd={'a': 0, 'b': 1, 'c': 2, 'd': 3, 'e': 4, 'f': 5, 'g': 6, 'h': 7, 'i': 8, 'j': 9, 'k': 10, 'l': 11, 'm': 12, 'n': 13, 'o': 14, 'p': 15, 'q': 16, 'r': 17, 's': 18, 't': 19, 'u': 20, 'v': 21, 'w': 22, 'x': 23, 'y': 24, 'z': 25}
mod=1000000007
#mod=998244353
inf = float("inf")
vow=['a','e','i','o','u']
dx,dy=[-1,1,0,0],[0,0,1,-1]
def getKey(item): return item[1] 
def sort2(l):return sorted(l, key=getKey,reverse=True)
def d2(n,m,num):return [[num for x in range(m)] for y in range(n)]
def isPowerOfTwo (x): return (x and (not(x & (x - 1))) )
def decimalToBinary(n): return bin(n).replace("0b","")
def ntl(n):return [int(i) for i in str(n)]
def ncr(n,r): return factorial(n)//(factorial(r)*factorial(max(n-r,1)))

def ceil(x,y):
    if x%y==0:
        return x//y
    else:
        return x//y+1

def powerMod(x,y,p):
    res = 1
    x %= p
    while y > 0:
        if y&1:
            res = (res*x)%p
        y = y>>1
        x = (x*x)%p
    return res

def gcd(x, y):
    while y:
        x, y = y, x % y
    return x
    
def isPrime(n) : # Check Prime Number or not 
    if (n <= 1) : return False
    if (n <= 3) : return True
    if (n % 2 == 0 or n % 3 == 0) : return False
    i = 5
    while(i * i <= n) : 
        if (n % i == 0 or n % (i + 2) == 0) : 
            return False
        i = i + 6
    return True



def read():
    sys.stdin  = open('input.txt', 'r')  
    sys.stdout = open('output.txt', 'w') 





def main():



    for _ in range(ii()):


        n,m=mi()


        a=[]
        for i in range(n):
            x=list(si())
            for j in range(m):
                x[j]=int(x[j])
            a.append(x)


        ans=[]

        f=0
        for i in range(n-1,0,-1):
            for j in range(m-1,0,-1):

                if i==1 and j==1:

                    f=1

                    i=0
                    j=0

                    c=a[i][j]+a[i+1][j]+a[i][j+1]+a[i+1][j+1]
                    

                    if c==0:
                        pass
     
                    elif c==1:
     
                        if a[i][j]==1:    
                            ans.append([[i,j],[i+1,j],[i+1,j+1]])
                            ans.append([[i,j],[i+1,j],[i,j+1]])
                            ans.append([[i,j],[i+1,j+1],[i,j+1]])
     
                        elif a[i+1][j]==1:
     
                            ans.append([[i+1,j],[i,j],[i,j+1]])
                            ans.append([[i+1,j],[i+1,j+1],[i,j]])
                            ans.append([[i+1,j],[i+1,j+1],[i,j+1]])
     
                        elif a[i][j+1]==1:
     
                            ans.append([[i,j],[i+1,j],[i,j+1]])
                            ans.append([[i,j],[i,j+1],[i+1,j+1]])
                            ans.append([[i+1,j],[i,j+1],[i+1,j+1]])
     
                        elif a[i+1][j+1]==1:
     
                            ans.append([[i,j],[i,j+1],[i+1,j+1]])
                            ans.append([[i,j],[i+1,j],[i+1,j+1]])
                            ans.append([[i+1,j],[i,j+1],[i+1,j+1]])
     
     
                    elif c==2:
     
                        if a[i][j]==a[i+1][j+1]:
     
                            if a[i][j]==1:
     
                                ans.append([[i+1,j],[i+1,j+1],[i,j+1]])
                                ans.append([[i,j],[i+1,j],[i,j+1]])
     
                            else:
     
                                ans.append([[i,j+1],[i+1,j+1],[i,j]])
                                ans.append([[i+1,j],[i,j],[i+1,j+1]])
     
                        else:
     
     
                            if a[i+1][j]==1 and a[i+1][j+1]==1:    
                                ans.append([[i,j],[i+1,j],[i,j+1]])
                                ans.append([[i,j],[i+1,j+1],[i,j+1]])
     
                            elif a[i][j]==1 and a[i][j+1]==1:
     
                                ans.append([[i+1,j],[i+1,j+1],[i,j]])
                                ans.append([[i+1,j],[i+1,j+1],[i,j+1]])
     
                            elif a[i][j]==1 and a[i+1][j]==1:
     
                                ans.append([[i,j],[i,j+1],[i+1,j+1]])
                                ans.append([[i+1,j],[i,j+1],[i+1,j+1]])
     
                            else:
     
                                ans.append([[i,j],[i+1,j],[i+1,j+1]])
                                ans.append([[i+1,j],[i,j+1],[i,j]])
     
                    elif c==3:
     
                        if a[i][j]==1 and a[i+1][j]==1 and a[i+1][j+1]==1:
                            ans.append([[i,j],[i+1,j],[i+1,j+1]])
                        elif a[i][j+1]==1 and a[i+1][j]==1 and a[i+1][j+1]==1:
                            ans.append([[i,j+1],[i+1,j],[i+1,j+1]])
     
                        elif a[i][j+1]==1 and a[i][j]==1 and a[i+1][j+1]==1:
                            ans.append([[i,j+1],[i,j],[i+1,j+1]])
                        else:
                            ans.append([[i,j],[i+1,j],[i,j+1]])
     
     
                    else:
     
     
                        ans.append([[i,j+1],[i+1,j],[i+1,j+1]])
                        ans.append([[i,j],[i+1,j],[i+1,j+1]])
                        ans.append([[i,j],[i+1,j],[i,j+1]])
                        ans.append([[i,j],[i+1,j+1],[i,j+1]])



                    break


                elif i==1:

                    if a[i][j]==1 and a[i-1][j]==1:
                        ans.append([[i,j],[i-1,j],[i,j-1]])
                        a[i][j]^=1
                        a[i-1][j]^=1
                        a[i][j-1]^=1
                    elif a[i][j]==1:
                        ans.append([[i,j],[i,j-1],[i-1,j-1]])
                        a[i][j]^=1
                        a[i-1][j-1]^=1
                        a[i][j-1]^=1
                    elif a[i-1][j]==1:
                        ans.append([[i-1,j],[i-1,j-1],[i,j-1]])
                        a[i-1][j]^=1
                        a[i-1][j-1]^=1
                        a[i][j-1]^=1
                    else:
                        pass



                elif j==1:

                    if a[i][j]==1 and a[i][j-1]==1:
                        ans.append([[i,j],[i-1,j],[i,j-1]])
                        a[i][j]^=1
                        a[i-1][j]^=1
                        a[i][j-1]^=1
                    elif a[i][j]==1:
                        ans.append([[i,j],[i-1,j-1],[i-1,j]])
                        a[i][j]^=1
                        a[i-1][j-1]^=1
                        a[i-1][j]^=1
                    elif a[i][j-1]==1:
                        ans.append([[i-1,j],[i-1,j-1],[i,j-1]])
                        a[i-1][j]^=1
                        a[i-1][j-1]^=1
                        a[i][j-1]^=1
                    else:
                        pass

                else:

                    if a[i][j]==1:
                        ans.append([[i,j],[i-1,j],[i,j-1]])
                        a[i][j]^=1
                        a[i-1][j]^=1
                        a[i][j-1]^=1
                    else:
                        pass

            if f:
                break




        print(len(ans))

        for i in ans:
            print(i[0][0]+1,i[0][1]+1,i[1][0]+1,i[1][1]+1,i[2][0]+1,i[2][1]+1)

















# region fastio
# template taken from https://github.com/cheran-senthil/PyRival/blob/master/templates/template.py

BUFSIZE = 8192


class FastIO(IOBase):
    newlines = 0

    def __init__(self, file):
        self._fd = file.fileno()
        self.buffer = BytesIO()
        self.writable = "x" in file.mode or "r" not in file.mode
        self.write = self.buffer.write if self.writable else None

    def read(self):
        while True:
            b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
            if not b:
                break
            ptr = self.buffer.tell()
            self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
        self.newlines = 0
        return self.buffer.read()

    def readline(self):
        while self.newlines == 0:
            b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
            self.newlines = b.count(b"\n") + (not b)
            ptr = self.buffer.tell()
            self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
        self.newlines -= 1
        return self.buffer.readline()

    def flush(self):
        if self.writable:
            os.write(self._fd, self.buffer.getvalue())
            self.buffer.truncate(0), self.buffer.seek(0)


class IOWrapper(IOBase):
    def __init__(self, file):
        self.buffer = FastIO(file)
        self.flush = self.buffer.flush
        self.writable = self.buffer.writable
        self.write = lambda s: self.buffer.write(s.encode("ascii"))
        self.read = lambda: self.buffer.read().decode("ascii")
        self.readline = lambda: self.buffer.readline().decode("ascii")


def print(*args, **kwargs):
    """Prints the values to a stream, or to sys.stdout by default."""
    sep, file = kwargs.pop("sep", " "), kwargs.pop("file", sys.stdout)
    at_start = True
    for x in args:
        if not at_start:
            file.write(sep)
        file.write(str(x))
        at_start = False
    file.write(kwargs.pop("end", "\n"))
    if kwargs.pop("flush", False):
        file.flush()


if sys.version_info[0] < 3:
    sys.stdin, sys.stdout = FastIO(sys.stdin), FastIO(sys.stdout)
else:
    sys.stdin, sys.stdout = IOWrapper(sys.stdin), IOWrapper(sys.stdout)

input = lambda: sys.stdin.readline().rstrip("\r\n")

# endregion


if __name__ == "__main__":
    #read()
    main()
    #dmain()

# Comment Read()