Python中最短的数独解算器-它是如何工作的?

ymzxtsji  于 2023-02-21  发布在  Python
关注(0)|答案(5)|浏览(120)

我正在玩我自己的数独解算器,并正在寻找一些指针,以良好和快速的设计时,我遇到了这个:

def r(a):i=a.find('0');~i or exit(a);[m
in[(i-j)%9*(i/9^j/9)*(i/27^j/27|i%9/3^j%9/3)or a[j]for
j in range(81)]or r(a[:i]+m+a[i+1:])for m in'%d'%5**18]
from sys import*;r(argv[1])

我自己的实现解决数独问题的方式和我在脑海中解决它们的方式一样,但这个神秘的算法是如何工作的呢?
http://scottkirkwood.blogspot.com/2006/07/shortest-sudoku-solver-in-python.html

sxpgvts3

sxpgvts31#

那么,您可以通过修改语法来使事情变得简单一些:

def r(a):
  i = a.find('0')
  ~i or exit(a)
  [m in[(i-j)%9*(i/9^j/9)*(i/27^j/27|i%9/3^j%9/3)or a[j]for j in range(81)] or r(a[:i]+m+a[i+1:])for m in'%d'%5**18]
from sys import *
r(argv[1])

稍微清理一下:

from sys import exit, argv
def r(a):
  i = a.find('0')
  if i == -1:
    exit(a)
  for m in '%d' % 5**18:
    m in[(i-j)%9*(i/9^j/9)*(i/27^j/27|i%9/3^j%9/3) or a[j] for j in range(81)] or r(a[:i]+m+a[i+1:])

r(argv[1])

