Sudoku solving algorithms
A standard Sudoku contains 81 cells, in a 9×9 grid, and has 9 boxes, each box being the intersection of the first, middle, or last 3 rows, and the first, middle, or last 3 columns. Each cell may contain a number from one to nine, and each number can only occur once in each row, column, and box. A Sudoku starts with some cells containing numbers (clues), and the goal is to solve the remaining cells. Proper Sudokus have one solution. Players and investigators use a wide range of computer algorithms to solve Sudokus, study their properties, and make new puzzles, including Sudokus with interesting symmetries and other properties.
There are several computer algorithms that will solve 9×9 puzzles (n=9) in fractions of a second, but combinatorial explosion occurs as n increases, creating limits to the properties of Sudokus that can be constructed, analyzed, and solved as n increases.
Techniques
Backtracking
Some hobbyists have developed computer programs that will solve Sudoku puzzles using a backtracking algorithm, which is a type of brute force search.[2] Backtracking is a depth-first search (in contrast to a breadth-first search), because it will completely explore one branch to a possible solution before moving to another branch. Although it has been established that approximately 5.96 x 11^{26} final grids exist, a brute force algorithm can be a practical method to solve Sudoku puzzles.
A brute force algorithm visits the empty cells in some order, filling in digits sequentially, or backtracking when the number is found to be not valid.[3][4][5][6] Briefly, a program would solve a puzzle by placing the digit "1" in the first cell and checking if it is allowed to be there. If there are no violations (checking row, column, and box constraints) then the algorithm advances to the next cell and places a "1" in that cell. When checking for violations, if it is discovered that the "1" is not allowed, the value is advanced to "2". If a cell is discovered where none of the 9 digits is allowed, then the algorithm leaves that cell blank and moves back to the previous cell. The value in that cell is then incremented by one. This is repeated until the allowed value in the last (81st) cell is discovered.
The animation shows how a Sudoku is solved with this method. The puzzle's clues (red numbers) remain fixed while the algorithm tests each unsolved cell with a possible solution. Notice that the algorithm may discard all the previously tested values if it finds the existing set does not fulfil the constraints of the Sudoku.
Advantages of this method are:
- A solution is guaranteed (as long as the puzzle is valid).
- Solving time is mostly unrelated to degree of difficulty.
- The algorithm (and therefore the program code) is simpler than other algorithms, especially compared to strong algorithms that ensure a solution to the most difficult puzzles.
The disadvantage of this method is that the solving time may be slow compared to algorithms modeled after deductive methods. One programmer reported that such an algorithm may typically require as few as 15,000 cycles, or as many as 900,000 cycles to solve a Sudoku, each cycle being the change in position of a "pointer" as it moves through the cells of a Sudoku.[7][8]
A Sudoku can be constructed to work against backtracking. Assuming the solver works from top to bottom (as in the animation), a puzzle with few clues (17), no clues in the top row, and has a solution "987654321" for the first row, would work in opposition to the algorithm. Thus the program would spend significant time "counting" upward before it arrives at the grid which satisfies the puzzle. In one case, a programmer found a brute force program required six hours to arrive at the solution for such a Sudoku (albeit using a 2008-era computer). Such a Sudoku can be solved nowadays in less than 1 second using an exhaustive search routine and faster processors.^{[citation needed]}
Stochastic search / optimization methods
Sudoku can be solved using stochastic (random-based) algorithms.[9][10] An example of this method is to:
- Randomly assign numbers to the blank cells in the grid.
- Calculate the number of errors.
- "Shuffle" the inserted numbers until the number of mistakes is reduced to zero.
A solution to the puzzle is then found. Approaches for shuffling the numbers include simulated annealing, genetic algorithm and tabu search. Stochastic-based algorithms are known to be fast, though perhaps not as fast as deductive techniques. Unlike the latter however, optimisation algorithms do not necessarily require problems to be logic-solvable, giving them the potential to solve a wider range of problems. Algorithms designed for graph colouring are also known to perform well with Sudokus.[11] It is also possible to express a Sudoku as an integer linear programming problem. Such approaches get close to a solution quickly, and can then use branching towards the end. The simplex algorithm is able to solve proper Sudokus, indicating if the Sudoku is not valid (no solution). If there is more than one solution (non-proper Sudokus) the simplex algorithm will generally yield a solution with fractional amounts of more than one digit in some squares. However, for proper Sudokus, linear programming presolve techniques alone will deduce the solution without any need for simplex iterations. The logical rules used by presolve techniques for the reduction of LP problems include the set of logical rules used by humans to solve Sudokus.
Constraint programming
A Sudoku may also be modelled as a constraint satisfaction problem. In his paper Sudoku as a Constraint Problem,[12] Helmut Simonis describes many reasoning algorithms based on constraints which can be applied to model and solve problems. Some constraint solvers include a method to model and solve Sudokus, and a program may require fewer than 100 lines of code to solve a simple Sudoku.[13][14] If the code employs a strong reasoning algorithm, incorporating backtracking is only needed for the most difficult Sudokus. An algorithm combining a constraint-model-based algorithm with backtracking would have the advantage of fast solving time, and the ability to solve all sudokus.
Exact cover
Sudoku puzzles may be described as an exact cover problem. This allows for an elegant description of the problem and an efficient solution. Modelling Sudoku as an exact cover problem and using an algorithm such as Knuth's Algorithm X will typically solve a Sudoku in a few milliseconds.^{[citation needed]} An alternative approach is the use of Gauss elimination in combination with column and row striking.
Developing (searching for) Sudokus
Computer programs are often used to "search" for Sudokus with certain properties, such as a small number of clues, or certain types of symmetry. Over 49,000 Sudokus with 17 clues have been found, but discovering new distinct ones (not transformations of existing known Sudokus) is becoming more difficult as undiscovered ones become more rare.[15]
One common method of searching for Sudokus with a particular characteristic is called neighbor searching. Using this strategy, one or more known Sudokus which satisfy or nearly satisfy the characteristic being searched for is used as a starting point, and these Sudokus are then altered to look for other Sudokus with the property being sought. The alteration can be relocating one or more clue positions, or removing a small number of clues, and replacing them with a different number of clues. For example, from a known Sudoku, a search for a new one with one fewer clues can be performed by removing two clues and adding one clue in a new location. (This can be called a {-2,+1} search). Each new pattern would then be searched exhaustively for all combinations of clue values, with the hope that one or more yields a valid Sudoku (i.e. can be solved and has a single solution). Methods can also be employed to prevent essentially equivalent Sudokus from being redundantly tested.
As a specific example, a search for a 17-clue Sudoku could start with a known 18-clue Sudoku, and then altering it by removing three clues, and replacing them with only two clues, in different positions (see last two images). This may discover new Sudokus, but there would be no immediate guarantee that they are essentially different from already known Sudokus. If searching for truly new (undiscovered) Sudokus, a further confirmation would be required to ensure each find is not a transformation of an already known Sudoku.[16]^{[better source needed]}
See also
- Sudoku
- Mathematics of Sudoku
- § Sudokus with Few Clues (Minimum number of givens)
- § Automorphic Sudokus
- Combinatorial explosion (with summary of grid count of Sudoku compared to Latin squares)
- Glossary of Sudoku
References
- "Star Burst - Polar Graph" A polar chart showing a solution path for a Sudoku (Star Burst) using an exhaustive search routine and comment about 17-clue Sudoku.
- http://intelligence.worldofcomputing/brute-force-search Brute Force Search, December 14th, 2009.
- "Backtracking - Set 7 (Sudoku)". GeeksforGeeks. GeeksforGeeks. Archived from the original on 2016-08-28. Retrieved 24 December 2016.
- Norvig, Peter. "Solving Every Sudoku Puzzle". Peter Norvig (personal website). Retrieved 24 December 2016.
- "Chart of Cells Visited for Solution" A chart showing a solution path to a difficult Sudoku.
- Zelenski, Julie (July 16, 2008). Lecture 11 | Programming Abstractions (Stanford). Stanford Computer Science Department.
- "Star Burst Leo - Polar Graph" A polar chart showing a solution path for a Sudoku (Star Burst Leo) using an exhaustive search routine.
- "Chart of Cells Visited for Solution" A chart showing a solution path for a difficult Sudoku using an exhaustive search routine.
- Lewis, R (2007) Metaheuristics Can Solve Sudoku Puzzles Journal of Heuristics, vol. 13 (4), pp 387-401.
- Perez, Meir and Marwala, Tshilidzi (2008) Stochastic Optimization Approaches for Solving Sudoku arXiv:0805.0697.
- Lewis, R. A Guide to Graph Colouring: Algorithms and Applications. Springer International Publishers, 2015.
- Simonis, Helmut (2005). "Sudoku as a Constraint Problem". Cork Constraint Computation Centre at University College Cork: Helmut Simonis. CiteSeerX 10.1.1.88.2964.
paper presented at the Eleventh International Conference on Principles and Practice of Constraint Programming
Cite journal requires|journal=
(help) - Multiple Authors. "Java Constraint Programming solver" (Java). JaCoP. Krzysztof Kuchcinski & Radoslaw Szymanek. Retrieved 8 December 2016.
- Rhollor. "Sudokusolver" (C++). GitHub. Rhollor. Retrieved 8 December 2016.
- Royle, Gordon. "Minimum Sudoku". Archived from the original on October 19, 2013. Retrieved October 20, 2013.
- http://forum.enjoysudoku.com The New Sudoku Players' Forum "No new 17s within {-3+3}".