Python数独求解器:从基础到回溯算法详解

Python数独求解器:从基础到回溯算法详解

本教程详细介绍了如何使用Python构建一个数独求解器。文章首先分析了数独求解中的常见问题,特别是文件操作和回溯逻辑的误区。随后,提供了两种核心解决方案:一种是基于回溯算法的通用数独求解器,能够解决任何有效数独;另一种是迭代式“单解”填充器,适用于仅需填充唯一确定单元格的简单数独。教程涵盖了代码实现、原理分析及关键注意事项,旨在帮助读者深入理解数独求解的算法思想。

1. 数独求解器基础:网格表示与有效性检查

数独求解的核心在于两个方面:一是如何表示数独网格,二是如何判断一个数字在特定位置是否有效。

1.1 网格表示

一个标准的9×9数独可以被表示为一个二维列表(或称列表的列表),其中0表示空单元格。

# 示例数独网格grid = [    [0, 4, 0, 0, 0, 0, 1, 7, 9],    [0, 0, 2, 0, 0, 8, 0, 5, 4],    [0, 0, 6, 0, 0, 5, 0, 0, 8],    [0, 8, 0, 0, 7, 0, 9, 1, 0],    [0, 5, 0, 0, 9, 0, 0, 3, 0],    [0, 1, 9, 0, 6, 0, 0, 4, 0],    [3, 0, 0, 4, 0, 0, 7, 0, 0],    [5, 7, 0, 1, 0, 0, 2, 0, 0],    [9, 2, 8, 0, 0, 0, 0, 6, 0]]

main 函数负责从文件读取输入并将其转换为这种网格格式:

import sysdef main():    # 从命令行参数获取输入文件路径    with open(sys.argv[1], 'r') as f:        s1 = f.read()        s2 = s1.split()        # 将字符串转换为整数列表        s_int_list = [int(x) for x in s2]        # 将一维列表转换为9x9的二维网格        grid = [s_int_list[i:i+9] for i in range(0, len(s_int_list), 9)]        # 接下来调用求解函数        # solve(grid) # 具体调用取决于采用哪种求解策略

1.2 有效性检查 (check 函数)

check 函数是数独求解的基石,它判断在给定行 r、列 c 的位置放置数字 k 是否符合数独规则。规则包括:

同行不能有重复数字。同列不能有重复数字。同一3×3宫格内不能有重复数字。

def check(grid, r, c, k):    # 检查行    for i in range(9):        if grid[r][i] == k:            return False    # 检查列    for i in range(9):        if grid[i][c] == k:            return False    # 检查3x3宫格    # 计算当前单元格所属3x3宫格的左上角坐标    x_area = (c // 3) * 3    y_area = (r // 3) * 3    for i in range(3):        for j in range(3):            if grid[y_area + i][x_area + j] == k:                return False    return True

2. 问题分析:原始代码的挑战与局限

原始代码尝试实现一个递归求解器,但存在几个关键问题,导致它只能输出第一步:

立即学习“Python免费学习笔记(深入)”;

文件重复打开与覆盖: 在每次递归调用 solve 函数时,都会执行 f = open(sys.argv[2],’w’)。’w’ 模式会清空文件内容,导致每次写入都覆盖之前的内容,最终文件中只保留最后一次写入的数据。文件未关闭: 文件句柄 f 在函数内部打开,但没有显式关闭。虽然 with open(…) 结构能确保文件关闭,但原始代码的 f = open(…) 没有使用 with 语句,可能导致资源泄露或数据未完全写入。poss 列表逻辑误解: 原始代码试图通过 len(poss) == 1 来判断是否只有一个可能值,并立即填充。这并没有真正实现“只填充唯一可能值”的逻辑,因为 poss 列表在循环中首次添加元素后,其长度就变为1。更重要的是,这种策略不足以解决所有数独,因为它没有尝试所有可能性,也没有回溯机制。缺乏回溯机制: 当某个尝试的数字无法导致最终解时,原始代码没有将该单元格重置为0(即回溯),而是直接返回 False,导致无法探索其他可能的路径。对于复杂的数独,必须要有回溯才能找到解。

3. 解决方案一:基于回溯法的通用数独求解器

回溯法是一种穷举搜索算法,适用于解决约束满足问题。其核心思想是:尝试一个可能的值,如果发现这条路径无法通向解,就“回溯”到上一步,撤销之前的尝试,然后尝试另一个可能的值。

3.1 原理阐述