好的,这个脚本需要一个命令行参数,并在其上调用函数r,如果字符串中没有零,r退出并打印出它的参数。
(If传递另一种类型的对象,None相当于传递零,任何其他对象都将打印到sys.stderr,并导致退出代码1。特别是,sys.exit(“某个错误消息”)是发生错误时退出程序的快速方法。请参见http://www.python.org/doc/2.5.2/lib/module-sys.html
我想这意味着零对应于开放空间,一个没有零的谜题就解决了,然后就是那个讨厌的递归表达式。
这个循环很有趣:for m in'%d'%5**18
为什么是5**18呢?结果'%d'%5**18的计算结果是'3814697265625'。这是一个字符串,其中每个数字1-9至少出现一次,所以它可能试图放置它们中的每一个。实际上,看起来r(a[:i]+m+a[i+1:])是这样做的:递归地调用r,用字符串中的数字填充第一个空格,但这只在前面的表达式为假时发生,让我们来看一下:
m in [(i-j)%9*(i/9^j/9)*(i/27^j/27|i%9/3^j%9/3) or a[j] for j in range(81)]
因此,只有当m不在那个怪物列表中时,才能进行替换。每个元素要么是数字(如果第一个表达式为非零),要么是字符(如果第一个表达式为零)。如果m以字符的形式出现,就排除了它作为可能的替换。只有当第一个表达式为零时,才可能发生这种情况。表达式何时为零?
它有三个相乘的部分:

  • (i-j)%9,如果i和j相隔9的倍数,即在同一列,则其为零。
  • (i/9^j/9),如果i/9 == j/9,即相同的行,则其为零。
  • (i/27^j/27|i%9/3^j%9/3),如果这两个都为零,则为零:
    • i/27^j^27,如果i/27 == j/27,即三行的相同块,则为零
    • i%9/3^j%9/3,如果i%9/3 == j%9/3,即三列的相同块,则为零

如果这三个部分中的任何一个为零,则整个表达式为零。换句话说,如果i和j共享一行、一列或3x 3块,则j的值不能用作i处空格的候选项。啊哈!

from sys import exit, argv
def r(a):
  i = a.find('0')
  if i == -1:
    exit(a)
  for m in '3814697265625':
    okay = True
    for j in range(81):
      if (i-j)%9 == 0 or (i/9 == j/9) or (i/27 == j/27 and i%9/3 == j%9/3):
        if a[j] == m:
          okay = False
          break
    if okay:
      # At this point, m is not excluded by any row, column, or block, so let's place it and recurse
      r(a[:i]+m+a[i+1:])

r(argv[1])

注意,如果没有一个放置成功,r将返回到可以选择其他东西的点,所以这是一个基本的深度优先算法。
没有使用任何启发式,它不是特别有效。我从维基百科(http://en.wikipedia.org/wiki/Sudoku)这个难题:

$ time python sudoku.py 530070000600195000098000060800060003400803001700020006060000280000419005000080079
534678912672195348198342567859761423426853791713924856961537284287419635345286179

real    0m47.881s
user    0m47.223s
sys 0m0.137s

附录:作为一名维护程序员,我将如何重写它(这个版本的加速比大约为93倍:)

import sys

def same_row(i,j): return (i/9 == j/9)
def same_col(i,j): return (i-j) % 9 == 0
def same_block(i,j): return (i/27 == j/27 and i%9/3 == j%9/3)

def r(a):
  i = a.find('0')
  if i == -1:
    sys.exit(a)

  excluded_numbers = set()
  for j in range(81):
    if same_row(i,j) or same_col(i,j) or same_block(i,j):
      excluded_numbers.add(a[j])

  for m in '123456789':
    if m not in excluded_numbers:
      # At this point, m is not excluded by any row, column, or block, so let's place it and recurse
      r(a[:i]+m+a[i+1:])

if __name__ == '__main__':
  if len(sys.argv) == 2 and len(sys.argv[1]) == 81:
    r(sys.argv[1])
  else:
    print 'Usage: python sudoku.py puzzle'
    print '  where puzzle is an 81 character string representing the puzzle read left-to-right, top-to-bottom, and 0 is a blank'
6yoyoihd

6yoyoihd2#

使其清晰:

def r(a):
    i = a.find('0') # returns -1 on fail, index otherwise
    ~i or exit(a) # ~(-1) == 0, anthing else is not 0
                  # thus: if i == -1: exit(a)
    inner_lexp = [ (i-j)%9*(i/9 ^ j/9)*(i/27 ^ j/27 | i%9/3 ^ j%9/3) or a[j] 
                   for j in range(81)]  # r appears to be a string of 81 
                                        # characters with 0 for empty and 1-9 
                                        # otherwise
    [m in inner_lexp or r(a[:i]+m+a[i+1:]) for m in'%d'%5**18] # recurse
                            # trying all possible digits for that empty field
                            # if m is not in the inner lexp

from sys import *
r(argv[1]) # thus, a is some string

所以,我们只需要计算出内部列表表达式。我知道它收集行中设置的数字--否则,它周围的代码就没有意义了。然而,我不知道它是如何做到这一点的(对不起,我现在太累了,无法计算出二进制的幻想)

nbysray5

nbysray53#

r(a)是一个递归函数,它尝试在每一步中填充电路板中的0
i=a.find('0');~i or exit(a)是成功时终止。如果板中不再存在0值,则我们完成。
m是我们将尝试填充0的当前值。
如果将m放入当前0明显不正确,则m in[(i-j)%9*(i/9^j/9)*(i/27^j/27|i%9/3^j%9/3)or a[j]for j in range(81)]的值为truthy。我们将其命名为“is_bad”。这是最棘手的一点。:)
is_bad or r(a[:i]+m+a[i+1:]是一个条件递归步骤。如果当前候选解看起来是正常的,它将递归地尝试评估棋盘中的下一个0
for m in '%d'%5**18枚举从1到9的所有数字(效率低下)。

twh00eeo

twh00eeo4#

很多的短数独解算器只是递归地尝试每一个可能的法律的数字,直到他们成功地填满了单元格,我没有把它拆开,只是略读一下,看起来它就是这样做的。

7uhlpewt

7uhlpewt5#

代码实际上并不工作。你可以自己测试一下。下面是一个未解数独谜题的例子:
807000003602080000000200900040005001000798000200100070004003000000040108300000506
你可以使用这个网站(http://www.sudokuwiki.org/sudoku.htm),点击import puzzle,简单地复制上面的字符串,python程序的输出是:817311213622482322131224934443535441555798655266156777774663869988847188399979596
这与解不对应。事实上你已经可以看到一个矛盾,第一行有两个1。

相关问题