Go: package main

import (
	"bufio"
	"fmt"
	"os"
	"strconv"
	"strings"
)

var scanner = bufio.NewScanner(os.Stdin)
var (
	table, ans = [12][12]int{}, [12][12]int{}
)

func main() {
	scanner.Scan()
	n, _ := strconv.Atoi(scanner.Text())
	for i := 0; i < n; i++ {
		for j := 1; j < 11; j++ {
			scanner.Scan()
			input := strings.Split(scanner.Text(), " ")
			tmp := [10]int{}
			for k := range input {
				tmp[k], _ = strconv.Atoi(input[k])
			}
			for l := 1; l < 11; l++ {
				table[j][l] = tmp[l-1]
			}
		}
		check(table, 1, ans)
	}
}

func check(table [12][12]int, i int, ans [12][12]int) {
	if i == 11 {
		solve(table, 2, 1, ans)
		return
	}
	ans[1][i] = 0
	check(table, i+1, ans)
	ans[1][i] = 1
	attack(table, 1, i)
	check(table, i+1, ans)
	attack(table, 1, i)
}

func solve(table [12][12]int, i, j int, ans [12][12]int) {
	if i == 11 {
		flag := true
		for k := 1; k < 11; k++ {
			if table[10][k] == 1 {
				flag = false
				break
			}
		}
		if flag {
			printans(ans)
		}
		return
	}
	if table[i-1][j] == 1 {
		ans[i][j] = 1
		attack(table, i, j)
		if j == 10 {
			solve(table, i+1, 1, ans)
		} else {
			solve(table, i, j+1, ans)
		}
		attack(table, i, j)
		ans[i][j] = 0
	} else {
		ans[i][j] = 0
		if j == 10 {
			solve(table, i+1, 1, ans)
		} else {
			solve(table, i, j+1, ans)
		}
	}
}

func printans(ans [12][12]int) {
	for i := 1; i < 11; i++ {
		for j := 1; j < 11; j++ {
			fmt.Print(ans[i][j])
			if j < 10 {
				fmt.Print(" ")
			}
		}
		fmt.Println()
	}
}

func attack(table [12][12]int, i, j int) {
	table[i][j] = 1 - table[i][j]
	table[i-1][j] = 1 - table[i-1][j]
	table[i+1][j] = 1 - table[i+1][j]
	table[i][j-1] = 1 - table[i][j-1]
	table[i][j+1] = 1 - table[i][j+1]
}

C++: #include<iostream>
#include<vector>
#include<queue>
#include<algorithm>
#include<cstring>
#include<cstdlib>

#define MAX 12
#define INV (-1)

using namespace std;

inline bool inRange(int t, int min, int max){
	return (min <= t && t < max);
}

class CMap{
public:
	int m[MAX][MAX];

	CMap()
	{
		for(int i = 0; i < MAX; ++i)
			for(int j = 0; j < MAX; ++j)
				m[i][j] = 0;
	}

	void PassParticle(int i, int j)
	{
		const int di[] = {0,-1,0,1,0};
		const int dj[] = {0,0,1,0,-1};

		for(int k = 0; k < sizeof(di)/sizeof(*di); ++k)
			if( inRange( i+di[k], 1, MAX - 1 ) && inRange( j + dj[k], 1, MAX - 1 ) )
				m[i+di[k]][j+dj[k]] = (m[i+di[k]][j+dj[k]] + 1) % 2;
	}

	bool NoLight() const
	{
		for(int i = 1; i < MAX - 1; ++i)
			for(int j = 1; j < MAX - 1; ++j)
				if( m[i][j] == 1 ) return false;
		return true;
	}

