zoukankan      html  css  js  c++  java
  • A*搜索 概念

    启发式搜索:启发式搜索就是在状态空间中的搜索对每一个搜索的位置进行评估,得到最好的位置,再从这个位置进行搜索直到目标。这样可以省略大量无畏的搜索路径,提到了效率。在启发式搜索中,对位置的估价是十分重要的。采用了不同的估价可以有不同的效果。

    估价函数:从当前节点移动到目标节点的预估费用;这个估计就是启发式的。在寻路问题和迷宫问题中,我们通常用曼哈顿(manhattan)估价函数(下文有介绍)预估费用。

    A*算法与BFS:可以这样说,BFS是A*算法的一个特例。对于一个BFS算法,从当前节点扩展出来的每一个节点(如果没有被访问过的话)都要放进队列进行进一步扩展。也就是说BFS的估计函数h永远等于0,没有一点启发式的信息,可以认为BFS是“最烂的”A*算法。

    选取最小估价:如果学过数据结构的话,应该可以知道,对于每次都要选取最小估价的节点,应该用到最小优先级队列(也叫最小二叉堆)。在C++的STL里有现成的数据结构priority_queue,可以直接使用。当然不要忘了重载自定义节点的比较操作符。

    A*算法的特点:A*算法在理论上是时间最优的,但是也有缺点:它的空间增长是指数级别的。

    启发函数:f=g+h;其中g是起点到当前结点的直线距离,h是当前结点到目的结点的某种度量函数,在本题中采用曼哈顿距离。

    选择路径中经过哪个方格的关键是下面这个等式:F = G + H这里:

    • G = 从起点A,沿着产生的路径,移动到网格上指定方格的移动耗费。
    • H = 从网格上那个方格移动到终点B的预估移动耗费。这经常被称为启发式的,可能会让你有点迷惑。这样叫的原因是因为它只是个猜测。我们没办法事先知道路径的长度,因为路上可能存在各种障碍(墙,水,等等)。

    4,4,1 A*算法步骤为:

      • 把起始格添加到开启列表。
      • 重复如下的工作:
        • 寻找开启列表中F值最低的格子。我们称它为当前格。
        • 把它切换到关闭列表。
        • 对相邻的格中的每一个?
          • 如果它不可通过或者已经在关闭列表中,略过它。反之如下。
          • 如果它不在开启列表中,把它添加进去。把当前格作为这一格的父节点。记录这一格的F,G,和H值。
          • 如果它已经在开启列表中,用G值为参考检查新的路径是否更好。更低的G值意味着更好的路径。如果是这样,就把这一格的父节点改成当前格,并且重新计算这一格的G和F值。如果你保持你的开启列表按F值排序,改变之后你可能需要重新对开启列表排序。
        • 停止,当你
          • 把目标格添加进了关闭列表,这时候路径被找到,或者
          • 没有找到目标格,开启列表已经空了。这时候,路径不存在。
      • 保存路径。从目标格开始,沿着每一格的父节点移动直到回到起始格。这就是你的路径。
    //Poj2243
    #include<iostream>
    #include<queue>
    #include<stdlib.h>
    using namespace std;
    char ss[3];
    char ee[3];
    typedef struct node
    {
        int x;
        int y;
        int steps;
        int g;
        int h;
        int f;
        friend bool operator < (const node & a,const node &b);
    }node;
    inline bool operator < (const node & a,const node &b)
    {
        return a.f>b.f;
    }
    int d[8][2]={{-2,1},{-2,-1},{-1,-2},{-1,2},{2,-1},{2,1},{1,-2},{1,2}};
    int visited[8][8];
    node s;
    node e;
    int in(node n)
    {
        if(n.x<0||n.y<0||n.x>7||n.y>7)
            return 0;
        return 1;
    }
    int Heuristic(const node &a){
        return (abs(a.x-e.x)+abs(a.y-e.y))*10;
    }//曼哈顿(manhattan)估价函数
    priority_queue<node> q; //最小优先级队列(开启列表) 这里有点优化策略,因为我发现如果把q
     //放在Astar函数里头的话,代码跑起来是157MS,放在外面的话是47MS,有显著的区别
    int Astar()
    {
        while(!q.empty())q.pop();
        memset(visited,0,sizeof(visited));
        q.push(s);
        while(!q.empty())
        {
            node front=q.top();
            node t;
            q.pop();
            visited[front.x][front.y]=1;
            if(front.x==e.x && front.y==e.y)
                return front.steps;
            for(int i=0;i<8;i++){
                t.x=front.x+d[i][0];
                t.y=front.y+d[i][1];
                if(in(t) && visited[t.x][t.y]==0){
                    t.g=23+front.g;
                    t.h=Heuristic(t);
                    t.f=t.g+t.h;
                    t.steps=front.steps+1;
                    q.push(t);
                }
            }
        }
    }
    int main(int argc, char *argv[])
    {
        //freopen("in.txt","r",stdin);
        while(scanf("%s %s",ss,ee)==2)
        {
            s.x=ss[0]-'a';
            s.y=ss[1]-'1';
            e.x=ee[0]-'a';
            e.y=ee[1]-'1';
            s.steps=0;
            s.g=0;
            s.h=Heuristic(s);
            s.f=s.g+s.h;
            if(s.x==e.x&&s.y==e.y)
                printf("To get from %s to %s takes 0 knight moves.
    ",ss,ee);
            else
                printf("To get from %s to %s takes %d knight moves.
    ",ss,ee,Astar());
        }
        return 0;
    }
  • 相关阅读:
    堆(优先队列)模板
    线段树(递归)模板
    快速幂(含二阶方阵类)模板
    ACM/CF赛制getstart模板
    [Andrew Stankevich's Contest#21] Lempel-Ziv Compression
    [NOIP2001]Car的旅行路线
    [NOIP2007] 矩阵取数游戏
    [NOIP2010] 关押罪犯
    [NOIP1999] 拦截导弹
    设计模式之第1章-工厂方法模式(Java实现)
  • 原文地址:https://www.cnblogs.com/FuTaimeng/p/5600762.html
Copyright © 2011-2022 走看看