7.4图的应用:骑士周游问题

在一个国际象棋棋盘上,一个棋子“马” (骑士),按照“马走日”的规则,从一个格子出发,要走遍所有棋盘格恰好一次。把一个这样的走棋序列称为一次“周游


7.4图的应用:骑士周游问题_第1张图片
image.png

在8×8的国际象棋棋盘上,合格的“周游”数量有1.305×1035这么多,走棋过程中
失败的周游就更多了
采用图搜索算法,是解决骑士周游问题最容易理解和编程的方案之一

解决方案还是分为两步:
❖首先将合法走棋次序表示为一个图
❖采用图搜索算法搜寻一个长度为(行×列-1)的路径,路径上包含每个顶点恰一次

将棋盘和走棋步骤构建为图的思路
将棋盘格作为顶点,按照“马走日”规则的走棋步骤作为连接边,建立每一个棋盘格的所有合法走棋步骤能够到达的棋盘格关系图


7.4图的应用:骑士周游问题_第2张图片
image.png

合法走棋位置函数

def genLegalMoves(x, y , bdSize):
    newMoves = []
    moveOffsets =[(-1,-2),(-1,2),(-2,-1),(-2,1),
                   ( 1,-2),( 1,2),( 2,-1),( 2,1)]
    for i in moveOffsets:
        newX = x + i[0]
        newY = y + i[1]
        if legalCoord(newX, bdSize) and legalCoord(newX, bdSize):
            newMoves.append((newX, newY))

        return newMoves

#是否超出棋盘范围
def legalCoord(x, bdSize):
    if x >= 0 and x < bdSize:
        return True
    else:
        return False

构建走棋关系图

def posToNodeID(row, col, bdSize):
    return row*bdSize + col

def knightGraph(bdSize):
    ktGraph = Graph()
    for row in range(bdSize):
        for col in range(bdSize):
            nodeid = posToNodeID(row, col, bdSize)
            newPositions = genLegalMoves(row, col, bdSize)
            for e in newPositions:
                nid = posToNodeID(e[0], e[1], bdSize)
                ktGraph.addEdge(nodeid, nid)
    return ktGraph

算法实现
❖用于解决骑士周游问题的图搜索算法是深度优先搜索(Depth First Search)
❖相比前述的广度优先搜索,其逐层建立搜索树的特点
❖深度优先搜索是沿着树的单支尽量深入向下搜索
(如果到无法继续的程度还未找到问题解就回溯上一层再搜索下一支)

DFS的两个实现算法
一个DFS算法用于解决骑士周游问题,其特点是每个顶点仅访问一次
另一个DFS算法更为通用,允许顶点被重复访问,可作为其它图算法的基础

深度优先搜索解决骑士周游的关键思路:
如果沿着单支深入搜索到无法继续(所有合法移动都已经被走过了)时,路径长度还没有达到预定值(8×8棋盘为63)那么就清除颜色标记,返回到上一层换一个分支继续深入搜索
引入一个栈来记录路径,并实施返回上一层的回溯操作

def knightTour(n,path,u,limit):#n:层次,path:路径,u:当前顶点,limit:搜索总深度
    u.setColor('gray')
    path.append(u)#当前顶点加入路径
    done = False

    if n < limit:
        nbrLst = list(u.getConnections())#对所有合法移动逐一深入
        i = 0
        while i < len(nbrLst) and not done:
            if nbrLst[i].getColor() == 'white':
                done = knightTour(n+1,path,nbrLst[i],limit)#nbrLst[i]递归深入
            i+=1

        if not done:
            path.pop()
            u.setColor('white')

    else:
        done = True

    return done
7.4图的应用:骑士周游问题_第3张图片
image.png

7.4图的应用:骑士周游问题_第4张图片
image.png

上述算法的性能高度依赖于棋盘大小:
就5×5棋盘,约1.5秒可以得到一个周游路径,但8×8棋盘,则要半个小时以上才能得到一个解
目前实现的算法,其复杂度为O(kn),其中n是棋盘格数目,这是一个指数时间复杂度的算法!其搜索过程表现为一个层次为n的树

骑士周游算法改进
对nbrList的灵巧构造,以特定方式排列顶点访问次序
可以使得8×8棋盘的周游路径搜索时间降低到秒级!
这个改进算法被特别以发明者名字命名:Warnsdorff算法
初始算法中nbrList,直接以原始顺序来确定深度优先搜索的分支次序
新的算法,仅修改了遍历下一格的次序,将u的合法移动目标棋盘格排序为:具有最少合
法移动目标的格子优先搜索