	void Print() const
	{
		for(int i = 1; i < MAX - 1; ++i){
			cout << m[i][1];
			for(int j = 2; j < MAX - 1; ++j)
				cout << ' ' << m[i][j];
			cout << '\n';
		}
	}
};
void DetermineShootingPoints(CMap &stm, CMap &ans){
	/* corner */
	{
		if( stm.m[1][1] == 1 ){
			if( stm.m[1][2] == 1 && stm.m[2][1] == 1 ){
				stm.PassParticle( 1, 1 );
				ans.m[1][1] = 1;
			}
		}
		if( stm.m[1][MAX-2] == 1 ){
			if( stm.m[1][MAX-3] == 1 && stm.m[2][MAX-2] == 1 ){
				stm.PassParticle( 1, MAX-2 );
				ans.m[1][MAX-2] = 1;
			}
		}
		if( stm.m[MAX-2][1] == 1 ){
			if( stm.m[MAX-3][1] == 1 && stm.m[MAX-2][2] == 1 ){
				stm.PassParticle( MAX-2, 1 );
				ans.m[MAX-2][1] = 1;
			}
		}
		if( stm.m[MAX-2][MAX-2] == 1 ){
			if( stm.m[MAX-2][MAX-3] == 1 && stm.m[MAX-3][MAX-2] == 1 ){
				stm.PassParticle( MAX-2, MAX-2 );
				ans.m[MAX-2][MAX-2] = 1;
			}
		}
	}

	while( !stm.NoLight() ){
		vector< pair<int, int > > vt;
		bool bSCont = true;

		/* deko */
		{
			const int ci[] = {1,1,1,1,4,5,6,7,4,5,6,7,MAX-2,MAX-2,MAX-2,MAX-2};
			const int cj[] = {4,5,6,7,1,1,1,1,MAX-2,MAX-2,MAX-2,MAX-2,4,5,6,7};

			for(int i = 0; i < sizeof(ci)/sizeof(*ci); ++i ){
				int cnt1 = 0;
				const int di[] = {0,-1,0,1,0};
				const int dj[] = {0,0,1,0,-1};

				for(int j = 0; j < sizeof(di)/sizeof(*di); ++j){
					if( stm.m[ci[i]+di[j]][cj[i]+dj[j]] == 1 &&
						inRange(ci[i]+di[j],1,MAX-1) && inRange(cj[i]+dj[j],1,MAX-1) )
						cnt1++;
				}
				if( cnt1 == 4 ){
					int ccnt1 = 0;
					const int ddi[] = {-1,-1,1,1};
					const int ddj[] = {1,1,-1,-1};

					for( int j = 0; j < sizeof(ddi)/sizeof(*ddi); ++j){
						if( stm.m[ci[i]+ddi[j]][cj[i]+ddj[j]] == 1 &&
							inRange(ci[i]+ddi[j],1,MAX-1) && inRange(ci[i]+ddi[j],1,MAX-1) )
							++ccnt1;
					}
					if( ccnt1 <= 1 ){
						stm.PassParticle( ci[i], cj[i] );
						ans.m[ci[i]][cj[i]] = 1;
					}
				}
			}
		}

		/* other */
		for(int i = 1; i < MAX - 1 && bSCont; ++i){
			for(int j = 1; j < MAX - 1 && bSCont; ++j){
				if( stm.m[i][j] == 1 ){
					
					const int di[][2] = { {0,-1}, {-1,0}, {0,1}, {1,0} };
					const int dj[][2] = { {-1,0}, {0,1}, {1,0}, {0,-1} };

					for(int k = 0; k < sizeof(di)/sizeof(*di); ++k){
						bool bCont = true;
						int ti1, tj1, ti2, tj2;

						for(int l = 0; l < sizeof(*di)/sizeof(**di) && bCont; ++l){
							if( stm.m[i+di[k][l]][j+dj[k][l]] != 0 ||
								!inRange(i+di[k][l],1,MAX-1) || !inRange(j+dj[k][l],1,MAX-1) ){
									ti1 = tj1 = ti2 = tj2 = INV;
									bCont = false;
							}else{
								switch( k ){
									case 0: ti1 = i + 1; tj1 = j; ti2 = i; tj2 = j + 1; break;
									case 1: ti1 = i; tj1 = j - 1; ti2 = i + 1; tj2 = j; break;
									case 2: ti1 = i - 1; tj1 = j; ti2 = i; tj2 = j - 1; break;
									case 3: ti1 = i; tj1 = j + 1; ti2 = i - 1; tj2 = j; break;
								}
							}
						}
						if( bCont ){
							if( inRange(ti1,1,MAX-1) && inRange(tj1,1,MAX-1) )
								vt.push_back( make_pair( ti1, tj1 ) );
							if( inRange(ti2,1,MAX-1) && inRange(tj2,1,MAX-1) )
								vt.push_back( make_pair( ti2, tj2 ) );
						}
					}

					
				}
			}
		}
		if( vt.size() > 0 ){
			int ti = INV, tj = INV;
			int max1 = 0;
			double maxEfficient = 0.499;
			
			for(unsigned int t = 0; t < vt.size(); ++t){
				int cnt1 = 0;
				int invCnt = 0;
				const int di[] = {0,-1,0,1,0};
				const int dj[] = {0,0,1,0,-1};

				for(int k = 0; k < sizeof(di)/sizeof(*di); ++k){
					if( !inRange(vt[t].first+di[k], 1, MAX-1) || !inRange(vt[t].second+dj[k],1,MAX-1) )
						invCnt++;
					if( stm.m[vt[t].first+di[k]][vt[t].second+dj[k]] == 1 )
						cnt1++;
				}

				if( cnt1 > max1 ){
					max1 = cnt1;
					if( (double)cnt1 / ( 5 - invCnt ) > maxEfficient ){
						maxEfficient = (double)cnt1 / (5 - invCnt);
						ti = vt[t].first;
						tj = vt[t].second;	
					}
				}
			}
			if( ti != INV && tj != INV ){
				//cout << ti << ' ' << tj << '\n';
				stm.PassParticle( ti, tj );
				ans.m[ti][tj] = 1;
				bSCont = false;
			}
		}
	}
}
int main(void){
	int n;

	cin >> n;
	for( int D = 0; D < n; ++D )
	{
		CMap M;
		CMap Answer;

		for(int i = 1; i < MAX - 1; ++i)
			for(int j = 1; j < MAX - 1 ; ++j)
				cin >> M.m[i][j];

		DetermineShootingPoints( M, Answer );
		Answer.Print();
	}
	return 0;
}

