A*算法与八数码问题(numpy)

news2024/11/17 7:37:01

在这里插入图片描述

努力是为了不平庸~

学习的最大理由是想摆脱平庸,早一天就多一份人生的精彩;迟一天就多一天平庸的困扰。

目录

一、引言 

二、思路

1. 确定问题和目标:

2. 确定算法和数据结构:

3. 编写代码框架

4. 实现辅助函数:

5. 实现主函数:

6. 实现搜索算法:

7. 实现状态操作函数:

8. 进行测试和调试:

 三、代码与函数

A、代码步骤的方法、目的及意义

①导入必要的库:

②定义open表与close表以及初始状态和目标状态:

③定义节点类:

④定义辅助函数,实现节点类:

⑤定义函数get_reverse_num(state)用于计算逆序数,判断解的存在性。

⑥定义函数display(cur_node)用于输出状态及深度信息。

⑦定义函数is_in_list(alist, state)用于检查状态是否在列表中。

⑧定义排序的权值函数delta(node)。

⑨输入初始状态和目标状态:

⑩判断是否存在可行解:

⑪启发式搜索的实现:

B、函数构造和作用

1. 函数get_reverse_num(state):

2. 函数display(cur_node):

3. 函数is_in_list(alist, state):

4. 函数delta(node):

四、state和parent的属性和作用

1.state属性:

2. parent属性:

五、解题代码


一、引言 

    A*算法是一种启发式图搜索算法,其特点在于对估价函数的定义上。对于一般的启发式图搜索,总是选择估价函数f值最小的节点作为扩展节点。因此,f是根据需要找到一条最小代价路径的观点来估算节点的,所以,可考虑每个节点n的估价函数值为两个分量:从起始节点到节点n的实际代价以及从节点n到达目标节点的估价代价。

 以8数码问题为例实现A*算法的求解程序。

 估价函数f(n) = g(n) + h(n)

g(n)=d(n)——结点n在搜索树中的深度

h(n)可选择h1(n)——结点n中“不在位”的数码个数 或 h2(n) =p(n)= 将牌“不在位”的距离和

二、思路

1. 确定问题和目标:

   八数码问题,即通过移动数字将初始状态转变为目标状态。目标是编写能够自动求解八数码问题的程序。

2. 确定算法和数据结构:

   选择A*算法,然后确定所需的数据结构,如节点类、状态表示方法等。

3. 编写代码框架

   ①定义必要的变量和数据结构,如open表、close表、初始状态和目标状态等。

   ②创建节点类,并定义节点类的属性和方法。

4. 实现辅助函数:

   根据代码中的需要,实现辅助函数,如find_pos函数用于找到数字在状态中的位置。

5. 实现主函数:

   ①编写主函数,处理用户输入和输出。

   ②获取用户输入的初始状态和目标状态。

   ③检查可行解的存在性,判断是否可以通过移动数字从初始状态到达目标状态。

   ④如果可行解存在,开始启发式搜索过程。

6. 实现搜索算法:

   ①在主函数中编写启发式搜索的逻辑。

   ②创建open表和close表,并初始化。

   ③使用循环迭代进行搜索,直到找到可行解或无解。

   ④在每次迭代中,根据算法规则进行状态扩展和判断,并更新open表和close表。

7. 实现状态操作函数:

   根据算法需求,实现状态操作函数,如移动空白格、判断状态是否在表中等。

8. 进行测试和调试:

   ①编写测试用例,包括不同的初始状态和目标状态。

   ②运行代码,并观察输出结果是否符合预期。

   ③根据需要进行调试和修正代码错误。

 三、代码与函数

A、代码步骤的方法、目的及意义

①导入必要的库:

   这里使用了NumPy库,它提供了Python中对多维数组进行高效操作的功能。

②定义open表与close表以及初始状态和目标状态:

   open表用于存放待扩展的节点,close表用于存放已扩展的节点。start_state和target_state分别表示初始状态和目标状态,都是3x3的二维数组。

③定义节点类:

   这里定义了一个名为Node的类,表示搜索过程中的节点。该类具有一些属性和方法,用于存储节点状态、计算代价等信息。

④定义辅助函数,实现节点类:

   find_pos(self, state, num):在状态state中查找值为num的元素的位置。

   __init__(self, state, prt=[]):节点类的构造函数,初始化节点的状态、父节点、代价等信息。

   moveto(self, x, y):将空白格移动到指定位置(x, y),生成新的状态。

⑤定义函数get_reverse_num(state)用于计算逆序数,判断解的存在性。

   意义:逆序数用于判断八数码问题是否有可行解。

