Z = int(input())

for _ in range(Z):
    N,K = list(map(int,input().split()))

    Bd = [input() for _ in range(N)]

    cX = "".join(Bd).count('x')
    cO = "".join(Bd).count('o')

    poses = {
        'x': [],
        'o': [],
        '.': []
    }


    winseqs = []
    for i in range(N):
        for j in range(N):
            poses[Bd[i][j]].append((i,j))
            if Bd[i][j] not in 'xo': continue
            for di,dj in ((1,0),(0,1),(1,-1),(1,1)):
                ichk = i
                jchk = j

                seq = set([(i,j)])
                for _ in range(K-1):
                    ichk += di
                    jchk += dj
                    if ichk < 0 or ichk >= N or jchk < 0 or jchk >= N:
                        break
                    if Bd[ichk][jchk] != Bd[i][j]:
                        break
                    seq.add((ichk, jchk))
                else:
                    winseqs.append((Bd[i][j], seq))
    
    if len(set(x for x,y in winseqs)) > 1:
        print('NIE')
        continue
    

    if len(winseqs) == 0:
        if cX + cO != N*N:
            print('NIE')
            continue
        if cX > cO:
            mode = 'x'
            first = poses['x'][0]
        else:
            mode = 'o'
            first = poses['o'][0]
    else:
        mode = list(set(x for x,y in winseqs))[0]
        breakers = set.intersection(*(y for x,y in winseqs))
        #print(winseqs)
        if not len(breakers):
            print('NIE')
            continue
        first = list(breakers)[0]
    

    omode = 'o' if mode == 'x' else 'x'
    if len(poses[mode]) < len(poses[omode]):
        print('NIE')
        continue
    if len(poses[mode]) > len(poses[omode])+1:
        print('NIE')
        continue

    out = []
    out.append(first)
    poses[mode].remove(first)
    while len(poses[mode])+len(poses[omode]):
        omode,mode=mode,omode
        out.append(poses[mode].pop())
    

    print('TAK')
    for x,y in out[::-1]:
        print(x+1,y+1)