A*算法:P5507 机关 题解

题目入口

这个题有很多解法,比如双向广搜。现在使用A*算法,这个算法属于启发式搜索的一类。

对于一个状态x设从起点到这个状态的代价是g(x),设h*(x)为理想情况下从x到终点的代价的估计值。注意这个值一定是理想状态下的,一定不会大于最终的实际距离。

设f*(x)=g(x)+h*(x)为状态的估值函数,在搜索时,找到f*值最小的状态(优先队列),优先对其进行搜索,配合最优性枝剪,可以有效地剪去较差分枝,提高算法效率。

对于这道题,g(x)就是已走的步数,h*(x)可设计为转到目标状态的步数除以2,因为在理想情况下转一次按钮有可能会让另一个按钮也想目标状态转一次。当然,在绝大部分时候,h(x)要远大于 h*(x),所以可以考虑给h*(x)加上系数,让它更逼近真实的步数,至于这个系数要有多大,就是比较玄学的事情了,太大有可能会出错,太小程序有可能会超时。

程序本身是BFS算法,不过使用的是优先队列维护优先级,f*(x)小的优先。从初始状态开始,扩展到下一个状态时发现目标状态,则可以直接输出结果;计算对应的f*(x),放入优先队列中。接着从优先队列中找到f*(x)最小的状态来进行扩展。

看上去挺复杂的,其实就比常规的BFS多了几行代码,代码如下:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=1<<24;
int g[N],nxt[14][6],fa[N],ans[30],choice[N];
struct node{
    int state;
    double F;
    node(int s):state(s){
        double h=0;
        F=0;
        for(int i=0;i<=11;i++)
            if((s>>(i<<1))&3)
                h+=4-((s>>(i<<1))&3);
        F=h*0.6+g[s];//给h加的系数,多次测试发现系数在0.6到1.1之间可以过全部检查点
    }
    bool operator<(const node &y)const{
        return F>y.F;
    }
};
priority_queue<node>q;
int main (){
    int button,Start=0;
    for(int i=0;i<=11;i++){
        cin>>button;
        Start|=(button-1)<<(i<<1);
        for(int j=0;j<=3;j++){
            int x;
            cin>>x;
            nxt[i][j]=x-1;
        }
    }
    q.push(node(Start));
    g[Start]=0;
    while(!q.empty()){
        int state=q.top().state;
        q.pop();
        if(state==0) break;//到达结果
        int si,sNxt,nx,nextState;
        for(int i=0;i<=11;i++){
            si=(state>>(i<<1))&3;//第i个按钮的状态
            nx=nxt[i][si];
            sNxt=(state>>(nx<<1))&3;//受牵连的按钮
            nextState = state ^ (si << (i << 1)) ^ (((si + 1) & 3) << (i << 1)) ^ (sNxt << (nx << 1)) ^ (((sNxt + 1) & 3) << (nx << 1));//旋转后的状态
            if(!g[nextState]){//如果没走过,更新状态
                g[nextState]=g[state]+1;
                fa[nextState]=state;
                choice[nextState]=i+1;
                q.push(node(nextState));
            }
        }
    }
    int cnt=0,state=0;
    while(state!=Start){//回溯
        ans[++cnt]=choice[state];
        state=fa[state];
    }
    cout<<cnt<<endl;
    for(int i=cnt;i;i--)
        cout<<ans[i]<<" ";
    return 0;
}
相关推荐
NAGNIP9 小时前
轻松搞懂全连接神经网络结构!
人工智能·算法·面试
NAGNIP10 小时前
一文搞懂激活函数!
算法·面试
董董灿是个攻城狮10 小时前
AI 视觉连载7:传统 CV 之高斯滤波实战
算法
爱理财的程序媛16 小时前
openclaw 盯盘实践
算法
MobotStone19 小时前
Google发布Nano Banana 2:更快更便宜,图片生成能力全面升级
算法
颜酱1 天前
队列练习系列:从基础到进阶的完整实现
javascript·后端·算法
用户5757303346241 天前
两数之和:从 JSON 对象到 Map,大厂面试官到底在考察什么?
算法
程序猿追1 天前
“马”上行动:手把手教你基于灵珠平台打造春节“全能数字管家”
算法
ZPC82102 天前
docker 镜像备份
人工智能·算法·fpga开发·机器人