Implementa un Risolutore di Sudoku Brute Force


20

Implementa il solutore Sudoku più breve usando le ipotesi. Da quando ho ricevuto alcune richieste ho aggiunto questa come una domanda alternativa per coloro che desiderano implementare un risolutore di sudoku a forza 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

Risposta:

 | 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

Regole:

  1. Supponiamo che tutti i labirinti siano risolvibili solo con la logica.
  2. Tutti gli input avranno una lunghezza di 81 caratteri. I caratteri mancanti saranno 0.
  3. Invia la soluzione come una singola stringa.
  4. La "griglia" può essere memorizzata internamente come desideri.
  5. La soluzione deve utilizzare una soluzione per indovinare la forza bruta.
  6. Le soluzioni dovrebbero risolversi entro un termine ragionevole.

Esempio I / O:

>sudoku.py "030001000006000050500000983080006302000050000903800060714000009020000800000400030"
832591674496387251571264983185746392267953418943812765714638529329175846658429137

In che modo l'input può contenere 27 caratteri? Deve essere lungo 81 caratteri - 9 righe x 9 colonne. Questo è anche il tuo esempio. Inoltre, suppongo che "i caratteri mancanti saranno 0" significa che se il numero di caratteri è inferiore a 81, gli zeri vanno alla fine?
Jonathan M Davis,

Oh aspetta. Ottengo che i caratteri mancanti saranno 0 bit. Duh. Quelli sono quelli che devono essere indovinati. In ogni caso, il numero di personaggi deve essere 81, non 27.
Jonathan M Davis

8
sembra che le regole 5 e 6 siano in conflitto ...
pseudonimo117

Risposte:


11

k (72 byte)

Il merito è di Arthur Whitney, ideatore della lingua k.

p,:3/:_(p:9\:!81)%3
s:{*(,x)(,/{@[x;y;:;]'&21=x[&|/p[;y]=p]?!10}')/&~x}

classico! Avrei pubblicato anche questo!
nightTrevors,

9

Python, 188 byte

Questa è un'ulteriore versione abbreviata della mia proposta vincente per CodeSprint Sudoku , modificata per l'input da riga di comando anziché da stdin (come da 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])

Se stai usando Python 2, '%d'%5**18può essere sostituito con `5**18`per salvare 3 byte.

Per farlo funzionare più velocemente, puoi sostituirlo '%d'%5**18con qualsiasi permutazione di'123456789' al costo di 1 byte.

Se invece vuoi accettare l'input su stdin, puoi sostituirlo import sys;f(sys.argv[1])con f(raw_input()), portandolo a 177 byte .

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())

EDIT: ecco un link a una procedura dettagliata più dettagliata.


Soluzione molto bella.
primo

8

Python, 197 caratteri

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())

6

Risposta in 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 l'input di esempio, sono necessari .033s sul mio Phenom II X6 1090T quando compilato dmd -w(cioè senza ottimizzazioni), e sono necessari .011s quando compilato dmd -w -O -inline -release(cioè con ottimizzazioni).


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

tempo di esecuzione previsto: O (gazillion miliardi di anni)


