yur.py

#!/usr/bin/env python
# -*- coding: utf-8 -*-

'''Module to solve Flow Free puzzles by reducing to SAT and invoking
pycosat's solver. The reduction to SAT does not automatically prevent
cycles from appearing which are disconnected from the pre-colored
cells; however, they are detected, and the solver will be run until
all cycles are eliminated.

'''

import os
import sys
import operator
import itertools
from datetime import datetime
from argparse import ArgumentParser
from collections import defaultdict
import pycosat
import functools


LEFT = 1
RIGHT = 2
TOP = 4
BOTTOM = 8

DELTAS = [(LEFT, 0, -1),
          (RIGHT, 0, 1),
          (TOP, -1, 0),
          (BOTTOM, 1, 0)]

LR = LEFT | RIGHT
TB = TOP | BOTTOM
TL = TOP | LEFT
TR = TOP | RIGHT
BL = BOTTOM | LEFT
BR = BOTTOM | RIGHT

DIR_TYPES = [LR, TB, TL, TR, BL, BR]

# DOES YOUR INDIVIDUAL CELL HAVE A CONNECTION TO THE LEFT & RIGHT OR DOES IT HAVE TO TOP AND RIGHT

DIR_FLIP = {
    LEFT: RIGHT,
    RIGHT: LEFT,
    TOP: BOTTOM,
    BOTTOM: TOP
    }

ANSI_LOOKUP = dict(R=101, B=104, Y=103, G=42,
                   O=43, C=106, M=105, m=41,
                   P=45, A=100, W=107, g=102,
                   T=47, b=44, c=46, p=35)

ANSI_RESET = '\033[0m'

ANSI_CELL_FORMAT = '\033[30;{}m'

DIR_LOOKUP = {
    LR: '─',
    TB: '│',
    TL: '┘',
    TR: '└',
    BL: '┐',
    BR: '┌'
    }

RESULT_STRINGS = dict(s='successful',
                      f='failed',
                      u='unsolvable')

######################################################################

def all_pairs(collection):
    '''Return all combinations of two items from a collection, useful for
making a large number of SAT variables mutually exclusive.

    '''

    return itertools.combinations(collection, 2)

######################################################################

def no_two(satvars):
    '''Given a collection of SAT variables, generates clauses specifying
that no two of them can be true at the same time.

    '''

    return ((-a, -b) for (a, b) in all_pairs(satvars))

######################################################################

def explode(puzzle):
    '''Iterator helper function to allow looping over 2D arrays without
nested 'for' loops.

    '''

    for i, row in enumerate(puzzle):
        for j, char in enumerate(row):
            yield i, j, char

######################################################################

def valid_pos(size, i, j):
    '''Check whether a position on a square grid is valid.'''
    return i >= 0 and i < size and j >= 0 and j < size

######################################################################

def all_neighbors(i, j):
    '''Return all neighbors of a grid square at row i, column j.'''
    return ((dir_bit, i+delta_i, j+delta_j)
            for (dir_bit, delta_i, delta_j) in DELTAS)

######################################################################

def valid_neighbors(size, i, j):
    '''Return all actual on-grid neighbors of a grid square at row i,
column j.'''

    return ((dir_bit, ni, nj) for (dir_bit, ni, nj)
            in all_neighbors(i, j)
            if valid_pos(size, ni, nj))

######################################################################
def repair_colors(puzzle, colors):

    '''If the puzzle file used "color labels" (A,B,C...), instead
    of color mnemonics (R,G,B...), convert the labels to mnemonics.
    Note: If a puzzle is in mnemonic format, it will contain the
    letter 'R' since Red is always the first color. Absence of an 'R'
    is the test for a color-labeled puzzle.
    '''
    if 'R' in colors.keys():
        return puzzle, colors

    color_lookup = 'RBYGOCMmPAWgTbcp'
    new_puzzle = []

    try:

        for row in puzzle:

            new_row = []

            for char in row:
                if char.isalnum():
                    char = color_lookup[ord(char)-ord('A')]
                new_row.append(char)

            new_puzzle.append(''.join(new_row))

        new_colors = dict((color_lookup[ord(char)-ord('A')], index)
                          for (char, index) in colors.items())

    except IndexError:
        return puzzle, colors

    return new_puzzle, new_colors

