【人工智能实验】A*算法求解8数码问题 golang

人工智能经典问题八数码求解

实际上是将求解转为寻找最优节点的问题,算法流程如下:

  1. 求非0元素的逆序数的和,判断是否有解
  2. 将开始状态放到节点集,并设置访问标识位为true
  3. 从节点集中取出h(x)+g(x)最小的节点
  4. 判断取出的节点的状态是不是最终状态,如果是的话则回溯打印
  5. 找出取出的节点的状态中的0的位置
  6. 对取出的节点进行move操作,包含up down left right
  7. 如果move后的状态的访问标识位为false,则添加。否则什么都不做

需要注意:节点的数据结构如下

  • 状态:int数组
  • h(x):当前节点的状态到目标状态的距离
  • g(x):当前节点的状态到初始状态的距离
  • 动作:到当前节点所进行的move类型
  • 父节点:记录上一个状态,方便回溯打印

使用go语言实现如下

  • main.go

    package main
    
    import (
    	"container/heap"
    	"github.com/gookit/color"
    	"log"
    	"os"
    	"os/signal"
    	"syscall"
    )
    
    var (
    	start  = []int{2, 8, 3, 1, 6, 4, 7, 0, 5}
    	target = []int{1, 2, 3, 8, 0, 4, 7, 6, 5}
    )
    var (
    	movables    = []string{"up", "down", "left", "right"}
    	moveOffsets = map[string]int{"up": -3, "down": 3, "left": -1, "right": 1}
    )
    var (
    	visited = make(map[string]bool)
    )
    
    func main() {
    	color.BgCyan.Println("Y02114562")
    	printFun := func(list []int) {
    		for _, i := range list {
    			color.BgLightCyan.Print(i, ",")
    		}
    		color.BgLightCyan.Print("\n")
    	}
    	printFun(start)
    	printFun(target)
    	if reverseSum(start) != reverseSum(target) {
    		log.Fatal("不可解")
    	}
    	path, steps := solve(start)
    	if steps == -1 {
    		log.Fatal("No solution")
    	}
    	color.BgGreen.Println("只需:", steps, "步")
    	color.BgGreen.Println("操作:", path)
    
    	quit := make(chan os.Signal, 1)
    	signal.Notify(quit, syscall.SIGINT, syscall.SIGTERM)
    	<-quit
    }
    
    // 启发函数:h(x) 从当前状态到目标的距离
    func manhattanDistance(state []int) int {
    	distance := 0
    
    	for i := 0; i < 9; i++ {
    		if state[i] != 0 {
    			row1, col1 := i/3, i%3
    			// 遍历所有不为0的点,计算他与他的目标位置的曼哈顿距离
    			for j := 0; j < 9; j++ {
    				if state[i] == target[j] {
    					row2, col2 := j/3, j%3
    					distance += abs(row1-row2) + abs(col1-col2)
    					break
    				}
    			}
    		}
    	}
    	return distance
    }
    
    // 启发式搜索:八数码问题求解
    func solve(start []int) ([]string, int) {
    	// 创建起始节点
    	startNode := &Node{
    		State:     start,
    		Heuristic: manhattanDistance(start),
    		G:         0,
    		PrevMove:  "",
    		PrevNode:  nil,
    	}
    
    	// 创建优先队列
    	pq := make(PriorityQueue, 0)
    	heap.Init(&pq)
    	heap.Push(&pq, startNode)
    
    	visited[listToString(startNode.State)] = true
    
    	// A*搜索
    	for pq.Len() > 0 {
    		currentNode := heap.Pop(&pq).(*Node)
    
    		// 到达目标状态,返回路径
    		if listToString(currentNode.State) == listToString(target) {
    			path := make([]string, 0)
    			for currentNode.PrevNode != nil {
    				path = append(path, currentNode.PrevMove)
    				currentNode = currentNode.PrevNode
    			}
    			return func(slice []string) ([]string, int) {
    				for i, j := 0, len(slice)-1; i < j; i, j = i+1, j-1 {
    					slice[i], slice[j] = slice[j], slice[i]
    				}
    				return slice, len(path)
    			}(path)
    		}
    
    		zeroIndex := func(state []int) int {
    			for i, num := range state {
    				if num == 0 {
    					return i
    				}
    			}
    			return -1
    		}(currentNode.State)
    		for _, move := range movables {
    			if canMove(move, zeroIndex) {
    				newState := make([]int, len(currentNode.State))
    				copy(newState, currentNode.State)
    
    				newZeroIndex := zeroIndex + moveOffsets[move]
    				newState[zeroIndex], newState[newZeroIndex] = newState[newZeroIndex], newState[zeroIndex]
    
    				// 创建新节点
    				newNode := &Node{
    					State:     newState,
    					Heuristic: manhattanDistance(newState),
    					G:         currentNode.G + 1,
    					PrevMove:  move,
    					PrevNode:  currentNode,
    				}
    
    				// 如果新状态未被访问,则加入优先队列和已访问集合
    					if !visited[listToString(newState)] {
    						heap.Push(&pq, newNode)
    						visited[listToString(newState)] = true
    					}
    			}
    		}
    	}
    
    	// 没有找到解
    	return nil, -1
    }
    
  • node.go

    package main
    
    // Node 节点结构体
    type Node struct {
    	State     []int  // 当前状态
    	Heuristic int    // 启发函数值
    	G         int    // 初始节点到当前节点
    	PrevMove  string // 上一步移动的方向
    	PrevNode  *Node  // 上一步的节点
    }
    
    // PriorityQueue 优先队列
    type PriorityQueue []*Node
    
    // Len 优先队列的方法:计算长度
    func (pq PriorityQueue) Len() int {
    	return len(pq)
    }
    
    // Less 优先队列的方法:比较优先级
    func (pq PriorityQueue) Less(i, j int) bool {
    	return pq[i].Heuristic+pq[i].G < pq[j].Heuristic+pq[j].G
    }
    
    // Swap 优先队列的方法:交换元素
    func (pq PriorityQueue) Swap(i, j int) {
    	pq[i], pq[j] = pq[j], pq[i]
    }
    
    // Push 优先队列的方法:向队列中插入元素
    func (pq *PriorityQueue) Push(x interface{}) {
    	node := x.(*Node)
    	*pq = append(*pq, node)
    }
    
    // Pop 优先队列的方法:从队列中弹出元素
    func (pq *PriorityQueue) Pop() interface{} {
    	old := *pq
    	n := len(old)
    	node := old[n-1]
    	*pq = old[0 : n-1]
    	return node
    }
    
  • tool.go

    package main
    
    import "fmt"
    
    // 辅助函数:判断是否可移动
    func canMove(move string, zeroIndex int) bool {
    	if move == "up" && zeroIndex >= 3 {
    		return true
    	}
    	if move == "down" && zeroIndex <= 5 {
    		return true
    	}
    	if move == "left" && zeroIndex%3 != 0 {
    		return true
    	}
    	if move == "right" && zeroIndex%3 != 2 {
    		return true
    	}
    	return false
    }
    
    // 辅助函数:将[]int转换为字符串
    func listToString(state []int) string {
    	str := ""
    	for _, num := range state {
    		str += fmt.Sprintf("%d", num)
    	}
    	return str
    }
    
    // 辅助函数:求除了0之外的逆序和
    func reverseSum(arr []int) bool {
    	sum := 0
    
    	for i := 1; i < len(arr); i++ {
    		if arr[i] != 0 {
    			for j := 0; j < i; j++ {
    				if arr[j] > arr[i] {
    					sum++
    				}
    			}
    		}
    	}
    
    	return sum%2 != 0
    }
    
    // 辅助函数:计算绝对值
    func abs(num int) int {
    	if num < 0 {
    		return -num
    	}
    	return num
    }
    