1
E perché è il tempo di esecuzione previsto "O (gazillion miliardi di anni)"? (non dovrebbe essere solo "gazillion miliardi di anni" senza la O?
Giustino,

1
Quando ho visto questa domanda ho capito subito che J avrebbe schiacciato questo. Deve esserci un modo per renderlo più breve di K.
koko,

1
@Quincunx, a rigor di termini, è un uso errato di big-O; lo "scherzo" doveva leggere "tempo di esecuzione costante, asintoticamente gazillion miliardi di anni".
Eelvex,

@koko, non sono riuscito a trovare qualcosa di meglio ma ci sto ancora lavorando.
Eelvex,

4

Perl, 120 byte

Oh, mi ricordo di aver giocato a golf nel 2008 ... E in effetti ha smesso di funzionare in perl 5.12 da quando è stata rimossa l'impostazione implicita di @_ per split. Quindi provalo solo su un perl sufficientemente vecchio.

Esegui con l'input su 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

2
È la terza legge di Clarke , ma al contrario!
Conor O'Brien,

3

Perl, 235 caratteri

$_=$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/

Questa è una versione da golf di qualcosa che ho pubblicato molti anni fa nella mailing list di Fun With Perl : una regexp per la risoluzione di sudoku.

Fondamentalmente, manipola l'input in 81 righe, ognuna contenente tutti i numeri che potrebbero verificarsi nel quadrato corrispondente. Costruisce quindi una regexp in modo che corrisponda a un numero per ciascuna riga, utilizzando riferimenti indietro e asserzioni di lookahead negative per rifiutare soluzioni che violano i vincoli di riga, colonna o area. Quindi abbina la stringa al regexp, lasciando che il motore regexp di Perl faccia il duro lavoro di prova e backtracking.

Sorprendentemente, è possibile creare un unico regexp che funzioni per qualsiasi input, come fa il mio programma originale. Sfortunatamente, è piuttosto lento, quindi ho basato qui il codice golfed sulla versione hardcoded-givens ( trovata più avanti nel thread FWP ), che modifica la regexp per rifiutare in anticipo qualsiasi soluzione che sappia in seguito violerebbe un vincolo. Ciò lo rende ragionevolmente veloce per i sudoku di livello da facile a moderato, anche se quelli particolarmente difficili possono ancora richiedere molto tempo per essere risolti.

Esegui il codice con perl -M5.010per abilitare la sayfunzione Perl 5.10+ . L'input deve essere fornito su input standard e la soluzione verrà stampata su output standard; esempio:

$ perl -M5.010 golf/sudoku.pl
030001000006000050500000983080006302000050000903800060714000009020000800000400030
832591674496387251571264983185746392267953418943812765714638529329175846658429137

2

1 tazza di caffè con sceneggiatura

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

Ecco la versione più grande con utilizzo di esempio :

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'

1
Potrebbe essere accorciato accorciando solve, rimuovendo un sacco di spazi bianchi (so che è significativo, ma in molti punti potrebbe essere rimosso), usando simboli invece di parole (come !=invece di isnt), usando rientro anziché thenparola chiave, sostituendo [0...9]con [0..8].
Konrad Borowski il

1

Clojure - 480 byte

Le dimensioni sono esplose, ma almeno è un bel numero. Penso che potrebbe essere migliorato molto usando solo il vettore 1D. Ad ogni modo, il test case impiega poco meno di quattro secondi sul mio laptop. Ho pensato che sarebbe opportuno definire una funzione, poiché dopo tutto è un linguaggio funzionale.

(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)))

Esempi:

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

Una versione leggermente non golfata (e più carina):

(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)))

1

PowerShell , 244 242 218 215 byte

$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}}

Provalo online!

Lo script trova tutte le soluzioni per un sudoku.

srotolato:

$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
    }

}

Casi test:

@(
    # 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
}

0

D (322 caratteri)

Per ogni quadrato irrisolto crea una serie di opzioni disponibili e quindi ci passa sopra.

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 spazi bianchi:

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;
}

0

Perl (195 caratteri)

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

Tutto il merito va al creatore qui e la spiegazione può essere trovata anche lì.


1
Se non l'hai scritto tu stesso, allora dovresti controllare il pulsante "Community Wiki".
Kyle Kanos,

Dopo alcune ricerche su cosa sia, non sembra possibile per me. A quanto pare sono necessari 100 rappresentanti per visualizzare la casella di controllo (vedere la sezione addendum al numero 2 di questo post )
Qwix,

Hmm, non ero a conoscenza di quel requisito.
Kyle Kanos,

0

J, 94 byte

Funziona esattamente allo stesso modo della versione K, in particolare con un BFS (quindi deve generare tutte le soluzioni). Stampa gli spazi tra le cifre di output, ma anche il programma K. Non sto contando "s =:" poiché si tratta solo di nominare la funzione (proprio come se non contassi il nome del file in un'altra lingua).

   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
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.