######################################################################

def parse_puzzle(options, file_or_str, filename='input'):

    '''Convert the given string or file object into a square array of
strings. Also return a dictionary which maps input characters to color
indices.

    '''

    if not isinstance(file_or_str, str):
        file_or_str = file_or_str.read()

    puzzle = file_or_str.splitlines()

    # assume size based on length of first line
    size = len(puzzle[0])

    # make sure enough lines
    if len(puzzle) < size:
        print('{}:{} unexpected EOF'.format(filename, len(puzzle)+1))
        return None, None

    # truncate extraneous lines
    puzzle = puzzle[:size]

    # count colors and build lookup
    colors = dict()
    color_count = []

    for i, row in enumerate(puzzle):
        if len(row) != size:
            print('{}:{} row size mismatch'.format(filename, i+1))
            return None, None
        for j, char in enumerate(row):
            if char.isalnum(): # flow endpoint
                if char in colors:
                    color = colors[char]
                    if color_count[color]:
                        print('{}:{}:{} too many {} already'.format(
                            filename, i+1, j, char))
                        return None, None
                    color_count[color] = 1
                else:
                    color = len(colors)
                    colors[char] = color
                    color_count.append(0)

    # check parity
    for char, color in colors.items():
        if not color_count[color]:
            print('color {} has start but no end!'.format(char))
            return None, None

    # print info
    if not options.quiet:
        print('read {}x{} puzzle with {} colors from {}'.format(
            size, size, len(colors), filename))
        print

    puzzle, colors = repair_colors(puzzle, colors)
    return puzzle, colors

######################################################################

def make_color_clauses(puzzle, colors, color_var):

    '''Generate CNF clauses entailing the N*M color SAT variables, where N
is the number of cells and M is the number of colors. Each cell
encodes a single color in a one-hot fashion.

    '''

    clauses = []
    num_colors = len(colors)
    size = len(puzzle)

    # for each cell
    for i, j, char in explode(puzzle):

        if char.isalnum(): # flow endpoint

            endpoint_color = colors[char]

            # color in this cell is this one
            clauses.append([color_var(i, j, endpoint_color)])

            # color in this cell is not the other ones
            for other_color in range(num_colors):
                if other_color != endpoint_color:
                    clauses.append([-color_var(i, j, other_color)])

            # gather neighbors' variables for this color
            neighbor_vars = [color_var(ni, nj, endpoint_color) for
                             _, ni, nj in valid_neighbors(size, i, j)]

            # one neighbor has this color
            clauses.append(neighbor_vars)

            # no two neighbors have this color
            clauses.extend(no_two(neighbor_vars))

        else:

            # one of the colors in this cell is set
            clauses.append([color_var(i, j, color)
                            for color in range(num_colors)])

            # no two of the colors in this cell are set
            cell_color_vars = (color_var(i, j, color) for
                               color in range(num_colors))

            clauses.extend(no_two(cell_color_vars))

    return clauses

######################################################################

def make_dir_vars(puzzle, start_var):

    '''Creates the direction-type SAT variables for each cell.'''

    size = len(puzzle)
    dir_vars = dict()
    num_dir_vars = 0

    for i, j, char in explode(puzzle):

        if char.isalnum(): # flow endpoint, no dir needed
            continue

        # collect bits for neighbors (TOP BOTTOM LEFT RIGHT)
        neighbor_bits = (dir_bit for (dir_bit, ni, nj)
                         in valid_neighbors(size, i, j))

        # OR them all together
        cell_flags = functools.reduce(operator.or_, neighbor_bits, 0)

        # create a lookup for dir type vars in this cell
        dir_vars[i, j] = dict()

        for code in DIR_TYPES:
            # only add var if cell has correct flags (i.e. if cell has
            # TOP, BOTTOM, RIGHT, don't add LR).

            # ONLY ADDES VALID COMBOS OF DIRECTIONS, SO IF A CELL HAS VALID NEIGHBORS UP DOWN AND LEFT THEN IT WONT ANT BOTTOM RIGHT, OR LEFT RIGHT, OR TOP RIGHT TO THE POTENTIAL DIRECTIONS
            if cell_flags & code == code:
                num_dir_vars += 1
                dir_vars[i, j][code] = start_var + num_dir_vars

    return dir_vars, num_dir_vars

