CF 1440 C2. Binary Table (Hard Version)

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 x1,y1,x2,y2,x3,y3 (1≤x1,x2,x3≤n,1≤y1,y2,y3≤m) describing the next operation. This operation will be made with three cells (x1,y1), (x2,y2), (x3,y3). 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

代码

#include <bits/stdc++.h>
typedef long long ll;

const ll mod =1e9 + 7;

using namespace std;
namespace fastIO {
    inline void input(int& res) {
        char c = getchar();res = 0;int f = 1;
        while (!isdigit(c)) { f ^= c == '-'; c = getchar(); }
        while (isdigit(c)) { res = (res << 3) + (res << 1) + (c ^ 48);c = getchar(); }
        res = f ? res : -res;
    }
    inline ll qpow(ll a, ll b) {
        ll ans = 1, base = a;
        while (b) {
            if (b & 1) ans = (ans * base % mod +mod )%mod;
            base = (base * base % mod + mod)%mod;
            b >>= 1;
        }
        return ans;
    }
}
using namespace fastIO;

const int N = 1e2+5;
int Case,n,m;
char x;
int a[105][105];

struct node{
	int x1;
	int y1;
	int x2;
	int y2;
	int x3;
	int y3;
};
vector<node> zr;

void cal(int x1,int y1,int x2,int y2,int x3,int y3){
	a[x1][y1]^=1,a[x2][y2]^=1,a[x3][y3]^=1;
	zr.push_back({x1,y1,x2,y2,x3,y3});
}

void solve(){
	input(n),input(m);
	zr.clear();
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			scanf("%c",&x);
			a[i][j]=x-'0';
		}
		getchar();
	}
	for(int i=2;i<n;i++){
		for(int j=2;j<=m;j++){
			if(a[i-1][j-1]&&a[i-1][j])
                cal(i-1,j-1,i-1,j,i,j-1);
            else if(a[i-1][j-1])
                cal(i-1,j-1,i,j-1,i,j);
            else if(a[i-1][j])
                cal(i-1,j,i,j-1,i,j);
		}
	}
	for (int j=2;j<=m;j++){
        if(a[n-1][j-1]&&a[n][j-1])
            cal(n-1,j-1,n,j-1,n,j);
        else if(a[n-1][j-1])
            cal(n-1,j-1,n-1,j,n,j);
        else if(a[n][j-1])
            cal(n,j-1,n-1,j,n,j);
    }
    if(a[n][m]&&a[n-1][m]){
        cal(n,m,n,m-1,n-1,m-1);
        cal(n,m-1,n-1,m-1,n-1,m);
    }
    else if(a[n][m]){
        cal(n,m,n,m-1,n-1,m);
        cal(n-1,m-1,n,m,n-1,m);
        cal(n,m,n-1,m-1,n,m-1);
    }
    else if(a[n-1][m]){
        cal(n-1,m,n-1,m-1,n,m);
        cal(n,m,n-1,m,n,m-1);
        cal(n-1,m-1,n,m-1,n-1,m);
    }
	printf("%d\n",zr.size());
	for(auto t:zr){
		printf("%d %d %d %d %d %d\n",t.x1,t.y1,t.x2,t.y2,t.x3,t.y3);
	}
}

int main(){
	Case=1;
	input(Case);
	while(Case--){
		solve(); 
	}
	return 0;
}
/*
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
*/