这是本蒟蒻在洛谷通过的第一道 “提高+/省选-” 的 水 蓝题。
八数码难题是一个非常经典的BFS的题,并且P1379这道题即使是最朴素的BFS也是能够通过的,但是结果是这样:
emmm这时长,这空间占用,貌似连朴素BFS都不是的样子
所以为了让程序更加高效,让技术更好的服务于人类,让自己的记录在最优解的排名上更靠前
我们能想到A*算法,也就是启发式搜索。(什么,你不懂啥是A*?赶紧去查资料啊)
利用STL的priority_queue和map来实现启发搜索和去重。
(STL大法好!)
最后实现的结果是
启发式搜索的效率真的高,朴素BFS就相当于估价一直都是0的A*,说白了就是爆搜。
97分那一次是因为我第一次使用的估价函数不够好,甚至导致了错误答案,于是我重新写了一个估价函数,不仅过了还让效率更高了~
好了,放上本蒟蒻的代码:
#include<bits/stdc++.h>
using namespace std;
struct node{
string s;//存储当前状态
int step;//存储到该状态所用的步数
int h;
int zpos;//记录一下0的位置
}indata,n;
bool operator < (node x,node y)
{
return x.h>y.h;
//重载一下运算符,不然优先队列不认识咱这结构体。
}
const string goal="123804765";
const int expos[9][4]={{0,0,1,3},
{1,0,2,4},
{2,1,2,5},
{0,3,4,6},
{1,3,5,7},
{2,4,5,8},
{3,6,7,6},
{4,6,8,7},
{5,7,8,8}};
/*
expos[i][0]的含义是:当前'0'位于i这个位置,可以与expos[i][0]的位置进行交换
并且我们将不能交换的位置设置为i本身,这样可以在下面进行判断
我们设置盘面坐标为
0 1 2
3 4 5
6 7 8
*/
const int gpos[9][2]={{1,1},{0,0},{0,1},{0,2},{1,2},{2,2},{2,1},{2,0},{1,0}};
/*
此处表示的是位置i上的数的正确位置在二维坐标系下的坐标
(0,0) (0,1) (0,2) 1 2 3
(1,0) (1,1) (1,2) => 8 0 4
(2,0) (2,1) (2,2) 7 6 5
*/
priority_queue<node> qu;
map<string,bool> m;//利用map查重
int getH(string s)
{
//非常重要的估价函数
int cnt=0;
for(int i=0;i<9;i++)
{
int x=i/3,y=i%3;
cnt+=abs(x-gpos[s[i]-'0'][0]);
cnt+=abs(y-gpos[s[i]-'0'][1]);
}
//通过计算当前位置上的数与目标位置的曼哈顿距离来作为估价
return cnt;
}
/*
这是我第一次交的时候用的估价函数
结果 #12 WA了,97分。
int getH(string s)
{
int cnt=0;
for(int i=0;i<9;i++)
if(s[i]!=goal[i])
cnt++;
return cnt;
}
*/
int main()
{
cin>>indata.s;//输入初始状态
indata.step=0;//步数置零
indata.h=indata.step+getH(indata.s);//求估价
for(int i=0;i<9;i++)
if(indata.s[i]=='0')
{indata.zpos=i;break;}//求0的位置
qu.push(indata);
m[indata.s]=true;//不要忘了把初始状态标记为true
//开始愉快的BFS
while(!qu.empty())
{
n=qu.top();
qu.pop();
if(n.s==goal)
{
cout<<n.step<<endl;
return 0;
}
int px=n.zpos,nx;//px表示当前0的位置
for(int k=0;k<4;k++)
{
nx=expos[px][k];//然后枚举要与之交换的位置nx
if(nx!=px)//这里就体现了上面将不能交换的位置设置成自己的好处了,nx!=px就说明交换是合法的
{
string s1=n.s;
swap(s1[px],s1[nx]);
if(m.count(s1)==0)//如果map中没有出现该状态就把该状态存入队列
{
node newn;
newn.s=s1;
newn.step=n.step+1;
newn.zpos=nx;//因为是交换了0的位置,所以新的0的位置就是nx
newn.h=newn.step+getH(newn.s);
qu.push(newn);
m[newn.s]=true;//不要忘记设置为已经访问过
}
}
}
}
}