def orderbyAvail(n):
    resultLst = []
    for v in n.getConnections():
        if v.getColor() == 'white':
            c = 0
            for w in v.getConnections():
                if w.getColor() == 'white':
                    c+=1
            resultLst.append((c,v))

    resultLst.sort(key=lambda x:x[0])
    return [y[1] for y in resultLst]

完整代码:
Graph类

import sys


class Vertex:
    def __init__(self,key):
        self.id = key
        self.connectedTo = {}
        self.color = 'white'
        self.dist = sys.maxsize
        self.pred = None

    def setColor(self, color):
        self.color = color

    def setDistance(self, d):
        self.dist = d

    def setPred(self, p):
        self.pred = p

    def getPred(self):
        return self.pred

    def getDistance(self):
        return self.dist

    def getColor(self):
        return self.color

    def addNeighbor(self,nbr,weight):
        self.connectedTo[nbr] = weight

    def __repr__(self):
        return f'{self.id} connectto:' + str([v.id for v in self.connectedTo])

    def getConnections(self):
        return self.connectedTo.keys()

    def getId(self):
        return self.id

    def getWeight(self,nbr):
        return self.connectedTo[nbr]



class Graph:
    def __init__(self):
        self.verList = {}
        self.num = 0

    def addVertex(self, key):
        self.num += 1
        v = Vertex(key)
        self.verList[key] = v
        return v

    def getVertex(self,key):
        if key in self.verList:
            return self.verList[key]
        else:
            return None

    def __contains__(self, key):
        return key in self.verList

    def addEdge(self,fromkey,tokey,cost=0):
        if fromkey not in self.verList:
            self.addVertex(fromkey)
        if tokey not in self.verList:
            self.addVertex(tokey)

        self.verList[fromkey].addNeighbor(self.verList[tokey],cost)

    def getVertices(self):
        return self.verList.keys()

    def __iter__(self):
        return iter(self.verList.values())


if __name__ == '__main__':
    g = Graph()
    for i in range(6):
        g.addVertex(i)

        print(g.addVertex(i))
from basicGraph import Graph


def posToNodeID(row, col, bdsize):
    return row*bdsize + col


def legalCoord(newX, bdsize):
    if newX >= 0 and newX < bdsize:
        return True
    else:
        return False


def genLegalMoves(x, y, bdsize):
    movesLst = []
    moveOffsets = [(-2,-1),(-1,-2),(-2,1),(1,-2),(2,-1),(-1,2),(1,2),(2,1)]
    for tup in moveOffsets:
        newX = x + tup[0]
        newY = y + tup[1]
        if legalCoord(newX,bdsize) and legalCoord(newY,bdsize):
            movesLst.append((newX,newY))

    return movesLst


def knightGraph(bdsize):
    g = Graph()
    for row in range(bdsize):
        for col in range(bdsize):
            nodeId = posToNodeID(row,col,bdsize)
            newPosLst = genLegalMoves(row,col,bdsize)
            for pos in newPosLst:
                targetNodeID = posToNodeID(pos[0],pos[1],bdsize)
                g.addEdge(nodeId,targetNodeID)
    return g

def orderbyAvail(n):
    resultLst = []
    for v in n.getConnections():
        if v.getColor() == 'white':
            c = 0
            for w in v.getConnections():
                if w.getColor() == 'white':
                    c+=1
            resultLst.append((c,v))

    resultLst.sort(key=lambda x:x[0])
    return [y[1] for y in resultLst]


def knightTour(n,path,u,limit):
    u.setColor('gray')
    path.append(u)
    done = False

    if n < limit:
        # nbrLst = list(u.getConnections())
        nbrLst = list(orderbyAvail(u))
        i = 0
        while i < len(nbrLst) and not done:
            if nbrLst[i].getColor() == 'white':
                done = knightTour(n+1,path,nbrLst[i],limit)
            i+=1

        if not done:
            path.pop()
            u.setColor('white')

    else:
        done = True

    return done



if __name__ == '__main__':
    g = knightGraph(8)
    # for i in g:
    #     print(i)

    path = []
    startVertex = g.getVertex(4)
    knightTour(0,path,startVertex,63)
    for node in path:
        print(node.getId(),end=' ')


你可能感兴趣的:(7.4图的应用:骑士周游问题)