Implementar un Solucionador de Sudoku de Fuerza Bruta

20

Implemente el solucionador de Sudoku más corto usando adivinanzas. Como he recibido algunas solicitudes, he agregado esto como una pregunta alternativa para aquellos que desean implementar un solucionador de sudoku de fuerza bruta.

Sudoku Puzzle:

 | 1 2 3 | 4 5 6 | 7 8 9
-+-----------------------
A|   3   |     1 |
B|     6 |       |   5
C| 5     |       | 9 8 3
-+-----------------------
D|   8   |     6 | 3   2
E|       |   5   |
F| 9   3 | 8     |   6
-+-----------------------
G| 7 1 4 |       |     9
H|   2   |       | 8
I|       | 4     |   3

Responder:

 | 1 2 3 | 4 5 6 | 7 8 9
-+-----------------------
A| 8 3 2 | 5 9 1 | 6 7 4
B| 4 9 6 | 3 8 7 | 2 5 1
C| 5 7 1 | 2 6 4 | 9 8 3
-+-----------------------
D| 1 8 5 | 7 4 6 | 3 9 2
E| 2 6 7 | 9 5 3 | 4 1 8
F| 9 4 3 | 8 1 2 | 7 6 5
-+-----------------------
G| 7 1 4 | 6 3 8 | 5 2 9
H| 3 2 9 | 1 7 5 | 8 4 6
I| 6 5 8 | 4 2 9 | 1 3 7

Reglas:

  1. Suponga que todos los laberintos se pueden resolver solo por lógica.
  2. Toda entrada tendrá 81 caracteres de longitud. Los caracteres que faltan serán 0.
  3. Salida de la solución como una sola cadena.
  4. La "cuadrícula" puede almacenarse internamente como lo desee.
  5. La solución debe usar una solución de adivinación de fuerza bruta.
  6. Las soluciones deben resolverse dentro de un plazo razonable.

Ejemplo de E / S:

>sudoku.py "030001000006000050500000983080006302000050000903800060714000009020000800000400030"
832591674496387251571264983185746392267953418943812765714638529329175846658429137
snmcdonald
fuente
¿Cómo puede la entrada tener 27 caracteres? Tiene que tener 81 caracteres de largo: 9 filas x 9 columnas. Eso es lo que hace tu ejemplo también. Además, supongo que "los caracteres faltantes serán 0" significa que si el número de caracteres es menor que 81, ¿entonces los ceros van al final?
Jonathan M Davis
Oh espera. Los caracteres que faltan serán de 0 bits. Duh Esos son los que hay que adivinar. En cualquier caso, el número de caracteres debe ser 81, no 27.
Jonathan M Davis
8
parece que las reglas 5 y 6 son un poco conflictivas ...
seudónimo117

Respuestas:

11

k (72 bytes)

El crédito por esto va para Arthur Whitney, creador del lenguaje k.