运行效果

Y02114562

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
include using namespace std; struct node{ int nodesun[4][4]; int pre; //上一步在队列中的位置 int flag ; //步数标识,表示当前的步数为有效的 int value; //与目标的差距 int x,y; //空格坐标 }queue[1000]; //移动方向数组 int zx[4]={-1,0,1,0}; int zy[4]={0,-1,0,1}; //当前步数 int top; int desti[4][4];//目标状态 int detect(struct node *p)//检查是否找到 {int i,j; for(i=1;i<4;i++) for(j=1;jnodesun[i][j]!=desti[i][j]) return 0; return 1; } //打印 void printlj() {int tempt; int i,j; tempt=top; while(tempt!=0) { for(i=1;i<4;i++) for(j=1;j<4;j++) {cout<<queue[tempt].nodesun[i][j]; if(j==3) cout<<" "<<endl; } tempt=queue[tempt].pre; } } //现在状态与目标状态有多少个不同位置 int VALUE(struct node *p) {int count=0; int i,j; for(i=1;i<4;i++) for(j=1;jnodesun[i][j]!=desti[i][j]) count++; return count; } void main() { //初始化 int i,j,m,n,f; int min=10; int temp,find=0,minnumber; top=1; for(i=1;i<4;i++) for(j=1;j<4;j++) {cout<<"请输入第"<<i<<"行"<<"第"<<j<<"列的值"<>temp; queue[1].nodesun[i][j]=temp; } cout<<"请输入初始状态的空格的位置(行)"<>temp; queue[1].x=temp; cout<<"请输入初始状态的空格的位置(列)"<>temp; queue[1].y=temp; queue[1].value=VALUE(&queue[1]); queue[1].pre=0; //上一步在队列中的位置 queue[1].flag=0; //目标状态 for(i=1;i<4;i++) for(j=1;j<4;j++) {cout<<"请输入目标状态第"<<i<<"行"<<"第"<<j<<"列的值"<>temp; desti[i][j]=temp; } //根据估价函数 while(!find&&top>0) { for(i=1;i<=top;i++) //////////////////////////////////////////// //min为上一图中与目标图有多少个元素不相同,queue[i]为当前图与目标图有多少个元素不相同通过这两个数的比较,就可以得出当前图较之上一图向目标图接近同时把当前的i记录下来进行下一步比较 {if(queue[i].value<min&&queue[i].flag==0) {minnumber=i;// min=queue[i].value; //还有多少不同的位数 } } queue[minnumber].flag=1; //表示此位有效 ////////////////////////////////////// // for(f=0;f=1&&i=1&&j<=3) {top++; ///////////////////////////////////////////// //位置交换 queue[top]=queue[minnumber]; queue[top].nodesun[m][n]=queue[minnumber].nodesun[i][j]; queue[top].nodesun[i][j]=0; /////////////////////////////////////// //空格移动方向 queue[top].x=i; queue[top].y=j; /////////////////////////////////////// queue[top].pre=minnumber; //上一步在队列中的位置 queue[top].value=VALUE(&queue[top]); //有多少位与目标不同 queue[top].flag=0; //标识位初始化 if(detect(&queue[top])) //检查是否为目标 {printlj(); //打印 find=1; //设找到标识位 break; } } } } }
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在两个状态中的偏移距离。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

有头发的琦玉

打点钱,我会再努力的

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值