Problem:


博士の光電管

博士 :   ピーター君。やったよ。

ピーター : またですか。今度はどんなくだらない発明ですか。

博士 :   あの幻の素粒子アキシオンの検出器を発明したんじゃ。

ピーター : アキシオンといえば、欧州合同原子核研究機構 (CERN) をはじめとする研究者たちが血眼で追っかけているという、あれですよね。本当ですかぁ?

博士 :   本当だってばよ。細かい説明は省くが、非常に強力な磁場を内蔵する特殊な光電管が光ることによって、通過するアキシオンを検出する。

ピーター : 他 に先んじて検出すれば、小柴先生のニュートリノ検出に匹敵するノーベル賞級の研究ですよ。
       これで役立たずの研究ばかりしている「ダメ研究室」などという汚名も返上できますね。

博士 :   そうだとも。小柴先生の「スーパーカミオカンデ」にあやかって、この装置を、(悪口言ったら)「タダジャオカンデ」と命名した。

ピーター : ちょっと苦しいって言うか、卑屈って言うか・・・。

博士 :   それはいいとして、この装置ちょっとした癖があるんじゃよ。
       アキシオン粒子がある光電管を通過すると、感度の関係でその光電管と隣接する上下左右の光電管が反応する。

図1図2

	

	

	
		

		
★●●●●

		
●●●★●

		
●●●●●

		
●●●●●

		
●●★●●

		

	
	
		

		
 

		
 

		
   &rarr;   

		
 

		
 

		

	
	
		

		
○○●○●

		
○●○○○

		
●●●○●

		
●●○●●

		
●○○○●

		

	
	

	

	

	

	
		

		
●○○●○

		
○★★○☆

		
●●○●●

		
○●●○●

		
●○○○●

		

	

	
		

		
 

		
 

		
   &rarr;   

		
 

		
 

		

	
	
		

		
●●●●●

		
●●●○●

		
●○●●○

		
○●●○●

		
●○○○●

		

	
	

	

ピーター : つまり、図1の左側★印の光電管を粒子が通過した場合、
       右側のように点灯すると言うことですね。
       (図は 5 &times; 5 の例。黒は消灯・白は点灯状態。以下同じ。)

博士 :   また反応とは、光電管の状態が反転する。つまり消えている光電管は光り、
       光っている光電管は消える。

ピーター : つまり、図2の左側の★印や☆印を粒子が通過した場合、右側のような状態になると。

博士 :   これを正方形になんと 100 個 (10 &times; 10) 配置して待ち受けるわけだ。

ピーター : これだけの大発明、ノーベル賞選考委員会も「ホッチャオカンデ」なんて。

博士 :   おぉピーター君、君も我が研究室の作風になじんできたようだね。いい感じだ。
       では、早速実験を始めよう。まずはこの装置、現在ランダムに光電管が点いているので、
       実験が始められるように全部消えた状態にリセットしてくれたまえ。
       なあに、どの光電管にアキシオン粒子を当てれば全部消えるか考えればいいだけじゃよ。
       簡単だろう?

ピーター : 考えるのはいいんですが博士。当てるためには、
       幻のアキシオン粒子を発生させて打ち込める装置がないといけないんじゃないですか。

博士 :   ・・・。

博士とピータ      (同時に)コリャアカンデ!だははは・・・・。
ー:

ということで、本日も和気あいあいと暮れてゆく博士の研究室ですが、例によって話は全く進みそうもありません。しょうがないので、ピーター君に代わってプログラムを作成してあげてください。
プログラムは次のようなものになります。

 A. 装置の光電管の状態を 10 &times; 10 の配列として入力します。0は消灯中、1は点灯中を表します。
    0 と 1 以外のデータは含まれません。

 B. 入力した装置の状態をすべて消灯にするために、アキシオン粒子を通過させる位置を計算して出力します。これは、入力と同じ 10 &times; 10の配列で光電管の位置を表します。「0 は通過させない」、「1 は通過させる」ことを表します。すべてを消灯する方法は必ず 1 通りだけ存在するものとします。

Input

複数のデータセットが与えられます。１行目にデータセットの数 n (n ≤ 20) が与えられます。各データセットは以下の形式で与えられます。

a1,1 a1,2 ... a1,10
a2,1 a2,2 ... a2,10
:
a10,1 a10,2 ... a10,10

ai,j は装置の i 行目 j 列目の光電管の状態を示す整数 (0 または 1) を表します。

Output

各データセットについて、粒子を通過させる位置を以下の形式で出力してください。

b1,1 b1,2 ... b1,10
b2,1 b2,2 ... b2,10
:
b10,1 b10,2 ... b10,10

bi,j は装置の i 行目 j 列目の光電管に粒子を通過させるかいなかを示す整数 (0 または 1) を表します。