######################################################################

def make_dir_clauses(puzzle, colors, color_var, dir_vars):

    '''Generate clauses involving the color and direction-type SAT
variables. Each free cell must be exactly one direction, and
directions imply color matching with neighbors.

    '''

    dir_clauses = []
    num_colors = len(colors)
    size = len(puzzle)

    # for each cell
    for i, j, char in explode(puzzle):

        if char.isalnum(): # flow endpoint, no dir needed
            continue

        cell_dir_dict = dir_vars[(i, j)]
        cell_dir_vars = cell_dir_dict.values()

        # at least one direction is set in this cell
        dir_clauses.append(cell_dir_vars)

        # no two directions are set in this cell
        dir_clauses.extend(no_two(cell_dir_vars))

        # for each color
        for color in range(num_colors):

            # get color var for this cell
            color_1 = color_var(i, j, color)

            # for each neighbor
            for dir_bit, n_i, n_j in all_neighbors(i, j):

                # get color var for other cell
                color_2 = color_var(n_i, n_j, color)

                # for each direction variable in this scell
                for dir_type, dir_var in cell_dir_dict.items():

                    # if neighbor is hit by this direction type
                    if dir_type & dir_bit:
                        # this direction type implies the colors are equal
                        dir_clauses.append([-dir_var, -color_1, color_2]) # COLOR 2 AND COLOR 1 HAVE TO BE TRUE IS ESSENTIALLY WHAT THIS DOES BUT IT EQUATES THAT TO ONLY HAPPEN IF THIS IS THE CORRECT DIRECTION
                        dir_clauses.append([-dir_var, color_1, -color_2]) # STATMENT GETS IGNORED IF THIS IS NOT DIRECTION
                    elif valid_pos(size, n_i, n_j):
                        # neighbor is not along this direction type,
                        # so this direction type implies the colors are not equal
                        dir_clauses.append([-dir_var, -color_1, -color_2]) # IF WE ARE THIS DIRECTION, THEN WE CAN'T BE COLOR 1 OR COLOR 2

    return dir_clauses

######################################################################

def reduce_to_sat(options, puzzle, colors):

    '''Reduces the given puzzle to a SAT problem specified in CNF. Returns
a list of clauses where each clause is a list of single SAT variables,
possibly negated.

    '''

    size = len(puzzle)
    num_colors = len(colors)

    num_cells = size**2
    num_color_vars = num_colors * num_cells

    def color_var(i, j, color):
        '''Return the index of the SAT variable for the given color in row i,
 column j.

        '''
        return (i*size + j)*num_colors + color + 1

    start = datetime.now()

    color_clauses = make_color_clauses(puzzle,
                                       colors,
                                       color_var)

    dir_vars, num_dir_vars = make_dir_vars(puzzle, num_color_vars)

    dir_clauses = make_dir_clauses(puzzle, colors,
                                   color_var, dir_vars)

    num_vars = num_color_vars + num_dir_vars
    clauses = color_clauses + dir_clauses

    reduce_time = (datetime.now() - start).total_seconds()

    if not options.quiet:
        print('generated {:,} clauses over {:,} color variables'.format(
            len(color_clauses), num_color_vars, grouping=True))

        print( 'generated {:,} dir clauses over {:,} dir variables'.format(
            len(dir_clauses), num_dir_vars))

        print( 'total {:,} clauses over {:,} variables'.format(
            len(clauses), num_vars))

        print( 'reduced to SAT in {:.3f} seconds'.format(reduce_time))
        print()

    return color_var, dir_vars, num_vars, clauses, reduce_time

######################################################################