p,:3/:_(p:9\:!81)%3
s:{*(,x)(,/{@[x;y;:;]'&21=x[&|/p[;y]=p]?!10}')/&~x}
skeevey
fuente
¡clásico! ¡Yo también iba a publicar esto!
nightTrevors
9

Python, 188 bytes

Esta es una versión más corta de mi presentación ganadora para CodeSprint Sudoku , modificada para entrada de línea de comando en lugar de stdin (según el OP):

def f(s):
 x=s.find('0')
 if x<0:print s;exit()
 [c in[(x-y)%9*(x/9^y/9)*(x/27^y/27|x%9/3^y%9/3)or s[y]for y in range(81)]or f(s[:x]+c+s[x+1:])for c in'%d'%5**18]
import sys
f(sys.argv[1])

Si está utilizando Python 2, '%d'%5**18puede reemplazarse `5**18`para ahorrar 3 bytes.

Para que funcione más rápido, puede reemplazarlo '%d'%5**18con cualquier permutación '123456789'a un costo de 1 byte.

Si desea que acepte la entrada en stdin, puede reemplazar import sys;f(sys.argv[1]) con f(raw_input()), reduciéndolo a 177 bytes .

def f(s):
 x=s.find('0')
 if x<0:print s;exit()
 [c in[(x-y)%9*(x/9^y/9)*(x/27^y/27|x%9/3^y%9/3)or s[y]for y in range(81)]or f(s[:x]+c+s[x+1:])for c in'%d'%5**18]
f(raw_input())

EDITAR: Aquí hay un enlace a un tutorial más detallado.

Cifrado
fuente
Muy buena solución.
primo
8

Python, 197 caracteres

def S(s):
 i=s.find('0')
 if i<0:print s;return
 for v in'123456789':
  if sum(v==s[j]and(i/9==j/9or i%9==j%9or(i%9/3==j%9/3and i/27==j/27))for j in range(81))==0:S(s[:i]+v+s[i+1:])
S(raw_input())
Keith Randall
fuente
6

Respuesta en D:

import std.algorithm;
import std.conv;
import std.ascii;
import std.exception;
import std.stdio;

void main(string[] args)
{
    enforce(args.length == 2, new Exception("Missing argument."));
    enforce(args[1].length == 81, new Exception("Invalid argument."));
    enforce(!canFind!((a){return !isDigit(to!dchar(a));})
                     (args[1]),
                      new Exception("Entire argument must be digits."));

    auto sudoku = new Sudoku(args[1]);
    sudoku.fillIn();

    writeln(sudoku);
}

class Sudoku
{
public:

    this(string str) nothrow
    {
        normal = new int[][](9, 9);

        for(size_t i = 0, k =0; i < 9; ++i)
        {
            for(size_t j = 0; j < 9; ++j)
                normal[i][j] = to!int(str[k++]) - '0';
        }

        reversed = new int*[][](9, 9);

        for(size_t i = 0; i < 9; ++i)
        {
            for(size_t j = 0; j < 9; ++j)
                reversed[j][i] = &normal[i][j];
        }

        boxes = new int*[][](9, 9);
        indexedBoxes = new int*[][][](9, 9);

        for(size_t boxRow = 0, boxNum = 0; boxRow < 3; ++boxRow)
        {
            for(size_t boxCol = 0; boxCol < 3; ++boxCol, ++boxNum)
            {
                for(size_t i = 3 * boxRow, square = 0; i < 3 * (boxRow + 1); ++i)
                {
                    for(size_t j = 3 * boxCol; j < 3 * (boxCol + 1); ++j)
                    {
                        boxes[boxNum][square++] = &normal[i][j];
                        indexedBoxes[i][j] = boxes[boxNum];
                    }
                }
            }
        }
    }

    void fillIn()
    {
        fillIn(0, 0);
    }

    @property bool valid()
    {
        assert(full);

        for(size_t i = 0; i < 9; ++i)
        {
            for(int n = 1; n < 10; ++n)
            {
                if(!canFind(normal[i], n) ||
                   !canFind!"*a == b"(reversed[i], n) ||
                   !canFind!"*a == b"(boxes[i], n))
                {
                    return false;
                }
            }
        }

        return true;
    }

    override string toString() const
    {
        char[81] retval;

        for(size_t i = 0, k =0; i < 9; ++i)
        {
            for(size_t j = 0; j < 9; ++j)
                retval[k++] = to!char(normal[i][j] + '0');
        }

        return to!string(retval);
    }

private:

    @property bool full()
    {
        for(size_t i = 0; i < 9; ++i)
        {
            if(canFind(normal[i], 0))
                return false;
        }

        return true;
    }

    bool fillIn(size_t row, size_t col)
    {
        if(row == 9)
            return valid;

        size_t nextRow = row;
        size_t nextCol = col + 1;

        if(nextCol == 9)
        {
            nextRow = row + 1;
            nextCol = 0;
        }

        if(normal[row][col] == 0)
        {
            for(int n = 1; n < 10; ++n)
            {
                if(canFind(normal[row], n) ||
                   canFind!"*a == b"(reversed[col], n) ||
                   canFind!"*a == b"(indexedBoxes[row][col], n))
                {
                    continue;
                }

                normal[row][col] = n;

                if(fillIn(nextRow, nextCol))
                    return true;
            }

            normal[row][col] = 0;

            return false;
        }
        else
            return fillIn(nextRow, nextCol);
    }

    int[][] normal;
    int*[][] reversed;
    int*[][] boxes;
    int*[][][] indexedBoxes;
}

Con la entrada de muestra, toma 0.033 s en mi Phenom II X6 1090T cuando se compila dmd -w(es decir, sin optimizaciones), y toma 0.011 s cuando se compila dmd -w -O -inline -release(es decir, con optimizaciones).

Jonathan M Davis
fuente
4

J, 103

'p n'=:(;#)I.0=a=:("."0)Y
((a p}~3 :'>:?n#9')^:([:(27~:[:+/[:(9=#@~.)"1[:,/(2 2$3),;.3],|:,])9 9$])^:_)a

tiempo de ejecución esperado: O (miles de millones de años)

Eelvex
fuente
1
¿Y por qué el tiempo de ejecución esperado es "O (miles de millones de años)"? (¿no debería ser solo "miles de millones de años" sin el O?
Justin
1
Cuando vi esta pregunta, inmediatamente supe que J iba a aplastarla. Tiene que haber una manera de hacer esto más corto que K.
koko
1
@Quincunx, estrictamente hablando, es un uso incorrecto de big-O; se suponía que el "chiste" decía "tiempo de ejecución constante, miles de millones de años asintóticamente".
Eelvex
@koko, no pude encontrar algo mejor pero todavía estoy trabajando en ello.
Eelvex
4

Perl, 120 bytes

Ah, recuerdo jugar al golf en 2008 ... Y de hecho dejó de funcionar en Perl 5.12 ya que la configuración implícita de @_ por división se eliminó en ese momento. Así que solo prueba esto en un perl suficientemente viejo.

Ejecutar con la entrada en STDIN:

sudoku.pl <<< "030001000006000050500000983080006302000050000903800060714000009020000800000400030"

sudoku.pl:

${/[@_[map{$i-($i="@-")%9+$_,9*$_+$i%9,9*$_%26+$i-$i%3+$i%9-$i%27}0..8%split""]]/o||do$0}for$_=$`.$_.$'.<>,/0/||print..9
Ton Hospel
fuente
2
Es la tercera ley de Clarke , ¡pero a la inversa!
Conor O'Brien
3

Perl, 235 caracteres

$_=$s=<>;$r=join$/,map{$n=$_;'.*(?!'.(join'|',map+($_%9==$n%9||int($_/9)==int($n/9)||int($_/27)==int($n/27)&&int($_/3%3)==int($n/3%3)and$_<$n?'\\'.($_+1):$_>$n&&substr$s,$_,1)||X,@a).')(.).*'}@a=0..80;s!.!($&||123456789).$/!eg;say/^$r/

Esta es una versión de golf de algo que publiqué hace muchos años en la lista de correo Fun With Perl : una expresión regular para resolver sudoku.

Básicamente, divide la entrada en 81 líneas, cada una de las cuales contiene todos los números que podrían aparecer en el cuadrado correspondiente. Luego construye una expresión regular para que coincida con un número de cada línea, utilizando referencias y afirmaciones negativas anticipadas para rechazar soluciones que violen las restricciones de fila, columna o región. Luego hace coincidir la cadena con la expresión regular, permitiendo que el motor de expresión regular de Perl haga el trabajo duro de prueba y retroceso.

Sorprendentemente, es posible crear una expresión regular única que funcione para cualquier entrada, como lo hace mi programa original. Desafortunadamente, es bastante lento, por lo que basé el código de golf aquí en la versión de datos codificados (que se encuentra más adelante en el hilo FWP ), que modifica la expresión regular para rechazar temprano cualquier solución que sepa que luego violará una restricción. Esto lo hace razonablemente rápido para sudokus de nivel fácil a moderado, aunque los particularmente difíciles aún pueden tardar bastante tiempo en resolverse.

Ejecute el código con perl -M5.010para habilitar la función Perl 5.10+ say. La entrada debe darse en la entrada estándar, y la solución se imprimirá en la salida estándar; ejemplo:

$ perl -M5.010 golf/sudoku.pl
030001000006000050500000983080006302000050000903800060714000009020000800000400030
832591674496387251571264983185746392267953418943812765714638529329175846658429137
Ilmari Karonen
fuente
2

Guión de café de 1 línea

solve = (s, c = 0) -> if c is 81 then s else if s[x = c/9|0][y = c%9] isnt 0 then solve s, c+1 else (([1..9].filter (g) -> ![0...9].some (i) -> g in [s[x][i], s[i][y], s[3*(x/3|0) + i/3|0][3*(y/3|0) + i%3]]).some (g) -> s[x][y] = g; solve s, c+1) or s[x][y] = 0

Aquí está la versión más grande con uso de muestra :

solve = (sudoku, cell = 0) ->
  if cell is 9*9 then return sudoku

  x = cell%9
  y = (cell - x)/9

  if sudoku[x][y] isnt 0 then return solve sudoku, cell+1

  row = (i) -> sudoku[x][i]
  col = (i) -> sudoku[i][y]
  box = (i) -> sudoku[x - x%3 + (i - i%3)/3][y - y%3 + i%3]

  good = (guess) -> [0...9].every (i) -> guess not in [row(i), col(i), box(i)]

  guesses = [1..9].filter good

  solves = (guess) -> sudoku[x][y] = guess; solve sudoku, cell+1

  (guesses.some solves) or sudoku[x][y] = 0

sudoku = [
  [1,0,0,0,0,7,0,9,0],
  [0,3,0,0,2,0,0,0,8],
  [0,0,9,6,0,0,5,0,0],
  [0,0,5,3,0,0,9,0,0],
  [0,1,0,0,8,0,0,0,2],
  [6,0,0,0,0,4,0,0,0],
  [3,0,0,0,0,0,0,1,0],
  [0,4,0,0,0,0,0,0,7],
  [0,0,7,0,0,0,3,0,0]
]
console.log if solve sudoku then sudoku else 'could not solve'
pathikrit
fuente
1
Podría acortarse acortando solve, eliminando muchos espacios en blanco (sé que es significativo, pero en muchos lugares podría eliminarse), usando símbolos en lugar de palabras (como en !=lugar de isnt), usando sangría en lugar de thenpalabra clave, reemplazando [0...9]con [0..8].
Konrad Borowski el
1

Clojure - 480 bytes

El tamaño explotó, pero al menos es un número bonito. Creo que podría mejorarse mucho usando solo 1D-vector. De todos modos, el caso de prueba tarda un poco menos de cuatro segundos en mi computadora portátil. Pensé que sería apropiado definir una función, ya que es un lenguaje funcional después de todo.

(defn f[o &[x y]](if x(if(> y 8)(apply str(map #(apply str %)o))(first(for[q[(o y)]v(if(=(q x)0)(range 1 10)[(q x)])d[(assoc o y(assoc(o y)x v))]s[(and(every? true?(concat(for[i(range 9)](and(or(not=((d y)i)v)(= i x))(or(not=((d i)x)v)(= i y))))(for[m[#(+ %2(- %(mod % 3)))]r[(range 3)]a r b r c[(m y b)]e[(m x a)]](or(and(= e x)(= c y))(not=((d y)x)((d c)e))))))(f d(mod(+ x 1)9)(if(= x 8)(+ 1 y)y)))]:when s]s)))(f(vec(for[a(partition 9 o)](vec(map #(Integer.(str %))a))))0 0)))

Ejemplos:

(f "030001000006000050500000983080006302000050000903800060714000009020000800000400030")
=> "832591674496387251571264983185746392267953418943812765714638529329175846658429137"
(f "004720900039008005001506004040010520028050170016030090400901300100300840007085600")
=> "654723981239148765871596234743819526928654173516237498482961357165372849397485612"

Una versión ungolfing (y más bonita):

(defn check-place [o x y v]
  (and (every? true? (for [i (range 9)]
                       (and (or (not= ((o y) i) v) (= i x))
                            (or (not= ((o i) x) v) (= i y)))))
       (every? true?
         (for [r [(range 3)]
               a r
               b r
               c [(+ b (- y (mod y 3)))]
               d [(+ a (- x (mod x 3)))]]
           (or (and (= d x) (= c y)) (not= ((o y) x) ((o c) d)))))))

(defn solve-sudoku [board & [x y]]
  (if x
    (if (> y 8)
      (apply str (map #(apply str %) board))
      (first
        (for [v (if (= ((board y) x) 0) (range 1 10) [((board y) x)])
              :let [a (mod (+ x 1) 9)
                    b (if (= x 8) (+ 1 y) y)
                    d (assoc board y (assoc (board y) x v))
                    s (and (check-place d x y v) (solve-sudoku d a b))]
              :when s]
          s)))
    (solve-sudoku (vec (for [a (partition 9 board)]
                         (vec (map #(Integer. (str %)) a)))) 0 0)))
seequ
fuente
1

PowerShell , 244 242 218 215 bytes

$a=(24,24,6)*3|%{,(0..8|%{($r++)});,(0..8|%{$c%81;$c+=9});$c++;,((1,1,7)*3|%{+$q;$q+=$_});$q-=$_}
$f={param($s)$l,$r=$s-split0,2;if($p=$a|?{$l.length-in$_}){1..9|?{"$_"-notin($p|%{$s[$_]})}|%{&$f "$l$_$r"}}else{$s}}

Pruébalo en línea!

El script encuentra todas las soluciones para un sudoku.

Desenrollado:

$a=(24,24,6)*3|%{                       # array of indexes for a sudoku...
    ,(0..8|%{($r++)})                   # rows
    ,(0..8|%{$c%81;$c+=9});$c++         # columns
    ,((1,1,7)*3|%{+$q;$q+=$_});$q-=$_   # and squares
}

$f = {
    param($s)

    # optional log. remove this statement in a release version.
    if($script:iter++ -lt 100 -or ($script:iter%100)-eq0){
        Write-Information ('{0}: {1,6}: {2}'-f (get-Date), $script:iter, ($s-replace0,' ')) -InformationAction Continue
    }

    $left,$right=$s-split0,2                # split by a first 0; $left.length is a position of this 0 if $s contains the 0
    if( $parts=$a|?{$left.length-in$_} ){   # get sudoku parts (rows, columns, squares) contain the position
        1..9|?{                             # try a digit
            "$_"-notin($parts|%{$s[$_]})    # all digits in these parts will be unique if parts do not contain the digit
        }|%{
            &$f "$left$_$right"             # recursive call with the digit
        } #|select -f 1                     # uncomment this to get a first result only
    }
    else{
        $s
    }

}

Casos de prueba:

@(
    # 5 iterations, my notebook: 00:00:00, all
    # 5 iterations, my notebook: 00:00:00, first only
    , ( "832591674496387251571264983185746392267953418943812765714638529329175846658400030",
        "832591674496387251571264983185746392267953418943812765714638529329175846658429137" )

    # ~29600 iterations, my notebook: 00:01:27, all
    #  ~2100 iterations, my notebook: 00:00:10, first only
    # , ( "830001000006000050500000983080006302000050000903800060714000009020000800000400030",
    #     "832591674496387251571264983185746392267953418943812765714638529329175846658429137" )

    # ~49900 iterations, my notebook: 00:02:39, all
    # ~22400 iterations, my notebook: 00:01:20, first only
    # , ( "030001000006000050500000983080006302000050000903800060714000009020000800000400030",
    #     "832591674496387251571264983185746392267953418943812765714638529329175846658429137" )

) | % {
    $sudoku, $expected = $_
    $time = Measure-Command {
        $result = &$f $sudoku
    }
    "$($result-contains$expected): $time"
    $result
}
mazzy
fuente
0

D (322 caracteres)

Para cada cuadrado sin resolver, crea una serie de opciones disponibles y luego recorre el mismo.

import std.algorithm,std.range,std.stdio;void main(char[][]args){T s(T)(T p){foreach(i,ref c;p)if(c<49){foreach(o;"123456789".setDifference(chain(p[i/9*9..i/9*9+9],p[i%9..$].stride(9),p[i/27*27+i%9/3*3..$][0..21].chunks(3).stride(3).joiner).array.sort)){c=o&63;if(s(p))return p;}c=48;return[];}return p;}s(args[1]).write;}

con espacios en blanco:

import std.algorithm, std.range, std.stdio;

void main(char[][] args) {
    T s(T)(T p) {
        foreach (i, ref c; p) if (c < 49) {
            foreach (o; "123456789".setDifference(chain(
                    p[i/9*9..i/9*9+9],
                    p[i%9..$].stride(9),
                    p[i/27*27+i%9/3*3..$][0..21].chunks(3).stride(3).joiner
                ).array.sort))
            {
                c = o&63;
                if (s(p)) return p;
            }
            c=48;
            return [];
        }
        return p;
    }
    s(args[1]).write;
}
mleise
fuente
0

Perl (195 caracteres)

use integer;@A=split//,<>;sub R{for$i(0..80){next if$A[$i];my%t=map{$_/9==$/9||$_%9==$i%9||$_/27==$i/27&&$_%9/3==$i%9/3?$A[$_]:0=>1}0..80;R($A[$i]=$_)for grep{!$t{$_}}1..9;return$A[$i]=0}die@A}R

Todo el crédito va al creador aquí , y la explicación también se puede encontrar allí.

Qwix
fuente
1
Si no lo escribió usted mismo, debe verificar el botón "Wiki de la comunidad".
Kyle Kanos
Después de investigar un poco sobre qué es eso, no me parece posible. Aparentemente se necesitan 100 repeticiones para que yo vea la casilla de verificación (vea la sección de apéndice en el n. ° 2 de esta publicación )
Qwix
Hmm, no estaba al tanto de ese requisito.
Kyle Kanos
0

J 94 bytes

Funciona exactamente de la misma manera que la versión K, es decir, con un BFS (por lo que generará todas las soluciones). Imprime espacios entre los dígitos de salida, pero también lo hace el programa K. No estoy contando "s =:" ya que esto es solo nombrar la función (al igual que no contaría el nombre de archivo en otro idioma).

   s=: [:<@((]i.0:)}"0 _~(>:i.9)-.{&((+./ .=|:)3(],.[,@#.<.@%~)9 9#:i.81)@i.&0#])"1^:(0 e.,)@;^:_"."0

   s'030001000006000050500000983080006302000050000903800060714000009020000800000400030'
8 3 2 5 9 1 6 7 4 4 9 6 3 8 7 2 5 1 5 7 1 2 6 4 9 8 3 1 8 5 7 4 6 3 9 2 2 6 7 9 5 3 4 1 8 9 4 3 8 1 2 7 6 5 7 1 4 6 3 8 5 2 9 3 2 9 1 7 5 8 4 6 6 5 8 4 2 9 1 3 7
Olius
fuente