Sample Input

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

Output for the Sample Input

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



Sample input:
1
0 1 0 0 0 0 0 0 0 0
1 1 1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
0 0 0 0 1 1 0 0 0 0
0 0 0 1 0 0 1 0 0 0
0 0 0 0 1 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1 1 1
0 0 0 0 0 0 0 0 1 0


Sample output:
0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 1 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0 0 0

Submission status: Wrong Answer

Solution code:
def f(x,y):
  B[y][x]^=1
  if y>0: B[y-1][x]^=1
  if y<9: B[y+1][x]^=1
  if x>0: B[y][x-1]^=1
  if x<9: B[y][x+1]^=1
  return
R=range(10)
n=input()
for _ in [0]*n:
  A=[map(int,raw_input().split()) for i in R]
  for p in range(1024):
    B=A[:]
    x=[]
    for j in R:
      if j==0: a=map(int,list(format(p,'010b')))
      else: a=B[j-1]
      for i in R:
        if a[i]==1:
          x.append((i,j))
          f(i,j)
    if sum(B[9])==0: break
  for i,j in x: A[j][i]=1
  for e in A: print " ".join(map(str,e))

Submission status: Accepted

Solution code:
def f(x,y):
  global B
  B[y][x]^=1
  if y>0: B[y-1][x]^=1
  if y<9: B[y+1][x]^=1
  if x>0: B[y][x-1]^=1
  if x<9: B[y][x+1]^=1
  return
R=range(10)
n=input()
for _ in [0]*n:
  A=[map(int,raw_input().split()) for _ in R]
  for p in range(1024):
    B=[A[i][:] for i in R]
    x=[]
    for j in R:
      if j==0: a=map(int,list(format(p,'010b')))
      else: a=B[j-1]
      for i in R:
        if a[i]==1:
          x.append((i,j))
          f(i,j)
    if sum(B[9])==0: break
  for i,j in x: B[j][i]=1
  for e in B: print " ".join(map(str,e))

Submission status: Accepted

Solution code:
def f(x,y):
  global B
  B[y][x]^=1
  if y>0: B[y-1][x]^=1
  if y<9: B[y+1][x]^=1
  if x>0: B[y][x-1]^=1
  if x<9: B[y][x+1]^=1
  return (x,y)
R=range(10)
n=input()
for _ in [0]*n:
  A=[map(int,raw_input().split()) for _ in R]
  for p in range(1024):
    B=[e[:] for e in A]
    x=[]
    for j in R:
      a=[map(int,list(format(p,'010b'))),B[j-1]][j!=0]
      for i in R:
        if a[i]==1: x.append(f(i,j))
    if sum(B[9])==0: break
  for i,j in x: B[j][i]=1
  for e in B: print " ".join(map(str,e))

Submission status: Accepted

Solution code:
def f(x,y):
  global B
  B[y][x]^=1
  if y>0: B[y-1][x]^=1
  if y<9: B[y+1][x]^=1
  if x>0: B[y][x-1]^=1
  if x<9: B[y][x+1]^=1
  return (x,y)
R=range(10)
n=input()
for _ in [0]*n:
  A=[map(int,raw_input().split()) for _ in R]
  for p in range(1024):
    B=[e[:] for e in A]
    x=[]
    for j in R:
      if j==0: a=map(int,list(format(p,'010b')))
      else: a=B[j-1]
      for i in R:
        if a[i]==1: x.append(f(i,j))
    if sum(B[9])==0: break
  for i,j in x: B[j][i]=1
  for e in B: print " ".join(map(str,e))

Submission status: Accepted

Solution code:
import itertools
def change(m, x, y):
    m[y][x] ^= 1
    if y>0: m[y-1][x] ^= 1
    if x>0: m[y][x-1] ^= 1
    if x<9: m[y][x+1] ^= 1
    if y<9: m[y+1][x] ^= 1
ans = [[0] * 10 for i in xrange(10)]
for t in xrange(input()):
    m = [map(int, raw_input().split()) for i in xrange(10)]
    for p in itertools.product((0,1), repeat=10):
        mm = [m[i][:] for i in xrange(10)]
        for i in xrange(10):
            if p[i]:
                change(mm, i, 0)
            ans[0][i] = p[i]
        for i in xrange(9):
            for j in xrange(10):
                if mm[i][j]:
                    change(mm, j, i+1)
                    ans[i+1][j] = 1
                else:
                    ans[i+1][j] = 0
        if not any(mm[9]):
            for i in xrange(10):
                print " ".join(map(str, ans[i]))
            break

Submission status: Accepted

Solution code:
def attack(table, i, j):
    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]

def printans(ans):
    for i in range(1, 11):
        for j in range(1, 11):
            print(ans[i][j], end="")
            if j < 10:
                print(" ", end="")
        print("")