def decode_solution(puzzle, colors, color_var, dir_vars, sol):

    '''Takes the solution set from SAT and decodes it by undoing the
one-hot encoding in each cell for color and direction-type. Returns a
2D array of (color, direction-type) pairs.

    '''

    sol = set(sol)
    num_colors = len(colors)

    decoded = []

    for i, row in enumerate(puzzle):

        decoded_row = []

        for j, char in enumerate(row):

            # find which color variable for this cell is in the
            # solution set
            cell_color = -1

            for color in range(num_colors):
                if color_var(i, j, color) in sol:
                    assert cell_color == -1
                    cell_color = color

            assert cell_color != -1

            cell_dir_type = -1

            if not char.isalnum(): # not a flow endpoint

                # find which dir type variable for this cell is in the
                # solution set
                for dir_type, dir_var in dir_vars[i, j].items():
                    if dir_var in sol:
                        assert cell_dir_type == -1
                        cell_dir_type = dir_type

                assert cell_dir_type != -1

            decoded_row.append((cell_color, cell_dir_type))

        decoded.append(decoded_row)

    return decoded

######################################################################

def make_path(decoded, visited, cur_i, cur_j):

    '''Follow a path starting from an arbitrary row, column location on
the grid until a non-path cell is detected, or a cycle is
found. Returns a list of (row, column) pairs on the path, as well as a
boolean flag indicating if a cycle was detected.

    '''

    size = len(decoded)

    run = []
    is_cycle = False
    prev_i, prev_j = -1, -1

    while True:

        advanced = False

        # get current cell, set visited, add to path
        color, dir_type = decoded[cur_i][cur_j]
        visited[cur_i][cur_j] = 1
        run.append((cur_i, cur_j))

        # loop over valid neighbors
        for dir_bit, n_i, n_j in valid_neighbors(size, cur_i, cur_j):

            # do not consider prev pos
            if (n_i, n_j) == (prev_i, prev_j):
                continue

            # get neighbor color & dir type
            n_color, n_dir_type = decoded[n_i][n_j]

            # these are connected if one of the two dir type variables
            # includes the (possibly flipped) direction bit.
            if ((dir_type >= 0 and (dir_type & dir_bit)) or
                    (dir_type == -1 and n_dir_type >= 0 and
                     n_dir_type & DIR_FLIP[dir_bit])):

                # if connected, they better be the same color
                assert color == n_color

                # detect cycle
                if visited[n_i][n_j]:
                    is_cycle = True
                else:
                    prev_i, prev_j = cur_i, cur_j
                    cur_i, cur_j = n_i, n_j
                    advanced = True

                # either cycle detected or path advanced, so stop
                # looking at neighbors
                break

        # if path not advanced, quit
        if not advanced:
            break

    return run, is_cycle

######################################################################

def detect_cycles(decoded, dir_vars):

    '''Examine the decoded SAT solution to see if any cycles exist; if so,
return the CNF clauses that need to be added to the problem in order
to prevent them.

    '''

    size = len(decoded)
    colors_seen = set()
    visited = [[0]*size for _ in range(size)]

    # for each cell
    for i, j, (color, dir_type) in explode(decoded):

        # if flow endpoint for color we haven't dealt with yet
        if dir_type == -1 and color not in colors_seen:

            # add it to set of colors dealt with
            assert not visited[i][j]
            colors_seen.add(color)

            # mark the path as visited
            run, is_cycle = make_path(decoded, visited, i, j)
            assert not is_cycle

    # see if there are any unvisited cells, if so they have cycles
    extra_clauses = []

    for i, j in itertools.product(range(size), range(size)):

        if not visited[i][j]:

            # get the path
            run, is_cycle = make_path(decoded, visited, i, j)
            assert is_cycle

            # generate a clause negating the conjunction of all
            # direction types along the cycle path.
            clause = []

            for r_i, r_j in run:
                _, dir_type = decoded[r_i][r_j]
                dir_var = dir_vars[r_i, r_j][dir_type]
                clause.append(-dir_var)

            extra_clauses.append(clause)

    # return whatever clauses we had to generate
    return extra_clauses