⑥定义函数display(cur_node)用于输出状态及深度信息。

   意义:用于展示搜索过程中每个节点的状态、搜索深度、代价等信息。

⑦定义函数is_in_list(alist, state)用于检查状态是否在列表中。

   意义:用于判断一个状态是否已经在open表或close表中。

⑧定义排序的权值函数delta(node)。

意义:用于指定节点的排序权值,这里使用节点的F值作为排序依据。

⑨输入初始状态和目标状态:

   意义:使用输入函数逐个获取用户输入的初始状态和目标状态。

⑩判断是否存在可行解:

    使用逆序数判断八数码问题是否有可行解。如果逆序数的奇偶性不同,则不存在可行解。

⑪启发式搜索的实现:

    在启发式搜索中,使用open表存放待扩展的节点,close表存放已扩展的节点。通过不断扩展节点,并将其加入open表和close表,最终找到可行解。搜索过程中,使用delta函数对open表中的节点按照F值进行排序

以上是对代码的关键步骤的分析。

B、函数构造和作用

1. 函数get_reverse_num(state):

   构成:

          ①参数:state,表示一个状态的二维数组。

          ②返回值:逆序数,表示状态中元素的逆序对数量。

   作用:该函数用于计算逆序数,判断八数码问题是否有可行解。通过将非空元素按顺序拼接成字符串s,然后统计在s中出现的逆序对数量,最终得到逆序数。

2. 函数display(cur_node):

   构成:

          ①参数:cur_node,表示当前节点对象。

          ②作用:该函数用于展示搜索过程中每个节点的状态、搜索深度、代价等信息。首先,它通过跟踪节点的父节点,将搜索路径中的所有节点存储在列表alist中。然后,它按照搜索路径的顺序遍历alist,打印每个节点的搜索深度G、状态state、代价信息G、H、F。

3. 函数is_in_list(alist, state):

   构成:

           ①参数:alist,表示一个节点列表;state,表示一个状态的二维数组。

           ②返回值:如果状态state在列表alist中存在,则返回对应的节点对象;否则返回-1。

   作用:该函数用于判断一个状态是否已经在open表或close表中。它遍历列表alist中的节点对象,将每个节点的状态与输入的state进行比较。如果找到匹配的状态,则返回对应的节点对象;否则返回-1,表示状态不在列表中。

4. 函数delta(node):

   构成:

            ①参数:node,表示一个节点对象。

            ②返回值:节点的代价F。

   作用:该函数用于作为排序的权值函数。在启发式搜索中,根据节点的代价F值对open表中的节点进行排序。

四、state和parent的属性和作用

state和parent是节点类中的两个重要属性

1.state属性:

   构成:state是一个二维数组,表示节点对应的状态。

   作用:state记录当前节点状态,即一个具体的排列。通过不同的状态组合,可以表示不同的八数码棋盘布局。

   例:在代码中的使用:通过current_node.state可以获取当前节点的状态。

2. parent属性:

   构成:parent是一个列表,存储指向父节点的引用。

   作用:parent记录了当前节点的父节点,即在搜索树中指向当前节点的上一层节点。

   例:在代码中的使用:通过current_node.parent可以获取当前节点的父节点。这对于回溯搜索路径和展示解决方案非常有用,因为可以通过跟踪父节点,从目标节点一直追溯到初始节点,形成搜索路径。

在启发式搜索中,每个节点都通过移动空白格生成新的状态,并将其作为子节点添加到搜索树中。通过state和parent属性的设置和更新,可以构建搜索树,跟踪搜索路径,并最终找到解决八数码问题的路径。

五、解题代码

import numpy as np

# 定义open表与close表
open_list = []  # 存放待扩展的节点
close_list = []  # 存放已扩展的节点
start_state = np.zeros((3, 3), dtype=int)  # 初始状态,3×3矩阵,先全设为零
target_state = np.zeros((3, 3), dtype=int)  # 目标状态,3×3矩阵,先全设为零