def solve(table, i, j, ans):

    #print(i,j)
    if i == 11:
        flag = True
        for k in range(1,11):
            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)


def check(table, i, ans):

    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)

N = int(input())

for l in range(N):
    table= [[0 for i in range(12)] for j in range(12)]
    ans= [[0 for i in range(12)] for j in range(12)]
    for i in range(1, 11):
        nums = [int(k) for k in input().split()]
        for j in range(1, 11):
            table[i][j] = nums[j-1]

    check(table, 1, ans)

Submission status: Accepted

Solution code:
from copy import deepcopy

def show(o):
    for i in range(10):
        print ' '.join(map(str, o[i]))

def e(o, i, j):
    px = [-1, 0, 0, 1, 0]
    py = [0, -1, 0, 0, 1]
    for k in range(5):
        if i + py[k] >= 0 and j + px[k] >= 0 and i + py[k] < 10 and j + px[k] < 10:
            o[i + py[k]][j + px[k]] = 1 - o[i + py[k]][j + px[k]]

def f(o, n):
    r = [[0] * 10 for i in range(10)]
    for i in range(10):
        if n >> i & 1 == 1:
            e(o, 0, i)
            r[0][i] = 1
    for i in range(1, 10):
        for j in range(10):
            if o[i - 1][j] == 1:
                e(o, i, j)
                r[i][j] = 1
    if 1 not in o[9]:
        show(r)
        return True
    return False

n = int(raw_input())
for _ in range(n):
    o = [map(int, raw_input().split()) for i in range(10)]
    for i in range(2 ** 10):
        if f(deepcopy(o), i):
            break

Submission status: Accepted

Solution code:
L = [(512,768,512)] + [(256>>i,896>>i,256>>i) for i in xrange(9)]

def put(f,x,y):
    for yy,b in zip(xrange(y-1,y+2),L[x]):
        f[yy] = f[yy] ^ b

def solve(F):
    for i in xrange(1<<10):
        f = [0] + F + [0]
        ret = []
        for x in (x for x in xrange(10) if 1<<x & i != 0):
            put(f,x,1)
            ret.append((x,0))
        f.pop(0)
        for y in xrange(9):
            for x in (x for x in xrange(10) if 1<<(9-x) & f[y] != 0):
                put(f,x,y+1)
                ret.append((x,y+1))

        if f[9] == 0:
            return ret
    
for _ in xrange(input()):
    F = [int(raw_input().replace(" ",""),2) for _ in xrange(10)]
    pos = solve(F)    
    print "\n".join([" ".join("1" if (x,y) in pos else "0" for x in xrange(10))\
                         for y in xrange(10)])

Submission status: Accepted

Solution code:
L = [(512,768,512)] + [(256>>i,896>>i,256>>i) for i in xrange(9)]

def put(f,x,y):
    for yy,b in zip(xrange(y-1,y+2),L[x]):
        f[yy] = f[yy] ^ b

def solve(F):
    for i in xrange(1<<10):
        f = [0] + F + [0]
        ret = []
        for x in (x for x in xrange(10) if (1<<x & i) != 0):
            put(f,x,1)
            ret.append((x,0))
        f.pop(0)
        for y in xrange(9):
            for x in (x for x in xrange(10) if (1<<(9-x) & f[y]) != 0):
                put(f,x,y+1)
                ret.append((x,y+1))
        if f[9] == 0:
            return ret
    
for _ in xrange(input()):
    F = [int(raw_input().replace(" ",""),2) for _ in xrange(10)]
    P = solve(F)    
    print "\n".join([" ".join("1" if (x,y) in P else "0" for x in xrange(10))\
                         for y in xrange(10)])

Submission status: Accepted

Solution code:
from itertools import product

def upd(line, fix):
  for i, v in enumerate(fix):
   if v:
     if i == 0:
       line[i] = not line[i]
       line[i + 1] = not line[i + 1]
     elif i == 9:
       line[i - 1] = not line[i - 1]
       line[i] = not line[i]
     else:
       line[i - 1] = not line[i - 1]
       line[i] = not line[i]
       line[i + 1] = not line[i + 1]
  return list(map(int, line))

def upd2(line, fix):
  return [int(i != j) for i, j in zip(line, fix)]

n = int(input())
for _ in range(n):
  base_mp = [list(map(int, input().split())) for _ in range(10)]
  for first_row in product([0, 1], repeat=10):
    mp = [list(first_row)] + [[i for i in base_mp[j]] for j in range(10)]
    save = []
    for i in range(10):
      save.append(mp[i])
      mp[i + 1] = upd(mp[i + 1], mp[i])
      if i == 9:continue
      mp[i + 2] = upd2(mp[i + 2], mp[i])
    
    if not any(mp[10]):
      for line in save:
        print(*line)
      break