######################################################################

def show_solution(options, colors, decoded):

    '''Print the puzzle solution to the terminal.'''

    # make an array to flip the key/value in the colors dict so we can
    # index characters numerically:

    color_chars = [None]*len(colors)

    do_color = options.display_color

    for char, color in colors.items():
        color_chars[color] = char
        do_color = do_color and char in ANSI_LOOKUP

    for decoded_row in decoded:
        for (color, dir_type) in decoded_row:

            assert color >= 0 and color < len(colors)

            color_char = color_chars[color]

            if dir_type == -1:
                if do_color:
                    display_char = 'O'
                else:
                    display_char = color_char
            else:
                display_char = DIR_LOOKUP[dir_type]

            if do_color:

                if color_char in ANSI_LOOKUP:
                    ansi_code = ANSI_CELL_FORMAT.format(
                        ANSI_LOOKUP[color_char])
                else:
                    ansi_code = ANSI_RESET

                sys.stdout.write(ansi_code)

            sys.stdout.write(display_char)

        if options.display_color:
            sys.stdout.write(ANSI_RESET)

        sys.stdout.write('\n')

######################################################################

def solve_sat(options, puzzle, colors, color_var, dir_vars, clauses):

    '''Solve the SAT now that it has been reduced to a list of clauses in
CNF.  This is an iterative process: first we try to solve a SAT, then
we detect cycles. If cycles are found, they are prevented from
recurring, and the next iteration begins. Returns the SAT solution
set, the decoded puzzle solution, and the number of cycle repairs
needed.

    '''

    start = datetime.now()

    decoded = None
    all_decoded = []
    repairs = 0

    while True:

        sol = pycosat.solve(clauses) # pylint: disable=E1101

        if not isinstance(sol, list):
            decoded = None
            all_decoded.append(decoded)
            break

        decoded = decode_solution(puzzle, colors, color_var, dir_vars, sol)
        all_decoded.append(decoded)

        extra_clauses = detect_cycles(decoded, dir_vars)

        if not extra_clauses:
            break

        clauses += extra_clauses
        repairs += 1

    solve_time = (datetime.now() - start).total_seconds()

    if not options.quiet:
        if options.display_cycles:
            for cycle_decoded in all_decoded[:-1]:
                print( 'intermediate solution with cycles:')
                print()
                show_solution(options, colors, cycle_decoded)
                print

        if decoded is None:
            print( 'solver returned {} after {:,} cycle '\
                'repairs and {:.3f} seconds'.format(
                    str(sol), repairs, solve_time))

        else:
            print( 'obtained solution after {:,} cycle repairs '\
                'and {:.3f} seconds:'.format(
                    repairs, solve_time))
            print
            show_solution(options, colors, decoded)
            print

    return sol, decoded, repairs, solve_time

######################################################################