C++: #include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

#define FOR(i,a,b) for(int i=(a);i<(b);i++)
#define REP(i,n) FOR(i,0,n)

void turn(int x, int y, int data[10][10]) {
	int dx[5] = {0,1,0,-1,0};
	int dy[5] = {0,0,1,0,-1};
	
	REP(i, 5) {
		int nx = x + dx[i], ny = y + dy[i];
		if (!(nx >= 0 && nx < 10 && ny >= 0 && ny < 10)) continue;
		data[ny][nx]++;
	}
}

bool check(int data[10][10]) {
	bool res = true;
	
	REP(i, 10) if (data[9][i] & 1) res = false;
	
	return res;
}

int main() {
	int n;
	cin >> n;
	while (n--) {
		int in[10][10];
		REP(i, 10) REP(j, 10) scanf("%d", &in[i][j]);
		
		bool ans[10][10];
		int data[10][10];
		
		REP(i, 1<<10) {
			memset(ans, 0, sizeof(ans));
			memcpy(data, in, sizeof(in));
			
			REP(j, 10) if ((i >> j) & 1) {
				ans[0][j] = true;
				turn(j, 0, data);
			}
			FOR(j, 1, 10) REP(k, 10) {	
				if (data[j - 1][k] & 1) {
					ans[j][k] = true;
					turn(k, j, data);
				}
			}
			
			if (check(data)) break;
		}
		
		REP(j, 10) REP(k, 10)
			printf("%d%c", ans[j][k] ? 1 : 0, k == 9 ? '\n' : ' ');
	}
	return 0;
}