选择空单元格: 从网格中找到一个空的单元格(值为0)。尝试数字: 遍历数字1到9,对于每个数字 k,使用 check 函数判断其是否可以在当前空单元格放置。递归填充: 如果 k 有效,则将其放入单元格,并递归调用求解函数来填充下一个空单元格。成功与回溯:如果递归调用返回 True(表示后续单元格都成功填充),则当前路径有效,返回 True。如果递归调用返回 False(表示当前 k 无法导致解),则将当前单元格重置为0(回溯),然后尝试下一个数字 k。无解: 如果所有数字1到9都尝试完毕,仍无法找到有效解,则说明当前路径无解,返回 False。终止条件: 如果所有单元格都被填充(即没有空单元格),则表示数独已解决,返回 True。

3.2 代码实现

为了解决文件重复打开和回溯问题,我们将文件操作和计数器变量提升到外部作用域,并使用一个嵌套函数 recur 来实现递归回溯逻辑。

import sysdef main():    with open(sys.argv[1], 'r') as f:        s1 = f.read()        s2 = s1.split()        s_int_list = [int(x) for x in s2]        grid = [s_int_list[i:i+9] for i in range(0, len(s_int_list), 9)]        solve(grid) # 调用新的 solve 函数def check(grid, r, c, k):    # 检查行    for i in range(9):        if grid[r][i] == k:            return False    # 检查列    for i in range(9):        if grid[i][c] == k:            return False    # 检查3x3宫格    x_area = (c // 3) * 3    y_area = (r // 3) * 3    for i in range(3):        for j in range(3):            if grid[y_area + i][x_area + j] == k:                return False    return Truedef solve(grid):    # 文件只打开一次,使用 'w' 模式会清空文件,但因为只打开一次,所以不会重复清空    # 更好的做法是使用 'a' 模式(追加)或者在 solve 外部处理文件    # 这里为了与原始需求保持一致,假设每次运行都生成新的输出文件    with open(sys.argv[2], 'w') as f: # 使用 with 语句确保文件关闭        counter = 0 # 计数器只初始化一次        # 嵌套函数实现递归,可以访问外部 solve 函数的 counter 和 f        def recur(r, c):            nonlocal counter # 声明 counter 为非局部变量,以便修改外部函数的 counter            # 基础情况:如果行索引达到9,表示所有行都已处理完毕,数独已解决            if r == 9:                return True            # 如果列索引达到9,表示当前行已处理完毕,移动到下一行的开头            elif c == 9:                return recur(r + 1, 0)            # 如果当前单元格不为0(已被填充),则跳过,处理下一个单元格            elif grid[r][c] != 0:                return recur(r, c + 1)            else:                # 尝试填充1到9的数字                for k in range(1, 10):                    if check(grid, r, c, k): # 检查当前数字 k 是否有效                        grid[r][c] = k # 放置数字                        counter += 1                        # 打印当前步骤到文件                        print("-" * 18,                               "Step " + str(counter) + " - " + str(k)                                       + " @ " + "R" + str(r + 1) + "C" + str(c + 1),                               "-" * 18,                              sep='n', file=f)                        for x in grid:                            print(" ".join(map(str, x)), file=f)                        print("-" * 18, file=f)                        # 递归调用以填充下一个单元格                        if recur(r, c + 1):                            return True # 如果后续填充成功,则当前路径有效                # 回溯:如果所有数字都尝试过,但没有一个能导致解,                # 则将当前单元格重置为0,并返回 False                grid[r][c] = 0                 return False        # 从 (0, 0) 开始递归求解        result = recur(0, 0)    return result # 返回最终求解结果(True/False)if __name__ == "__main__":    main()

关键改进点:

文件句柄管理: f = open(sys.argv[2],’w’) 被移动到 solve 函数的开头,并使用 with 语句,确保文件只打开一次且在函数结束时自动关闭。计数器 counter: counter 也被移到 solve 函数内部,并使用 nonlocal 关键字在嵌套函数 recur 中进行修改,确保其累加正确。回溯逻辑: 在 for 循环结束后,如果没有任何数字能成功填充当前单元格并导致解,grid[r][c] = 0 会将单元格重置,实现回溯。poss 列表移除: 回溯法不需要预先计算所有可能值,而是直接尝试1-9的每个数字。

4. 解决方案二:迭代式“单解”数独填充器

原始需求中提到“如果有一个单元格只有一个可能的数字,就填入这个数字并打印表格,然后重复这些步骤”。这是一种简化版的数独填充策略,不涉及回溯,只适用于那些可以通过逻辑推理逐步填充的“简单”数独。如果数独需要猜测并回溯,这种方法将无法解决。