def print_summary(options, stats):

    '''Print out stats for all solutions.'''

    max_width = max(len(f) for f in options.filenames)

    solution_types = stats.keys()

    all_stats = defaultdict(float)

    for result_char in solution_types:
        for k in stats[result_char]:
            all_stats[k] += stats[result_char][k]

    if all_stats['count'] > 1:

        if not options.quiet:

            print( '\n'+('*'*70)+'\n')

            for result_char in solution_types:

                print( '{:d} {:s} searches took:\n'\
                    '  {:,.3f} sec. to reduce '\
                    '(with {:,d} variables and {:,d} clauses)\n'\
                    '  {:,.3f} sec. to solve (with {:d} repairs)\n'\
                    '  {:,.3f} sec. total\n'.format(
                        stats[result_char]['count'], result_char,
                        stats[result_char]['reduce_time'],
                        stats[result_char]['num_vars'],
                        stats[result_char]['num_clauses'],
                        stats[result_char]['solve_time'],
                        stats[result_char]['repairs'],
                        stats[result_char]['total_time']))

            if len(solution_types) > 1:

                print( 'overall, {:d} searches took:\n'\
                    '  {:,.3f} sec. to reduce '\
                    '(with {:,d} variables and {:,d} clauses)\n'\
                    '  {:,.3f} sec. to solve (with {:d} repairs)\n'\
                    '  {:,.3f} sec. total\n'.format(
                        int(all_stats['count']),
                        all_stats['reduce_time'],
                        int(all_stats['num_vars']),
                        int(all_stats['num_clauses']),
                        all_stats['solve_time'],
                        int(all_stats['repairs']),
                        all_stats['total_time']))

        else:

            print

            for result_char in solution_types:

                print( '{:s}{:3d} total {:s} {:9,d} {:9,d} {:12,.3f} '\
                    '{:3d} {:12,.3f} {:12,.3f}'.format(
                        ' '*(max_width-9), stats[result_char]['count'],
                        result_char,
                        stats[result_char]['num_vars'],
                        stats[result_char]['num_clauses'],
                        stats[result_char]['reduce_time'],
                        stats[result_char]['repairs'],
                        stats[result_char]['solve_time'],
                        stats[result_char]['total_time']))

            if len(solution_types) > 1:

                print( '{:s}{:3d} overall {:9,d} {:9,d} {:12,.3f} '\
                    '{:3d} {:12,.3f} {:12,.3f}'.format(
                        ' '*(max_width-9), int(all_stats['count']),
                        int(all_stats['num_vars']), int(all_stats['num_clauses']),
                        all_stats['reduce_time'], int(all_stats['repairs']),
                        all_stats['solve_time'], all_stats['total_time']))

######################################################################

def pyflow_solver_main():

    '''Main loop if module run as script.'''

    color_capable = (sys.platform != 'win32' and os.isatty(1))

    parser = ArgumentParser(
        description='Solve Flow Free puzzles via reduction to SAT')

    parser.add_argument('filenames', metavar='PUZZLE', nargs='+',
                        help='puzzle file to load')

    parser.add_argument('-q', dest='quiet', default=False,
                        action='store_true',
                        help='quiet mode (reduce output)')

    parser.add_argument('-c', dest='display_cycles', default=False,
                        action='store_true',
                        help='display intermediate solutions with cycles')

    parser.add_argument('-C', dest='display_color', default=color_capable,
                        action='store_true',
                        help='always display color')

    options = parser.parse_args()

    max_width = max(len(f) for f in options.filenames)

    puzzle_count = 0

    stats = dict()

    for filename in options.filenames:

        if not options.quiet and puzzle_count:
            print( '\n'+('*'*70)+'\n')

        # open file
        try:
            with open(filename, 'r') as infile:
                puzzle, colors = parse_puzzle(options, infile, filename)
        except IOError:
            print( '{}: error opening file'.format(filename))
            continue

        if colors is None:
            continue

        puzzle_count += 1

        color_var, dir_vars, num_vars, clauses, reduce_time = \
            reduce_to_sat(options, puzzle, colors)

        sol, _, repairs, solve_time = solve_sat(options, puzzle, colors,
                                                color_var, dir_vars, clauses)

        total_time = reduce_time + solve_time

        if isinstance(sol, list):
            result_char = 's'
        elif str(sol) == 'UNSAT':
            result_char = 'u'
        else:
            result_char = 'f'

        cur_stats = dict(repairs=repairs,
                         reduce_time=reduce_time,
                         solve_time=solve_time,
                         total_time=total_time,
                         num_vars=num_vars,
                         num_clauses=len(clauses),
                         count=1)

        if result_char not in stats:
            stats[result_char] = cur_stats
        else:
            for key in cur_stats.keys():
                stats[result_char][key] += cur_stats[key]

        if not options.quiet:
            print( 'finished in total of {:.3f} seconds'.format(
                total_time))
        else:

            print( '{:>{}s} {} {:9,d} {:9,d} {:12,.3f} '\
                '{:3d} {:12,.3f} {:12,.3f}'.format(
                    filename, max_width, result_char,
                    num_vars, len(clauses), reduce_time,
                    repairs, solve_time, total_time))


    print_summary(options, stats)


######################################################################

if __name__ == '__main__':

    pyflow_solver_main()