I am working on an optimized Sudoku solver. The basic idea of backtracking remains, though I have added forward checking. The domain (what numbers can be set into a cell) is stored for every cell. Before I try to set a number into a cell and backtrack if it doesn’t work, I check all the related cells to see if the number is the only one that could be set there. If a conflict occurs, I try the next number, and so on. Optimally, this reduces the amount of backtracking necessary.
However, my code is unfortunately not working properly. I would greatly appreciate your help.
Best regards
Eric
Code (Comments are in German):
#include <iostream>
#include <bitset>
#include <array>
using namespace std;
#define n 9
void print_board(int board[n][n]);
bool solve(array<bitset<n>, n * n> &domains, int board[n][n], int i = 0, int j = 0);
void initialize_domains(array<bitset<n>, n * n> &domains, int board[n][n]);
bool update_domains(array<bitset<n>, n * n> &domains, int board[n][n], int i, int j, int k);
std::bitset<n> reverse_bits(const std::bitset<n> &bitset)
{
std::bitset<n> reversed;
for (int i = 0; i < n; ++i)
{
reversed[i] = bitset[n - 1 - i]; // Setzt das i-te Bit in reversed auf das (N-1-i)-te Bit in bitset
}
return reversed;
}
void printDomains(const std::array<std::bitset<n>, n * n> &domains)
{
for (int i = 0; i < n; ++i)
{
for (int j = 0; j < n; ++j)
{
std::cout << "Domaene fuer Zelle (" << i << ", " << j << "): " << reverse_bits(domains[i * n + j]) << std::endl;
}
}
}
int main() // hauptfunktion handelt den Ablauf des Programms
{
int board[n][n] = {// definition anfänglichen des Spielfeldes
{5, 3, 0, 0, 7, 0, 0, 0, 0},
{6, 0, 0, 1, 9, 5, 0, 0, 0},
{0, 9, 8, 0, 0, 0, 0, 6, 0},
{8, 0, 0, 0, 6, 0, 0, 0, 3},
{4, 0, 0, 8, 0, 3, 0, 0, 1},
{7, 0, 0, 0, 2, 0, 0, 0, 6},
{0, 6, 0, 0, 0, 0, 2, 8, 0},
{0, 0, 0, 4, 1, 9, 0, 0, 5},
{0, 0, 0, 0, 8, 0, 0, 7, 9}};
std::array<std::bitset<n>, n * n> domains;
initialize_domains(domains, board);
print_board(board);
solve(domains, board);
cout << "n";
print_board(board);
}
void print_board(int board[n][n]) // zeigt das Sudoku Feld an
{
for (int i = 0; i < n; ++i)
{
for (int j = 0; j < n; ++j)
{
cout << board[i][j] << " ";
}
cout << "n";
}
}
// ausgehend von der Zelle einer Zelle werden die Domänmen der beeinflussten Zellen aktualisiert
bool update_domains(std::array<std::bitset<n>, n * n> &domains, int board[n][n], int i, int j, int k)
{
for (int a = 0; a < n; ++a)
{
domains[i * n + a].set(k - 1); // reihe
domains[j + n * a].set(k - 1); // splate
// Wenn es nicht mehr möglich ist eine Zahl an die Stelle zu setzten
if (domains[i * n + a].all() && board[i][a] == 0)
return false;
if (domains[j + n * a].all() && board[a][j] == 0)
return false;
}
for (int row = (i / 3) * 3; row < ((i / 3) * 3) + 3; ++row)
{
for (int col = (j / 3) * 3; col < ((j / 3) * 3) + 3; ++col)
{
domains[row * n + col].set(k - 1); // Subgrid
if (domains[row * n + col].all() && board[row][col] == 0)
return false;
}
}
return true;
}
// 0 = möglich die Zahl zu setzten
// 1 = nicht möglich die Zahl zu setzten
void initialize_domains(array<bitset<n>, n * n> &domains, int board[n][n])
{
for (int i = 0; i < n; ++i)
{
for (int j = 0; j < n; ++j)
{
if (board[i][j] != 0)
update_domains(domains, board, i, j, board[i][j]);
}
}
}
// i gibt die Zeile an
// j gibt die Spalte an
bool solve(array<bitset<n>, n * n> &domains, int board[n][n], int i, int j) // Ist für das lösen des Sudokos durch Bruteforce Backtracking tätig
{
if (i == n) // wenn alle Zeilen und damit auch alle Splaten erfolgreich ausgefüllt wurden, ist das Sudoku Spiel gelößt
return true;
if (j == n) // Wenn die Splatenangabe über die Diemension geht, muss die erste Zelle der nächsten Zeile betrachtet werden
return solve(domains, board, i + 1, 0);
if (board[i][j] != 0) // Wenn keine leere Zelle an der aktullen Position vorhanden ist, wird die nächste Position betrachtet
return solve(domains, board, i, j + 1);
for (int k = 1; k <= n; ++k) // Alle Validen Zahlen von 1-9 werden linear ausprobiert
{
if (domains[i * n + j].test(k - 1)) // Wenn die Zahl nicht in der Domäne der aktuellen Zelle vorhanden ist, kann sie übersprungen werden
continue;
array<bitset<n>, n * n> pre_updated_domain = domains;
if (update_domains(domains, board, i, j, k)) // Kein Konfilkt
{
board[i][j] = k; // Die vermeintlich valide Zahl wird ausprobiert
if (solve(domains, board, i, j + 1)) // kommen keine weiteren Konfilkte auf, ist dies die richtie lösung
return true;
board[i][j] = 0;
}
// Wenn die Zahl zu einem Konfilkt führt (in einer beeinflussten Zelle, kann keine Zahl mehr eingesetzt werden) wird diese nicht verwendet
domains = pre_updated_domain;
}
return false; // Falls keine der n Zahlen valide ist, ist ein Konfilkt aufgetreten
}
My current is outcome i the following:
5 3 0 0 7 0 0 0 0
6 0 0 1 9 5 0 0 0
0 9 8 0 0 0 0 6 0
8 0 0 0 6 0 0 0 3
4 0 0 8 0 3 0 0 1
7 0 0 0 2 0 0 0 6
0 6 0 0 0 0 2 8 0
0 0 0 4 1 9 0 0 5
0 0 0 0 8 0 0 7 9
5 3 0 0 7 0 0 0 0
6 0 0 1 9 5 0 0 0
0 9 8 0 0 0 0 6 0
8 0 0 0 6 0 0 0 3
4 0 0 8 0 3 0 0 1
7 0 0 0 2 0 0 0 6
0 6 0 0 0 0 2 8 0
0 0 0 4 1 9 0 0 5
0 0 0 0 8 0 0 7 9
I have expected:
5 3 4 6 7 8 9 1 2
6 7 2 1 9 5 3 4 8
1 9 8 3 4 2 5 6 7
8 5 9 7 6 1 4 2 3
4 2 6 8 5 3 7 9 1
7 1 3 9 2 4 8 5 6
9 6 1 5 3 7 2 8 4
2 8 7 4 1 9 6 3 5
3 4 5 2 8 6 1 7 9
Eric is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.