4.1 原理阐述

循环迭代: 不断循环,直到无法再填充任何单元格或数独被解决。查找唯一解单元格: 在每次迭代中,遍历所有空单元格。对于每个空单元格,尝试数字1到9,找出所有可能的有效数字。填充唯一解: 如果某个空单元格只有一个可能的有效数字,则填充该数字,并打印当前数独状态。终止条件: 如果一轮遍历下来,没有找到任何具有唯一可能解的单元格,则停止。此时,数独可能已解决,或者是一个“复杂”数独,无法通过这种简单方法解决。

4.2 代码实现

import sys# check 函数与 main 函数保持不变,此处省略重复代码def solve_simple_sudoku(grid):    def count_empty_cells():        count = 0        for r in range(9):            for c in range(9):                if grid[r][c] == 0:                    count +=1        return count    def find_cell_with_one_solution():        """        在网格中寻找一个空单元格,该单元格只有一个可能的有效填充数字。        如果找到,返回 (行, 列, 唯一数字);否则返回 None。        """        for r in range(9):            for c in range(9):                if grid[r][c] == 0: # 找到空单元格                    poss = [] # 存储当前单元格的所有可能数字                    for k in range(1, 10):                        if check(grid, r, c, k):                            poss.append(k)                    if len(poss) == 1: # 如果只有一个可能数字                        return r, c, poss[0]        return None # 没有找到具有唯一解的空单元格    with open(sys.argv[2], 'w') as f:        # 循环次数最多为所有空单元格的数量,因为每次成功填充一个        # 实际可能更少,因为如果无法找到唯一解,会提前退出        for counter in range(count_empty_cells()):             result = find_cell_with_one_solution()            if not result:  # 如果找不到具有唯一解的空单元格                # 如果所有单元格都已填充,则数独已解决                if count_empty_cells() == 0:                    print("Sudoku solved successfully!", file=f)                else: # 否则,这个数独对这种方法来说太复杂了                    print("This is not a simple Sudoku puzzle! Requires backtracking.", file=f)                return False # 无法通过此方法解决            r, c, k = result # 获取找到的唯一解单元格信息            grid[r][c] = k   # 填充数字            # 打印当前步骤和网格状态            print("-" * 18,                   "Step " + str(counter+1) + " - " + str(k)                           + " @ " + "R" + str(r + 1) + "C" + str(c + 1),                   "-" * 18,                  sep='n', file=f)            for x in grid:                print(" ".join(map(str, x)), file=f)            print("-" * 18, file=f)        # 如果循环结束且所有单元格都已填充        if count_empty_cells() == 0:            print("Sudoku solved successfully!", file=f)            return True        else:            print("Finished simple filling, but puzzle not fully solved.", file=f)            return False# 在 main 函数中调用此 solve_simple_sudoku# if __name__ == "__main__":#     # 假设 main 函数已经处理了文件读取和 grid 初始化#     # 然后可以这样调用:#     # main()#     # 并在 main 函数内部调用 solve_simple_sudoku(grid)

注意事项:

此方法不适用于所有数独。如果一个数独需要通过试错和回溯才能解决,find_cell_with_one_solution 将会返回 None,导致程序提前终止,并指出其“不是一个简单数独”。counter 在这里表示成功填充的步数,而不是递归调用的次数。

5. 总结与最佳实践

本教程展示了两种不同策略的数独求解器:

基于回溯的通用求解器: 适用于解决任何有效数独难题。它通过递归尝试所有可能性并在遇到死胡同时回溯,是解决这类问题的标准算法。文件操作和计数器变量的正确管理是实现的关键。迭代式“单解”填充器: 满足特定需求,即只填充那些具有唯一确定解的单元格。这种方法更直观,但其局限性在于无法解决需要猜测和回溯的复杂数独。

在实际开发中,应根据数独的复杂度和预期行为选择合适的求解策略。对于通用的数独求解需求,回溯法是首选。同时,良好的文件I/O习惯(如使用 with open 语句)和清晰的变量作用域管理(如 nonlocal 关键字)是编写健壮Python代码的重要实践。

以上就是Python数独求解器:从基础到回溯算法详解的详细内容,更多请关注创想鸟其它相关文章!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 chuangxiangniao@163.com 举报,一经查实,本站将立刻删除。
发布者:程序猿,转转请注明出处:https://www.chuangxiangniao.com/p/1366832.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
Python数独求解器:从基础回溯到单解填充策略
上一篇 2025年12月14日 06:54:44
Python 数独求解教程:详解回溯与迭代策略
下一篇 2025年12月14日 06:54:58