# 定义节点类
# 定义了一个名为Node的类,表示搜索过程中的节点。该类具有一些属性和方法,用于存储节点状态、计算代价等信息。
class Node:
    G = 0  # 起点到当前节点的最短路径代价值
    H = 0  # 当前节点到目标节点的最短路径代价值
    F = 0  # F = G + H,起点经过当前节点到目标节点的最短路径的总代价值
    state = np.zeros((3, 3), dtype=int)  # 一个二维数组,记录了当前节点的状态,即八数码问题中的一个具体排列。通过不同的状态组合,可以表示不同的八数码棋盘布局。
    parent = []  # 一个列表,记录了当前节点的父节点,即在搜索树中指向当前节点的上一层节点。


    #find_pos(self, state, num)是节点类中的一个方法,用于找到指定数字在状态中的位置
    # self:表示节点对象自身。  state:表示一个状态的二维数组。  num:表示要查找的数字。  返回一个包含两个元素的元组(i, j),表示数字num在状态state中的行索引和列索引。
    def find_pos(self, state, num):             # 在状态state中查找值为num的元素的位置x,y
        for i in range(len(state)):
            for j in range(len(state[i])):
                if state[i][j] == num:
                    return i, j

    # 节点类的构造函数,初始化节点的状态、父节点、代价等信息。
    def __init__(self, state, prt=[]):
        self.state = state
        if prt:
            self.parent = prt
            self.G = prt.G + 1
        for i in range(len(state)):
            for j in range(len(state[i])):
                x, y = self.find_pos(target_state, state[i][j])
                self.H = self.H + abs(x - i) + abs(y - j)
        self.F = self.G * 1 + self.H * 10

    # 将空白格移动到指定位置(x, y),生成新的状态。
    def moveto(self, x, y):
        x0, y0 = self.find_pos(self.state, 0)
        newstate = self.state.copy()
        tmp = newstate[x0][y0]
        newstate[x0][y0] = newstate[x][y]
        newstate[x][y] = tmp
        return newstate

# 得到逆序数,用于判断解的存在性
# 逆序数用于判断八数码问题是否有可行解
def get_reverse_num(state):
    ans = 0
    s = ""
    for i in range(len(state)):
        for j in range(len(state[i])):
            if not state[i][j] == 0:
                s += str(state[i][j])

    for i in range(len(s)):
        for j in range(i):
            if s[j] > s[i]:
                ans += 1
    return ans

# 输出状态及深度信息
# 用于展示搜索过程中每个节点的状态、搜索深度、代价等信息
def display(cur_node):
    alist = []
    tmp = cur_node
    while tmp:
        alist.append(tmp)
        tmp = tmp.parent
    alist.reverse()
    for node in alist:
        print("搜索深度%d" % node.G)
        print(node.state)
        print("G={0},H={1},F={2}".format(node.G, node.H, node.F))
        print()
    print("总共探索了%d次" % int(node.G + 1))

# 检查state状态是否在list中(可能是open或close表)
# 用于判断一个状态是否已经在open表或close表中
def is_in_list(alist, state):
    for stat in alist:
        if (stat.state == state).all():
            return stat
    return -1

# 排序的权值函数
# 用于指定节点的排序权值,这里使用节点的F值作为排序依据
def delta(node):
    return node.F

# 输入初始与目标状态
# 使用输入函数逐个获取用户输入的初始状态和目标状态
for i in range(len(start_state)):
    for j in range(len(start_state[i])):
        start_state[i][j] = input("start_state" + "(" + str(i + 1) + "," + str(j + 1) + "):")

print("the start state is:")
print(start_state)

for i in range(len(target_state)):
    for j in range(len(target_state)):
        target_state[i][j] = input("target_state" + "(" + str(i + 1) + "," + str(j + 1) + "):")

print("the target state is:")
print(target_state)

# 可行解判断
# 使用逆序数判断八数码问题是否有可行解。如果逆序数的奇偶性不同,则不存在可行解
if get_reverse_num(target_state) % 2 != get_reverse_num(start_state) % 2:
    print(get_reverse_num(target_state))
    print(get_reverse_num(start_state))
    print("找不到可行解!")
    exit(-1)

# 可行解存在时,开始启发搜索
# 在启发式搜索中,使用open表存放待扩展的节点,close表存放已扩展的节点。
# 通过不断扩展节点,并将其加入open表和close表,最终找到可行解。搜索过程中,使用delta函数对open表中的节点按照F值进行排序
open_list.append(Node(start_state))
while open_list:
    current_node = open_list.pop(0)  # 从open表中取出F值最小的节点进行扩展
    close_list.append(current_node)  # 将扩展的节点加入close表
    # 当open表中取出的恰好为目标状态时,找到可行解
    if (current_node.state == target_state).all():
        print("可行解已找到!")
        display(current_node)
        exit(0)
    # 否则对当前节点进行拓展
    x, y = current_node.find_pos(current_node.state, 0)  # 找到空白格的位置
    for [x_, y_] in [[x + 1, y], [x - 1, y], [x, y + 1], [x, y - 1]]:
        if 0 <= x_ < len(start_state) and 0 <= y_ < len(start_state):
            new_state = current_node.moveto(x_, y_)  # 移动空白格得到新状态
            # 判断新状态是否在close表
            if is_in_list(close_list, new_state) == -1:
                # 如果不在close表
                if is_in_list(open_list, new_state) == -1:
                    # 如果也不在open表,则将其加入open表
                    open_list.append(Node(new_state, current_node))
                else:
                    # 如果open表中已存在这种状态,则选取G值较小的
                    index = is_in_list(open_list, new_state)
                    if current_node.G + 1 < open_list[index].G:
                        # 如果新路线更好,则放弃旧路线而选择新路线
                        open_list.pop(index)
                        open_list.append(Node(new_state, current_node))
                    # 否则忽略
    # 对open表按F值从小到大进行排序
    open_list.sort(key=delta)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/642250.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

