Problem:
H - 植林

問題文

とある大学の魔法学部の生徒であるK君は広大な屋敷に住んでいる．
屋敷の裏には林が広がっていたが，所々木が生えていない箇所があった．
このままでは屋敷からの見た目が悪いと感じたので，K君は木が生えていない箇所に木を植えることに決めた．
魔法使いの卵であるK君は使い魔を召喚することが出来るので，植林の作業を使い魔にやらせようと考えた．

林は H × W 個のセルを持つ長方形領域とみなすことが出来る．
各セルは木が1本生えているか1本も生えていないかのどちらかである．
K君はこの H × W 個のセルのどれか 1 つを指定し，そこに使い魔を召喚する．
しかしK君はまだ未熟者なので，魔力を調節することが出来ず，使い魔を召喚するときは必ず 5 匹召喚してしまう．
さらに悪いことに，K君が召喚する使い魔はひねくれもので，訪れたセルに木が生えていない場合はそこに木を植えるが，訪れたセルに既に木が生えている場合はそのセルの木を消してしまう．
召喚された使い魔のうちの 4 匹は，指定されたセルを始点として東西南北に散らばり1直線上を進み，1つ1つセルを訪れていく．これらの使い魔は林の外に出ると消える．
残りの使い魔は指定されたセルのみを訪問し，その後直ちに消える．

より正確に言えば，K君がセル (i, j) に使い魔を召喚すると，i 行目か或いは
j 列目にある H+W-1 個のセルに対し，
木が生えていなければ木が植えられ，木が生えていれば木が消されるという操作が行われる．

召喚には多くの魔力を必要とするので，出来るだけ少ない召喚回数で林を木で覆いつくしたい．
K君はどのセルに使い魔を召喚すれば最小の召喚回数で林を木で覆いつくすことが出来るだろうか．

入力形式

入力は以下の形式で与えられる．

H W
a1,1 ... a1,W
...
aH,1 ... aH,W

ai,j が 1 ならセル (i,j) に木が生えていることを意味し,
0なら生えていないことを意味する．

出力形式

もし林を木で覆いつくすことが不可能なら，1行に Impossible を出力せよ．
そうでなければ，召喚回数を最小にするような召喚手順を以下の形式で H 行に出力せよ．

b1,1 ... b1,W
...
bH,1 ... bH,W

bi,j は 0 もしくは 1 でなければならず,
1 ならセル (i,j) に使い魔を召喚する事を意味し，0 なら召喚しない事を意味する．

2回同じ場所に使い魔を召喚しても意味が無いことに注意せよ．召喚回数を最小にするような召喚手順が複数ある場合はどれを出力しても良い．

制約

2 ≤ H，W ≤ 1000

H，W は偶数

ai,j ∈ {0, 1}

この問題の判定には，3 点分のテストケースのグループが設定されている． このグループに含まれるテストケースは上記の制約に加えて下記の制約も満たす．

H × W ≤ 20

入出力例

入力例 1

4 4
0 0 0 1
0 1 1 0
0 1 1 0
1 0 0 0

出力例 1

1 0 0 0
0 0 0 0
0 0 0 0
0 0 0 1

セル (1,1) とセル (4,4) に使い魔を召喚することで林を木で覆いつくすことが出来る．

Writer: 田村和範
Tester: 花田裕一朗

Solution in C++14 (judge verdict: Accepted)

#include <bits/stdc++.h>

#define _overload(_1,_2,_3,name,...) name
#define _rep(i,n) _range(i,0,n)
#define _range(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload(__VA_ARGS__,_range,_rep,)(__VA_ARGS__)

#define _rrep(i,n) _rrange(i,n,0)
#define _rrange(i,a,b) for(int i=int(a)-1;i>=int(b);--i)
#define rrep(...) _overload(__VA_ARGS__,_rrange,_rrep,)(__VA_ARGS__)

#define _all(arg) begin(arg),end(arg)
#define uniq(arg) sort(_all(arg)),(arg).erase(unique(_all(arg)),end(arg))
#define getidx(ary,key) lower_bound(_all(ary),key)-begin(ary)
#define clr(a,b) memset((a),(b),sizeof(a))
#define bit(n) (1LL<<(n))
#define popcount(n) (__builtin_popcountll(n))

using namespace std;

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

using ll=long long;
using R=long double;
const R EPS=1e-9L; // [-1000,1000]->EPS=1e-8 [-10000,10000]->EPS=1e-7
inline int sgn(const R& r){return(r > EPS)-(r < -EPS);}
inline R sq(R x){return sqrt(max(x,0.0L));}

const int dx[8]={1,0,-1,0,1,-1,-1,1};
const int dy[8]={0,1,0,-1,1,1,-1,-1};

// Problem Specific Parameter:

int h,w;
int a[1010][1010];
int fh[1010],fw[1010];

int main(void){
	cin >> h >> w;
	rep(i,h)rep(j,w) cin >> a[i][j];

	rep(i,h)rep(j,w){
		if(a[i][j]==0){
			fh[i]++;
			fw[j]++;
		}
	}
	
	rep(i,h){
		rep(j,w) cout << (j?" ":"") << (fh[i]+fw[j]+(a[i][j]==0))%2;
		cout << endl;
	}

	return 0;
}