相关推荐

  • Matplotlib 地图中多类型图例的创建与优化

    Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化

    本教程旨在解决matplotlib地图可视化中,如何在一个图例中同时展示颜色块(如区域分类)和自定义标记(如特定兴趣点)的问题。文章详细介绍了当传统`patch`对象无法正确显示标记时,如何利用`matplotlib.lines.line2d`创建标记图例句柄,并将其与颜色块图例句柄合并,从而生成一…

    2026年5月10日 用户投稿
    100
  • Golang JSON序列化:控制敏感字段暴露的最佳实践

    本教程探讨golang中如何高效控制结构体字段在json序列化时的可见性。当需要将包含敏感信息的结构体数组转换为json响应时,通过利用`encoding/json`包提供的结构体标签,特别是`json:”-“`,可以轻松实现对特定字段的忽略,从而避免敏感数据泄露,确保api…

    2026年5月10日
    000
  • 利用海象运算符简化条件赋值:Python教程与最佳实践

    本文旨在探讨Python中海象运算符(:=)在条件赋值场景下的应用。通过对比传统if/else语句与海象运算符,以及条件表达式,分析海象运算符在简化代码、提高可读性方面的优势与局限性。并通过具体示例,展示如何在列表推导式等场景下合理使用海象运算符,同时强调其潜在的复杂性及替代方案,帮助开发者更好地掌…

    2026年5月10日
    100
  • 比特币新手教程 比特币交易平台有哪些

    比特币是一种去中心化的数字货币,基于区块链技术实现点对点交易,具有匿名性、有限发行和不可篡改等特点;新手可通过交易所购买,P2P交易获得比特币,常用平台包括Binance、OKX和Huobi;交易流程包括注册账户、实名认证、绑定支付方式、充值法币并下单购买,可选择市价单或限价单;比特币存储方式有交易…

    2026年5月10日
    000
  • c++中的SFINAE技术是什么_c++模板编程中的SFINAE原理与应用

    SFINAE 是“替换失败不是错误”的原则,指模板实例化时若参数替换导致错误,只要存在其他合法候选,编译器不报错而是继续重载决议。它用于条件启用模板、类型检测等场景,如通过 decltype 或 enable_if 控制函数重载,实现类型特征判断。尽管 C++20 引入 Concepts 简化了部分…

    2026年5月10日
    000
  • Go语言mgo查询构建:深入理解bson.M与日期范围查询的正确实践

    本文旨在解决go语言mgo库中构建复杂查询时,特别是涉及嵌套`bson.m`和日期范围筛选的常见错误。我们将深入剖析`bson.m`的类型特性,解释为何直接索引`interface{}`会导致“invalid operation”错误,并提供一种推荐的、结构清晰的代码重构方案,以确保查询条件能够正确…

    2026年5月10日
    100
  • RichHandler与Rich Progress集成:解决显示冲突的教程

    在使用rich库的`richhandler`进行日志输出并同时使用`progress`组件时,可能会遇到显示错乱或溢出问题。这通常是由于为`richhandler`和`progress`分别创建了独立的`console`实例导致的。解决方案是确保日志处理器和进度条组件共享同一个`console`实例…

    2026年5月10日
    000
  • Golang goroutine与channel调试技巧

    使用go run -race检测数据竞争,结合runtime.NumGoroutine监控协程数量,通过pprof分析阻塞调用栈,利用select超时避免永久阻塞,有效排查goroutine泄漏、死锁和数据竞争问题。 Go语言的goroutine和channel是并发编程的核心,但它们也带来了调试上…

    2026年5月10日
    000
  • 使用 Jupyter Notebook 进行探索性数据分析

    Jupyter Notebook通过单元格实现代码与Markdown结合,支持数据导入(pandas)、清洗(fillna)、探索(matplotlib/seaborn可视化)、统计分析(describe/corr)和特征工程,便于记录与分享分析过程。 Jupyter Notebook 是进行探索性…

    2026年5月10日
    000
  • 《魔兽世界》将于6月11日开启国服回归技术测试

    《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试《魔兽世界》将于6月11日开启国服回归技术测试

    《%ign%ignore_a_1%re_a_1%》官方宣布,将于6月11日开启国服回归技术测试,时间为7天,并称可以在6月内正式开服,玩家们可以访问官网下载战网客户端并预下载“巫妖王之怒”客户端,技术测试详情见下图。 WordAi WordAI是一个AI驱动的内容重写平台 53 查看详情 以上就是《…

    2026年5月10日 用户投稿
    200
  • php常量怎么用_PHP常量(define/const)定义与使用方法

    PHP中可通过define函数和const关键字定义常量,用于存储不可变值。define适用于全局作用域,支持动态名称和条件定义,如define(‘SITE_NAME’, ‘MyWebsite’);const在编译时生效,语法简洁但限制多,只能在类或全…

    2026年5月10日
    000
  • 如何在HTML中插入表单元素_HTML表单控件与输入类型使用指南

    HTML表单通过标签构建,包含action和method属性定义数据提交目标与方式,常用input类型如text、password、email等适配不同输入需求,配合label、required、placeholder提升可用性,结合textarea、select、button等控件实现完整交互,是…

    2026年5月10日
    100
  • 创建指定大小并填充特定数据的Golang文件教程

    本文将介绍如何使用Golang创建一个指定大小的文件,并用特定数据填充它。我们将使用 `os` 包提供的函数来创建和截断文件,从而实现快速生成大文件的目的。示例代码展示了如何创建一个10MB的文件,并将其填充为全零数据。掌握这些方法,可以方便地在例如日志系统或磁盘队列等场景中,预先创建测试文件或初始…

    2026年5月10日
    000
  • Python命令怎样使用profile分析脚本性能 Python命令性能分析的基础教程

    使用Python的cProfile模块分析脚本性能最直接的方式是通过命令行执行python -m cProfile your_script.py,它会输出每个函数的调用次数、总耗时、累积耗时等关键指标,帮助定位性能瓶颈;为进一步分析,可将结果保存为文件python -m cProfile -o ou…

    2026年5月10日
    000
  • 使用 WebCodecs VideoDecoder 实现精确逐帧回退

    本文档旨在解决在使用 WebCodecs VideoDecoder 进行视频解码时,实现精确逐帧回退的问题。通过比较帧的时间戳与目标帧的时间戳,可以避免渲染中间帧,从而提高用户体验。本文将提供详细的解决方案和示例代码,帮助开发者实现精确的视频帧控制。 在使用 WebCodecs VideoDecod…

    2026年5月10日
    000
  • 如何插入查询结果数据_SQL插入Select查询结果方法

    如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法如何插入查询结果数据_SQL插入Select查询结果方法

    使用INSERT INTO…SELECT语句可高效插入数据,通过NOT EXISTS、LEFT JOIN、MERGE语句或唯一约束避免重复;表结构不一致时可通过别名、类型转换、默认值或计算字段处理;结合存储过程可提升可维护性,支持参数化与动态SQL。 将查询结果数据插入到另一个表中,可以…

    2026年5月10日 用户投稿
    000
  • Discord.py 交互按钮超时与持久化解决方案

    本教程旨在解决Discord.py中交互按钮在一段时间后出现“This Interaction Failed”错误的问题。我们将深入探讨视图(View)的超时机制,并提供通过正确设置timeout参数以及利用bot.add_view()方法实现按钮持久化的具体方案,确保您的机器人交互功能稳定可靠,即…

    2026年5月10日
    000
  • Debian Copilot的社区活跃度如何

    debian copilot是codeberg社区维护的ai助手,旨在为debian用户提供服务。尽管搜索结果中没有直接提供关于debian copilot社区支持活跃度的具体数据,但我们可以通过debian社区的整体活跃度和特点来推断其活跃性。 Debian社区的一般情况: Debian拥有详尽的…

    2026年5月10日
    000
  • Python递归函数追踪与性能考量:以序列打印为例

    本文深入探讨了Python中一种递归打印序列元素的方法,并着重演示了如何通过引入缩进参数来有效追踪递归函数的执行流程和参数变化。通过实际代码示例,文章揭示了递归调用可能带来的潜在性能开销,特别是对调用栈空间的需求,以及Python默认递归深度限制可能导致的错误,为读者提供了理解和优化递归算法的实用见…

    2026年5月10日
    000
  • python中zip函数详解 python多序列压缩zip函数应用场景

    zip函数的应用场景包括:1) 同时遍历多个序列,2) 合并多个列表的数据,3) 数据分析和科学计算中的元素运算,4) 处理csv文件,5) 性能优化。zip函数是一个强大的工具,能够简化代码并提高处理多个序列时的效率。 在Python中,zip函数是一个非常有用的工具,它能够将多个可迭代对象打包成…

    2026年5月10日
    000

发表回复

登录后才能评论
关注微信