CreateML 使用以及机器学习基础概念

1. 前言 在学习 CreateML 之前&#xff0c;我们先了解一下什么是机器学习&#xff1f;目前还不存在被广泛认可的定义来准确定义机器学习是什么。第一个机器学习的定义来自于Arthur Samuel。他定义机器学习为&#xff0c;在进行特定编程的情况下&#xff0c;给予计算机学习能力…

C语言-类型转换

数据有不同的类型&#xff0c;不同类型数据之间进行混合运算时必然涉及到类型的转换问题. 转换的方法有两种&#xff1a; 自动转换: 遵循一定的规则,由编译系统自动完成. 强制类型转换&#xff1a; 把表达式的运算结果强制转换成所需的数据类型 1 自动转换 自动转换原则 …

阿里人手一本的Java性能调优手册,几乎涵盖了性能优化的所有操作

说起性能调优&#xff0c;想必大家都知道&#xff0c;但是就是没怎么用过&#xff0c;所以在Java性能优化上面临着很多的困扰&#xff0c;比如&#xff1a; 能力修炼中&#xff0c;由于常年接触 CRUD&#xff0c;缺乏高并发这一实践环境&#xff0c;对“性能优化”只能通过理论…

剑指 Offer 59 - I: 滑动窗口的最大值

第一眼看这个感觉思路没啥大问题&#xff0c;就是一个大循环不断后移&#xff0c;然后小循环维护一个k长度的队列。注意peekFirst和peekLast的使用&#xff0c;双端队列可以打破先进先出或者先进后出的局限性&#xff0c;因此peek没有默认值&#xff0c;得加上First或者Last来进…

干货 | 腾讯云李滨:云时代数据安全治理方法与实践

作者&#xff1a;李滨‍‍‍‍本文约2600字&#xff0c;建议阅读8分钟 本文与你分享腾旭在数据安全和人工智能的安全治理方法&#xff0c;以及实践上的一些经验。 今天很荣幸与大家分享在过去几年以内&#xff0c;腾讯在数据安全和人工智能的安全治理方法&#xff0c;以及实践上…

使用服务器的妙妙工具

为什么使用服务器&#xff1f; 一台电脑的性能终究是有限的&#xff0c;所以当计算量太大时&#xff0c;我们可以借助服务器&#xff0c;它的运算速度快&#xff0c;存储空间大。 和个人电脑不同&#xff0c;服务器通常由一个管理员进行管理&#xff0c;而有多用户同时使用。…

一行代码实现网页直接跳转抖音,2023年6月新方法

之前 发的文章 方法已经失效&#xff0c;以下为新版本方案&#xff01; 背景 大家可能都发现在一些APP中无法直接打开抖音的链接&#xff0c;就连抖音官方的分享链接在浏览器打开都无法直接唤醒APP。需要我们手动点击&#xff0c;如果我们在自己网站等推广希望直接拉起抖音可…

工厂方法模式(五)

过气的&#xff0c;终究是过气了 上一章简单介绍了工厂模式(四), 如果没有看过,请观看上一章 一.工厂方法模式 工厂方法模式&#xff0c;通过定义工厂父类负责定义创建对象的公共接口&#xff0c;而子类则负责生成具体的对象。 将类的实例化&#xff08;具体产品的创建&…

全网最火爆,性能测试基本原则方法总结,再不卷就被淘汰了...

目录&#xff1a;导读 前言一、Python编程入门到精通二、接口自动化项目实战三、Web自动化项目实战四、App自动化项目实战五、一线大厂简历六、测试开发DevOps体系七、常用自动化测试工具八、JMeter性能测试九、总结&#xff08;尾部小惊喜&#xff09; 前言 性能测试关键点 …

