A*算法求解八数码问题的python实现

实验四 A算法求解八数码问题实验

一、实验目的

熟悉和掌握启发式搜索的定义、估价函数和算法过程,并利用 A*算法求解 N 数码难题,理解求解流程和搜索顺序。

二、实验原理

八数码问题是在3×3的九宫格棋盘上,摆有8个刻有1~8数码的将牌。棋盘中有一个空格,允许紧邻空格的某一将牌可以移到空格中,这样通过平移将牌可以将某一将牌布局变换为另一布局。针对给定的一种初始布局或结构(目标状态),问如何移动将牌,实现从初始状态到目标状态的转变。

三、实验条件

Python版本:Python 3及以上
所需要的依赖包:无
可采用的软件:PyCharm
知识库是教材中的动物识别系统中的15条规则。

四、实验内容

以8 数码问题和 15 数码问题为例实现 A算法的求解程序(编程语言不限,要求设计两种不同的估价函数。

五、实验步骤

具体步骤:
1.定义四个函数:move、create、show 和 hx,及一段主程序代码。
2.move 函数定义了移动操作,根据给定的方向和当前位置,计算新的位置。
3.create 函数用于生成新的节点,并记录路径。它创建一个新节点 p,将其插入到路径数组,然后检查新节点是否已经在 close 列表中,如果不存在则添加到 open 列表中。
4.show 函数打印八数码拼图的状态。
5.hx 函数是估价函数,用于评估当前状态与目标状态之间的差距。计算当前状态中与目标状态相同位置的数码个数。
6.主程序在 if name==“main”: 下开始执行。
7.初始化初始状态 start 和目标状态 goal,并创建空的 open 和 close 列表。
8.进入一个主循环,用于执行A*搜索。
9.首先检查初始状态是否已经是目标状态,如果是,直接输出初始状态即为解,并结束搜索。如果 open 列表为空,且仍未找到解,输出未找到解并结束搜索。否则,从 open 列表中取出第一个节点 this,并将其添加到 close 列表中表示已经探索。如果当前节点是目标节点,输出搜索成功,并打印出解路径。然后结束搜索。程序根据当前节点的位置,生成可能的下一步节点,包括向上、向下、向左和向右移动操作。对于每个生成的新节点,检查是否已经在 close 列表中,如果不在则将新节点添加到 open 列表中。
10.最后,根据估价函数 hx 来对 open 列表中的节点按估价值从高到低排序

六、实验结果及分析


图1:初始八数码与目标八数码列表


图2:实验结果

七、实验心得体会

通过本次实验,发现选用不同的启发函数,对于实验的结果有较大的影响。正如表3-1所示,选用第一或第二种(也就是采用A*算法)远远优于普通的广度优先搜索,同时,明显的感觉到第二种启发函数效率更高,更快的找到最优解。但是,在实验过程中,也遇到了一些问题,比如初始值的八数码初始值的选择对于实验结果的影响很大,在选取一些样例时,比如1,3,0,2,8,4,7,6,5,实验结果达到20000次依然没有停止,无法比较两种启发函数的优越性,鉴于时间原因,选取一些迭代次数较小就可以达到目标状态的样例进行验证,发现第二种结果优于第一种启发函数的结果。
import random
#修改索引,方便接下来交换
def move(str,index):
if str==‘up’:
return index-3
if str==‘dowm’:
return index+3
if str==‘left’:
return index-1
if str==‘right’:
return index+1

#生成新的结点并且记录路径
def create(array1,array2):
p=array1[:]
p.insert(0,array2)
for i in close:
if i[0] == p[0]:
return False
open.append§
return True
#打印
def show(list):
for j in range(len(list)):
if j % 3 == 0:
print(‘\t’)
print(list[j], end=’ ‘)
print(’\t’)

#估计函数,与目标结点对比,返回当前数码与目标数码相同位置的个数
def hx(current):
count=0
goal = [1, 2, 3, 8, 0, 4, 7, 6, 5]
for i in range(len(current)):
if current[i]==goal[i]:
count+=1
return count

if name==“main”:

# print(start)
start=[2,8,3,1,0,4,7,6,5]
goal=[1,2,3,8,0,4,7,6,5]

open=[]
close=[]
step=0
open.append([start])
while 1:
    if start==goal:
        print('初始状态即为解!')
        break
    if len(open)==0:
        print('未找到解')
        break
    else:
        this=open.pop(0)
      
        close.append(this)
        if this[0]==goal:
            print('搜索成功')
            print('共{}步达成目标'.format(len(this)-1))
            for i in this[::-1]:
                show(i)
            exit()

        #上
        if this[0].index(0)>2 :
            node=this[0].copy()
            a=this[0].index(0)
            b=move('up',a)
            node[a],node[b]=node[b],node[a]
            create(this, node)
           

        #下
        if this[0].index(0) < 6 :
            node = this[0].copy()
            a = this[0].index(0)
            b = move('dowm', a)
            node[a],node[b]=node[b],node[a]
            create(this, node)
           

        #左
        if this[0].index(0) != 0 and this[0].index(0) != 3 and this[0].index(0) != 6:
            node = this[0].copy()
            a = this[0].index(0)
            b = move('left', a)
            node[a],node[b]=node[b],node[a]
            create(this, node)
         

        #右
        if this[0].index(0) != 2 and this[0].index(0) != 5 and this[0].index(0) != 8:
            node = this[0].copy()
            a = this[0].index(0)
            b = move('right', a)
            node[a],node[b]=node[b],node[a]
            create(this, node)
           
        #对与目标数码结点相同个数最多的数码结点放到open表最前面
        for i in range(len(open)-1):
            for j in range(i+1,len(open)):
                if hx(open[i][0])<hx(open[j][0]):
                    open[i],open[j]=open[j],open[i]
  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
A*算法求解八数码问题 1、A*算法基本思想: 1)建立一个队列,计算初始结点的估价函数f,并将初始结点入队,设置队列头和尾指针。 2)取出队列头(队列头指针所指)的结点,如果该结点是目标结点,则输出路径,程序结束。否则对结点进行扩展。 3)检查扩展出的新结点是否与队列中的结点重复,若与不能再扩展的结点重复(位于队列头指针之前),则将它抛弃;若新结点与待扩展的结点重复(位于队列头指针之后),则比较两个结点的估价函数中g的大小,保留较小g值的结点。跳至第五步。 4)如果扩展出的新结点与队列中的结点不重复,则按照它的估价函数f大小将它插入队列中的头结点后待扩展结点的适当位置,使它们按从小到大的顺序排列,最后更新队列尾指针。 5)如果队列头的结点还可以扩展,直接返回第二步。否则将队列头指针指向下一结点,再返回第二步。 2、程序运行基本环境: 源程序所使用编程语言:C# 编译环境:VS2010,.net framework 4.0 运行环境:.net framework 4.0 3、程序运行界面 可使用程序中的test来随机生成源状态与目标状态 此停顿过程中按Enter即可使程序开始运行W(n)部分; 此停顿部分按Enter后程序退出; 4、无解问题运行情况 这里源程序中是先计算源状态与目标状态的逆序对的奇偶性是否一致来判断是否有解的。下面是无解时的运行画面: 输入无解的一组源状态到目标状态,例如: 1 2 3 4 5 6 7 8 0 1 2 3 4 5 6 8 7 0 运行画面如下: 5、性能比较 对于任一给定可解初始状态,状态空间有9!/2=181440个状态;当采用不在位棋子数作为启发函数时,深度超过20时,算法求解速度较慢; 其中启发函数P(n)与W(n)的含义如下: P(n): 任意节点与目标结点之间的距离; W(n): 不在位的将牌数; 源状态 目标状态 P(n) 生成节点数 W(n) 生成节点数 P(n) 扩展节点数 W(n) 扩展节点数 2 8 3 1 6 4 7 0 5 1 2 3 8 0 4 7 6 5 11 13 5 6 1 2 3 8 0 4 7 6 5 0 1 3 8 2 4 7 6 5 6 6 2 2 4 8 2 5 1 6 7 0 3 7 4 2 8 5 6 1 3 0 41 79 22 46 6 2 5 8 7 0 3 1 4 0 3 6 7 1 8 4 5 2 359 10530 220 6769 7 6 3 1 0 4 8 5 2 2 8 7 1 3 4 6 5 0 486 8138 312 5295 下图是解决随机生成的100中状态中,P(n)生成函数的生成节点与扩展节点统计图: 由上图可知,P(n)作为启发函数,平均生成节点数大约在1000左右,平均扩展节点数大约在600左右; 下图是解决随机生成的100中状态中,W(n)生成函数的生成节点与扩展节点统计图: 由上图可知,W (n)作为启发函数,平均生成节点数大约在15000左右,是P(n)作为启发函数时的平均生成节点的15倍;W (n)作为启发函数,平均扩展节点数大约在10000左右,是P(n)作为启发函数时的平均扩展节点的15倍; 下图是解决随机生成的100中状态中,两个生成函数的生成节点与扩展节点统计图: 由上述图表可以看到,将P(n)作为启发函数比将W(n)作为启发函数时,生成节点数与扩展节点数更稳定,相比较来说,采用P(n)作为启发函数的性能比采用W(n)作为启发函数的性能好。 6、源代码说明 1)AStar-EightDigital-Statistics文件夹:用来随机生成100个状态,并对这100个状态分别用P(n)与W(n)分别作为启发函数算出生成节点以及扩展节点,以供生成图表使用;运行界面如下: 2)Test文件夹:将0-8这9个数字随机排序,用来随机生成源状态以及目标状态的;运行界面如下: 3)AStar-EightDigital文件夹:输入源状态和目标状态,程序搜索出P(n)与W(n)分别作为启发函数时的生成节点数以及扩展节点数,并给出从源状态到目标状态的移动步骤;运行界面如下: 提高了运行速度的几处编码思想: 1、 在维护open以及close列表的同时,也维护一个类型为hashtable的open以及close列表,主要用来提高判断当前节点是否在open列表以及close列表中出现时的性能; 2、 对于每个状态,按照从左到右,从上到下,依次将数字拼接起来,形成一个唯一标识identify,通过该标识,可以直接判断两个状态是否是同一个状态,而不需要循环判断每个位置上的数字是否相等 3、 在生成每个状态的唯一标识identify时,同时计算了该状态的空格所在位置,通过空格所在位置,可以直接判断能否进行上移、下移、左移、右移等动作; 4、 只计算初始节点的h值,其它生成的节点的h值是根据当前状态的h值、移动的操作等计算后得出的,规则如下: a) 采用W(n)这种方式,不在位置的将牌数,共有以下3中情况: i. 该数字原不在最终位置上,移动后,在其最终位置上 这种情况下,生成的子节点的h值= 父节点的h值-1 ii. 该数字原在最终位置上,移动后,不在其最终位置上 这种情况下,生成的子节点的h值= 父节点的h值 +1 iii. 该数字原不在最终位置上,移动后,还是不在其最终位置上 这种情况下,生成的子节点的h值= 父节点的h值 iv. 该数字原在最终位置上,移动后,还在其最终位置 这种情况不存在 b) 采用P(n)这种方式,节点与目标距离,可通过下面3步完成 i. 首先计算在原位置时,与目标位置的距离,命名为Distance1 ii. 移动后,计算当前位置与目标位置的距离,命名为Distance2 iii. 计算子节点的h值: 子节点的h值 = 父节点的h值- Distance1+ Distance2 5、 在任意状态中的每个数字和目标状态中同一数字的相对距离就有9*9种,可以先将这些相对距离算出来,用一个矩阵存储,这样只要知道两个状态中同一个数字的位置,就可查出它们的相对距离,也就是该数字的偏移距离;例如在一个状态中,数字8的位置是3,在另一状态中位置是7,那么从矩阵的3行7列可找到2,它就是8在两个状态中的偏移距离。
以下是使用Python实现A*算法求解8数码问题的代码: ```python from queue import PriorityQueue class Puzzle: def __init__(self, start, goal): self.start = start self.goal = goal def get_manhattan_distance(self, state): distance = 0 for i in range(3): for j in range(3): value = state[i][j] if value != 0: goal_position = self.get_goal_position(value) distance += abs(i - goal_position[0]) + abs(j - goal_position[1]) return distance def get_goal_position(self, value): if value == 0: return (2, 2) else: return divmod(value - 1, 3) def get_successors(self, state): successors = [] zero_position = self.get_zero_position(state) for move in ["up", "down", "left", "right"]: new_state = self.get_new_state(state, zero_position, move) if new_state is not None: successors.append(new_state) return successors def get_zero_position(self, state): for i in range(3): for j in range(3): if state[i][j] == 0: return (i, j) def get_new_state(self, state, zero_position, move): new_state = [row[:] for row in state] i, j = zero_position if move == "up" and i > 0: new_state[i][j], new_state[i - 1][j] = new_state[i - 1][j], new_state[i][j] return new_state elif move == "down" and i < 2: new_state[i][j], new_state[i + 1][j] = new_state[i + 1][j], new_state[i][j] return new_state elif move == "left" and j > 0: new_state[i][j], new_state[i][j - 1] = new_state[i][j - 1], new_state[i][j] return new_state elif move == "right" and j < 2: new_state[i][j], new_state[i][j + 1] = new_state[i][j + 1], new_state[i][j] return new_state else: return None def solve(self): start_node = Node(self.start, None, None, 0, self.get_manhattan_distance(self.start)) open_list = PriorityQueue() open_list.put(start_node) closed_list = set() while not open_list.empty(): current_node = open_list.get() if current_node.state == self.goal: return self.get_path(current_node) closed_list.add(current_node.state) for successor in self.get_successors(current_node.state): if successor not in closed_list: successor_node = Node(successor, current_node, None, current_node.g + 1, self.get_manhattan_distance(successor)) if self.add_to_open_list(successor_node, open_list): open_list.put(successor_node) return None def add_to_open_list(self, successor_node, open_list): for node in open_list.queue: if node.state == successor_node.state and node.f <= successor_node.f: return False return True def get_path(self, node): path = [] while node is not None: path.append(node.state) node = node.parent return path[::-1] class Node: def __init__(self, state, parent, move, g, h): self.state = state self.parent = parent self.move = move self.g = g self.h = h self.f = g + h def __lt__(self, other): return self.f < other.f start_state = [[1, 2, 3], [4, 5, 6], [7, 8, 0]] goal_state = [[1, 2, 3], [4, 5, 6], [7, 8, 0]] puzzle = Puzzle(start_state, goal_state) path = puzzle.solve() if path is None: print("No solution found.") else: print("Solution found in", len(path) - 1, "moves:") for state in path: print(state) ``` 在此代码中,`Puzzle`类表示一个8数码问题,`Node`类表示搜索中的一个节点。`get_manhattan_distance`方法计算给定状态的曼哈顿距离,`get_goal_position`方法返回给定数字在目标状态中的位置,`get_successors`方法返回给定状态的所有后继状态,`get_zero_position`方法返回给定状态中数字0的位置,`get_new_state`方法返回给定状态中进行给定移动后的新状态,`solve`方法使用A*算法解决问题,`add_to_open_list`方法用于将节点添加到开放列表中,`get_path`方法返回从起始状态到给定节点状态的路径。在主函数中,我们创建一个`Puzzle`对象并调用`solve`方法来解决问题。如果找到了解决方案,则打印出路径;否则,打印“No solution found.”。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值