2023 年牛客网最新版大厂 Java 八股文面试题总结(覆盖所有面试题考点)

序员真的是需要将终生学习贯彻到底的职业&#xff0c;一旦停止学习&#xff0c;离被淘汰&#xff0c;也就不远了。 金九银十跳槽季&#xff0c;这是一个千年不变的话题&#xff0c;每到这个时候&#xff0c;很多人都会临阵磨枪&#xff0c;相信不快也光。于是&#xff0c;大家…

什么是Vue的数据流(单向数据流)?如何进行数据流管理

什么是Vue的数据流&#xff08;单向数据流&#xff09;&#xff1f;如何进行数据流管理 在Vue中&#xff0c;数据流是指数据的传递和管理方式。Vue采用的是单向数据流&#xff0c;也就是说&#xff0c;数据是从父组件流向子组件&#xff0c;子组件不能直接修改父组件的数据。本…

夏季来临居民如何预防电气火灾隐患

安科瑞虞佳豪 居住场所预防电气火患 请注意以下几点 1.用具有生产许可证或CCC证书并与居住场所的环境相适应的电气线路以及插座、插线板、充电器、用电设备等电器产品&#xff1b; 2.不私自拉接电气线路&#xff1b;不随意更改电气线路保护装置&#xff1b; 3.不随意更换原…

泛娱乐社交(一)直播产品商业化解决方案

摘要 在过去几年的直播行业创业风口期中&#xff0c;直播的用户关注度疯狂增长&#xff0c;但用户质量却参差不齐。随着用户新鲜感一过&#xff0c;流失率变得相当严重&#xff0c;各大平台都在竭尽全力防御。然而&#xff0c;留住“凑热闹”的非直播受众用户并不是最关键的问…

SpringBoot(基础篇)

SpringBoot基础篇 入门案例 在创建SpringBoot项目时&#xff0c;会出现以下不需要的文件&#xff0c;如果每次都手动删除的话&#xff0c;就会很麻烦。 教你一招 在setting设置中找到Editor&#xff0c;选择File Types–>Ignored Files and Folders–>点击号&#xff…

pycharm使用之torch_geometric安装

正式安装之前要先查看一下torch的版本 一、查看torch版本 1、winR &#xff0c;输入cmd 2、输入python 3、 输入import torch&#xff0c;然后输入torch.__version__&#xff0c;最后回车 可以看到我的torch版本是1.10.0 二、下载合适的版本 1、打开链接 https://pytorch-…

千万级入口服务[Gateway]框架设计(一)

本文将以技术调研模式编写&#xff0c;非技术同学可跳过。 文章目录 背景问题[不涉及具体业务]目标技术选型语言框架模式实现一&#xff1a;go 原生组件Demo 实现Benchwork 基准性能小结实现二&#xff1a;开源 go-plugin 附录入口服务演变 背景 在历史架构的迭代中&#xff…

Elasticsearch:部署 ELSER - Elastic Learned Sparse EncoderR

警告&#xff1a;此功能处于技术预览阶段&#xff0c;可能会在未来的版本中更改或删除。 Elastic 将尽最大努力修复任何问题&#xff0c;但技术预览中的功能不受官方 GA 功能的支持 SLA 约束。 Elastic Learned Sparse EncodeR - 或 ELSER - 是由 Elastic 训练的检索模型&#…

Python程序员必会技能:如何在Python中连接数据库获取数据

前言 嗨喽&#xff0c;大家好呀~这里是爱看美女的茜茜呐 想要成为一名数据分析师&#xff0c;一个最基础的技能就是&#xff1a;提取数据。 数据是数据分析的基础。 因此&#xff0c;掌握各种获取数据的技能一定是一个合格的数据分析师的必备技能。 而我们的需要进行提取的…

总结排查服务器上传下载慢的几种手段与查看服务器带宽的具体方法

一、排查服务器上传下载 最近出现的一个情况&#xff0c;服务器上传和下载比较慢&#xff0c;因此我排查了种种手段&#xff0c;特此记录下几种常见的手段。 1、使用speedtest-cli 测试网速&#xff1a; 该方法是测试网速的速度怎么样&#xff0c;看看是否真的慢&#xff1f; …

python中os库用法详解(总结)

os库主要是对文件和文件夹进行操作&#xff0c;在Python中对⽂件和⽂件夹的操作要借助os模块⾥⾯的相关功能。 具体步骤如下&#xff1a; 1. 导⼊os模块 import os 2. 使⽤ os 模块相关功能 os.函数名() 1、⽂件重命名 os.rename(⽬标⽂件名, 新⽂件名) 示例代码&#…