八数码

http://www.dakaren.com/blog/post/3.html

八数码问题

一.八数码问题
八数码问题也称为九宫问题。在3×3的棋盘，摆有八个棋子，每个棋子上标有1至8的某一数字，不同棋子上标的数字不相同。棋盘上还有一个空格，与空格相邻的棋子可以移到空格中。要求解决的问题是：给出一个初始状态和一个目标状态，找出一种从初始转变成目标状态的移动棋子步数最少的移动步骤。
所谓问题的一个状态就是棋子在棋盘上的一种摆法。棋子移动后，状态就会发生改变。解八数码问题实际上就是找出从初始状态到达目标状态所经过的一系列中间过渡状态。
八数码问题一般使用搜索法来解。
搜索法有广度优先搜索法、深度优先搜索法、A*算法等。这里通过用不同方法解八数码问题来比较一下不同搜索法的效果。
二.搜索算法基类
1.八数码问题的状态表示
八数码问题的一个状态就是八个数字在棋盘上的一种放法。每个棋子用它上面所标的数字表示，并用0表示空格，这样就可以将棋盘上棋子的一个状态存储在一个一维数组p[9]中，存储的顺序是从左上角开始，自左至右，从上到下。也可以用一个二维数组来存放。
2.结点
搜索算法中，问题的状态用结点描述。结点中除了描述状态的数组p[9]外，还有一个父结点指针last，它记录了当前结点的父结点编号，如果一个结点v是从结点u经状态变化而产生的，则结点u就是结点v的父结点，结点v的last记录的就是结点u的编号。在到达目标结点后，通过last 可以找出搜索的路径。
3.类的结构
在C++中用类来表示结点，类将结点有关的数据操作封装在一起。
不同的搜索算法具有一定共性，也有各自的个性，因此这里将不同搜索算法的共有的数据和功能封装在一个基类中，再通过继承方式实现不同的搜索算法。
4.结点扩展规则
搜索就是按照一定规则扩展已知结点，直到找到目标结点或所有结点都不能扩展为止。
八数码问题的结点扩展应当遵守棋子的移动规则。按照棋子移动的规则，每一次可以将一个与空格相邻棋子移动到空格中，实际上可以看作是空格作相反移动。空格移动的方向可以是右、下、左、上，当然不能移出边界。
棋子的位置，也就是保存状态的数组元素的下标。空格移动后，它的位置发生变化，在不移出界时，空格向右、下、左和上移动后，新位置是原位置分别加上1、3、-1、-3，如果将空格向右、下、左和上移动分别用0、1、2、3表示，并将-3、3、-1、1放在静态数组d[4]中，空格位置用spac表示，那么空格向方向i移动后，它的位置变为spac+d[i]。
空格移动所产生的状态变化，反映出来则是将数组p[]中，0的新位置处的数与0交换位置。
5.八数码问题的基类
八数码问题的基类及其成员函数的实现如下：

[cpp] view plaincopyprint?

1. #define Num 9  
2.   
3. class TEight  
4. {  
5. public:  
6.     TEight(){}  
7.     TEight(char *fname);  
8.     virtual void Search()=0;  
9. protected:  
10.     int p[Num];  
11.     int last,spac;  
12.     static int q[Num],d[],total;  
13.     void Printf();  
14.     bool operator==(const TEight &T);  
15.     bool Extend(int i);  
16. };  
17. int TEight::q[Num];  
18. int TEight::d[]={1,3,-1,-3};  
19. int TEight::total=0;  
20.   
21. TEight::TEight(char *fname)  
22. {  
23.     ifstream fin;  
24.     fin.open(fname,ios::in | ios::nocreate);  
25.     if(!fin)  
26.     {  
27.         cout<<"不能打开数据文件!"<<endl;  
28.         return;  
29.     }  
30.     for(int i=0;i<num;)  
31.         fin>>p[i++];  
32.     fin>>spac;  
33.     for(i=0;i<num;)  
34.         fin>>q[i++];  
35.     fin.close();  
36.     last=-1;  
37.     total=0;  
38. }  
39.   
40. void TEight::Printf()  
41. {  
42.     ofstream fout;  
43.     fout.open(".\Eightr.txt",ios::ate);  
44.     fout<<setw(4)<<total++<<"t";  
45.     for(int i=0;i<num;)  
46.         fout<<setw(2)<<p[i++];  
47.     fout<<endl;  
48.     fout.close();  
49. }  
50.   
51. bool TEight::operator==(const TEight &T)  
52. {  
53.     for(int i=0;i<num;)  
54.         if(T.p[i]!=p[i++])  
55.             return 0;  
56.     return 1;  
57. }  
58.   
59. bool TEight::Extend(int i)  
60. {  
61.     if(i==0 && spac%3==2 || i==1 && spac>5   
62.         || i==2 && spac%3==0 || i==3 && spac<3)  
63.         return 0;  
64.     int temp=spac;  
65.     spac+=d[i];  
66.     p[temp]=p[spac];  
67.     p[spac]=0;  
68.     return 1;  
69. }  

数据文件的结构：
一共三行，第一行是用空格隔开的九个数字0~8，这是初始状态。第二行是一个数字，空格（数字0）的位置，第三行也是用空格隔开的九个数字0~8，这是目标状态。

三.线性表
搜索法在搜索过程中，需要使用一个队列存储搜索的中间结点，为了在找到目标结点后，能够找到从初始结点到目标结点的路径，需要保留所有搜索过的结点。另一方面，不同问题甚至同一问题的不同搜索方法中，需要存储的结点数量相差很大，所以这里采用链式线性表作为存储结构，同时，为适应不同问题，线性表设计成类模板形式。

[cpp] view plaincopyprint?

1. template class TList; //线性表前视定义  
2.   
3. template class TNode //线性表结点类模板  
4. {  
5.     friend class TList;  
6. public:  
7.     TNode(){}  
8.     TNode(const Type& dat);  
9. private:  
10.     TNode* Next;  
11.     Type Data;  
12. };  
13.   
14. template class TList  
15. {  
16. public:  
17.     TList(){Last=First=0;Length=0;} //构造函数  
18.     int Getlen()const{return Length;} //成员函数，返回线性表长度  
19.     int Append(const Type& T); //成员函数，从表尾加入结点  
20.     int Insert(const Type& T,int k); //成员函数，插入结点  
21.     Type GetData(int i); //成员函数，返回结点数据成员  
22.     void SetData(const Type& T,int k); //成员函数，设置结点数据成员   
23. private:  
24.     TNode *First,*Last; //数据成员，线性表首、尾指针   
25.     int Length; //数据成员，线性表长度  
26. };  
27.   
28. template int TList::Append(const Type& T)  
29. {  
30.     Insert(T,Length);  
31.     return 1;  
32. }  
33.   
34. template int TList::Insert(const Type& T,int k)  
35. {  
36.     TNode *p=new TNode;  
37.     p->Data=T;  
38.     if(First)  
39.     {  
40.         if(k<=0)  
41.         {  
42.             p->Next=First;  
43.             First=p;  
44.         }  
45.         if(k>Length-1)  
46.         {  
47.             Last->Next=p;   
48.             Last=Last->Next;  
49.             Last->Next=0;  
50.         }  
51.         if(k>0 && k<length)  
52.         {  
53.             k--;  
54.             TNode *q=First;  
55.             while(k-->0)  
56.                 q=q->Next;  
57.             p->Next=q->Next;  
58.             q->Next=p;  
59.         }  
60.     }  
61.     else  
62.     {  
63.         First=Last=p;  
64.         First->Next=Last->Next=0;  
65.     }  
66.     Length++;  
67.     return 1;  
68. }  
69.   
70. template Type TList::GetData(int k)  
71. {  
72.     TNode *p=First;  
73.     while(k-->0)  
74.         p=p->Next;  
75.     return p->Data;  
76. }  
77.   
78. template void TList::SetData(const Type& T,int k)  
79. {  
80.     TNode *p=First;  
81.     while(k-->0)  
82.         p=p->Next;  
83.     p->Data=T;  
84. }  

线性表单独以头文件形式存放。

四.广度优先搜索法
在搜索法中，广度优先搜索法是寻找最短路经的首选。
1.广度优先搜索算法的基本步骤
1）建立一个队列，将初始结点入队，并设置队列头和尾指针
2）取出队列头（头指针所指）的结点进行扩展，从它扩展出子结点，并将这些结点按扩展的顺序加入队列。 
3）如果扩展出的新结点与队列中的结点重复，则抛弃新结点，跳至第六步。
4）如果扩展出的新结点与队列中的结点不重复，则记录其父结点，并将它加入队列，更新队列尾指针。
5）如果扩展出的结点是目标结点，则输出路径，程序结束。否则继续下一步。
6）如果队列头的结点还可以扩展，直接返回第二步。否则将队列头指针指向下一结点，再返回第二步。
2.搜索路径的输出
搜索到目标结点后，需要输出搜索的路径。每个结点有一个数据域last，它记录了结点的父结点，因此输出搜索路径时，就是从目标结点Q出发，根据last找到它的父结点，再根据这个结点的last找到它的父结点，....,最后找到初始结点。搜索的路径就是从初始结点循相反方向到达目标结点的路径。
3.广度优先搜索法TBFS类的结构
广度优先搜索法TBFS类是作为TEight类的一个子类。其类的结构和成员函数的实现如下：

[cpp] view plaincopyprint?

1. class TBFS:public TEight  
2. {  
3. public:  
4.     TBFS(){}  
5.     TBFS(char *fname):TEight(fname){}  
6.     virtual void Search();  
7. private:  
8.     void Printl(TList &L);  
9.     int Repeat(TList &L);  
10.     int Find();  
11. };  
12.   
13. void TBFS::Printl(TList &L)  
14. {  
15.     TBFS T=*this;  
16.     if(T.last==-1)  
17.         return;  
18.     else  
19.     {  
20.         T=L.GetData(T.last);  
21.         T.Printl(L);  
22.         T.Printf();  
23.     }  
24. }  
25.   
26. int TBFS::Repeat(TList &L)  
27. {  
28.     int n=L.Getlen();  
29.     for(int i=0;i<n;i++)  
30.         if(L.GetData(i)==*this)  
31.             break;  
32.     return i;  
33. }  
34.   
35. int TBFS::Find()  
36. {  
37.     for(int i=0;i<num;)  
38.         if(p[i]!=q[i++])  
39.             return 0;  
40.     return 1;  
41. }  
42.   
43. void TBFS::Search()  
44. {  
45.     TBFS T=*this;  
46.     TList L;  
47.     L.Append(T);  
48.     int head=0,tail=0;  
49.     while(head<=tail)  
50.     {  
51.         for(int i=0;i<4;i++)  
52.         {  
53.             T=L.GetData(head);  
54.             if(T.Extend(i) && T.Repeat(L)>tail)  
55.             {  
56.                 T.last=head;  
57.                 L.Append(T);  
58.                 tail++;  
59.             }  
60.             if(T.Find())  
61.             {  
62.                 T.Printl(L);  
63.                 T.Printf();  
64.                 return;  
65.             }  
66.         }  
67.         head++;  
68.     }  
69. }  

4.广度优先搜索法的缺点
广度优先搜索法在有解的情形总能保证搜索到最短路经，也就是移动最少步数的路径。但广度优先搜索法的最大问题在于搜索的结点数量太多，因为在广度优先搜索法中，每一个可能扩展出的结点都是搜索的对象。随着结点在搜索树上的深度增大，搜索的结点数会很快增长，并以指数形式扩张，从而所需的存储空间和搜索花费的时间也会成倍增长。

五.双向广度优先搜索法
1.双向广度优先搜索法
八数码问题具有可逆性，也就是说，如果可以从一个状态A扩展出状态B，那么同样可以从状态B扩展出状态A，这种问题既可以从初始状态出发，搜索目标状态，也可以从目标状态出发，搜索初始状态。对这类问题如果采用双向广度优先搜索法，将可以大大节省搜索的时间。
所谓双向广度优先搜索法，是同时从初始状态和目标状态出发，采用广度优先搜索的策略，向对方搜索，如果问题存在解，则两个方向的搜索会在中途相遇，即搜索到同一个结点。将两个方向的搜索路径连接起来，就可以得到从初始结点到目标结点的搜索路径。
2.双向广度优先搜索算法
双向广度优先搜索算法的基本步骤如下：
1）建立两个队列，一个是正向搜索的队列，另一个是反向搜索的队列。将初始结点放入正向队列，将目标结点放入反向队列，并设置两个队列的头和尾指针。
2）从正向队列取出队列头（头指针所指）的结点进行扩展。
3）如果扩展出的新结点与队列中的结点重复，则抛弃新结点，跳至第六步。
4）如果扩展出的新结点与队列中的结点不重复，则记录其父结点，并将它加入队列，更新队列尾指针。
5）检查扩展出的结点是否在另一方向的队列中，如果是则两个方向的搜索相遇，显示搜索路径，程序结束。否则继续下一步。
6）如果队列头的结点还可以扩展，直接返回第二步。否则将队列头指针指向下一结点，然后对另一方向搜索的队列，按照第二步开始的同样步骤处理。
3.双向广度优先搜索法的优势
广度优先搜索法搜索时，结点不断扩张，深度越大，结点数越多。如果从两个方向向对方搜索，就会在路径中间某个地方相会，这样，双方的搜索的深度都不大，所搜索过的结点数就少得多，搜索时间也就节省不少。
从理论上说，如果每一结点可扩展的子结点数为m，广度优先搜索的搜索树就是一颗m叉树，也就是每个结点都由m个分支。按完全m叉树计算，如果目标结点在第n层，广度优先搜索就必须在搜索树上扩展完n-1层的所有结点，扩展的结点数为m(mn-1)/(m-1)。对于双向广度优先搜索来说，如果两个方向的搜索在第i层生成同一子结点，那么正向搜索扩展的结点数为m(mi-1)/(m-1)，反向搜索扩展的结点数为m(mn-i-1)/(m-1)，搜索的结点总数为m(mi+mn-i-1)/(m-1)（其中n是最优解路径长度，i=(m+1) div 2，）。设n为偶数(n=2*i)，广度优先双向搜索扩展的结点数约是广度优先搜索的2/(mi/2+1)*100%，相对减少(mi/2-1)/(mi/2+1)*100%。
4.判断两个方向的搜索相遇
在双向广度优先搜索法中，如何判断两个方向的搜索相遇呢？只要我们在生成结点的同时，判断该结点是否出现在相反方向的搜索树上即可，也就是说，在某个方向搜索中扩展出一个新结点，如果它与另一个方向已扩展出的结点重复，也就找到了解。
5.双向广度优先搜索法的TDBFS类结构
双向广度优先搜索法的TDBFS和广度优先搜索法类似，也是TEight类的子类，类结构及其成员函数的实现如下：

[cpp] view plaincopyprint?

1. class TDBFS:public TEight  
2. {  
3. public:  
4.     TDBFS(){}  
5.     TDBFS(char *fname):TEight(fname){}  
6.     virtual void Search();  
7. private:  
8.     void Printp(TList &L);  
9.     void Printb(TList &L);  
10.     int Repeat(TList &L);   
11. };  
12.   
13. void TDBFS::Printp(TList &L)  
14. {  
15.     TDBFS T=*this;  
16.     if(T.last==-1)  
17.         return;  
18.     else  
19.     {  
20.         T=L.GetData(T.last);  
21.         T.Printp(L);  
22.         T.Printf();  
23.     }  
24. }  
25.   
26. void TDBFS::Printb(TList &L)  
27. {  
28.     TDBFS T=*this;  
29.     while(T.last>-1)  
30.     {  
31.         T=L.GetData(T.last);  
32.         T.Printf();  
33.     }   
34. }  
35.   
36. int TDBFS::Repeat(TList &L)  
37. {  
38.     int n=L.Getlen();  
39.     for(int i=0;i<n;)  
40.         if(L.GetData(i++)==*this)  
41.             break;  
42.     return i;  
43. }  
44.   
45. void TDBFS::Search()  
46. {  
47.     TDBFS T1=*this;  
48.     TDBFS T2;  
49.     for(int i=0;i<num;i++)  
50.     {  
51.         T2.p[i]=q[i];  
52.         if(q[i]==0)  
53.             T2.spac=i;  
54.     }  
55.     T2.last=-1;  
56.     TList L1,L2;  
57.     L1.Append(T1);  
58.     L2.Append(T2);  
59.     int head1=0,tail1=0,head2=0,tail2=0;  
60.     while(head1<=tail1 || head2<=tail2)  
61.     {  
62.         for(int i=0;i<4;i++)  
63.         {  
64.             T1=L1.GetData(head1);  
65.             if(T1.Extend(i) && T1.Repeat(L1)>tail1)  
66.             {  
67.                 T1.last=head1;  
68.                 L1.Append(T1);  
69.                 tail1++;  
70.             }  
71.             int m=T1.Repeat(L2);  
72.             if(m<tail2)  
73.             {  
74.                 T1.Printp(L1);  
75.                 T1.Printf();  
76.                 T2=L2.GetData(m);  
77.                 T2.Printb(L2);  
78.                 return;  
79.             }  
80.         }  
81.         head1++;  
82.         for(i=0;i<4;i++)  
83.         {  
84.             T2=L2.GetData(head2);  
85.             if(T2.Extend(i) && T2.Repeat(L2)>tail2)  
86.             {  
87.                 T2.last=head2;  
88.                 L2.Append(T2);  
89.                 tail2++;  
90.             }  
91.             int m=T2.Repeat(L1);  
92.             if(m<tail1)  
93.             {  
94.                 T1=L1.GetData(m);  
95.                 T1.Printb(L1);  
96.                 T1.Printf();  
97.                 T2.Printp(L2);  
98.   
99.                 return;  
100.             }  
101.         }  
102.         head2++;  
103.     }  
104. }  

六.A*算法
1.启发式搜索
广度优先搜索和双向广度优先搜索都属于盲目搜索，这在状态空间不大的情况下是很合适的算法，可是当状态空间十分庞大时，它们的效率实在太低，往往都是在搜索了大量无关的状态结点后才碰到解答，甚至更本不能碰到解答。
搜索是一种试探性的查寻过程，为了减少搜索的盲目性引，增加试探的准确性，就要采用启发式搜索了。所谓启发式搜索就是在搜索中要对每一个搜索的位置进行评估，从中选择最好、可能容易到达目标的位置，再从这个位置向前进行搜索，这样就可以在搜索中省略大量无关的结点，提高了效率。
2.A*算法
A*算法是一种常用的启发式搜索算法。
在A*算法中，一个结点位置的好坏用估价函数来对它进行评估。A*算法的估价函数可表示为： 
f'(n) = g'(n) + h'(n) 
这里，f'(n)是估价函数，g'(n)是起点到终点的最短路径值（也称为最小耗费或最小代价），h'(n)是n到目标的最短路经的启发值。由于这个f'(n)其实是无法预先知道的，所以实际上使用的是下面的估价函数：
f(n) = g(n) + h(n) 
其中g(n)是从初始结点到节点n的实际代价，h(n)是从结点n到目标结点的最佳路径的估计代价。在这里主要是h(n)体现了搜索的启发信息，因为g(n)是已知的。用f(n)作为f'(n)的近似，也就是用g(n)代替g'(n)，h(n)代替h'(n)。这样必须满足两个条件：（1）g(n)>=g'(n)（大多数情况下都是满足的，可以不用考虑），且f必须保持单调递增。（2）h必须小于等于实际的从当前节点到达目标节点的最小耗费h(n)<=h'(n)。第二点特别的重要。可以证明应用这样的估价函数是可以找到最短路径的。
3.A*算法的步骤
A*算法基本上与广度优先算法相同，但是在扩展出一个结点后，要计算它的估价函数，并根据估价函数对待扩展的结点排序，从而保证每次扩展的结点都是估价函数最小的结点。
A*算法的步骤如下：
1）建立一个队列，计算初始结点的估价函数f，并将初始结点入队，设置队列头和尾指针。
2）取出队列头（队列头指针所指）的结点，如果该结点是目标结点，则输出路径，程序结束。否则对结点进行扩展。 
3）检查扩展出的新结点是否与队列中的结点重复，若与不能再扩展的结点重复（位于队列头指针之前），则将它抛弃；若新结点与待扩展的结点重复（位于队列头指针之后），则比较两个结点的估价函数中g的大小，保留较小g值的结点。跳至第五步。
4）如果扩展出的新结点与队列中的结点不重复，则按照它的估价函数f大小将它插入队列中的头结点后待扩展结点的适当位置，使它们按从小到大的顺序排列，最后更新队列尾指针。
5）如果队列头的结点还可以扩展，直接返回第二步。否则将队列头指针指向下一结点，再返回第二步。
4.八数码问题的A*算法的估价函数
估价函数中，主要是计算h，对于不同的问题，h有不同的含义。那么在八数码问题中，h的含意是各什么？八数码问题的一个状态实际上是数字0~8的一个排列，用一个数组p[9]来存储它，数组中每个元素的下标，就是该数在排列中的位置。例如，在一个状态中，p[3]=7，则数字7的位置是3。如果目标状态数字3的位置是8，那么数字7对目标状态的偏移距离就是3，因为它要移动3步才可以回到目标状态的位置。
八数码问题中，每个数字可以有9个不同的位置，因此，在任意状态中的每个数字和目标状态中同一数字的相对距离就有9*9种，可以先将这些相对距离算出来，用一个矩阵存储，这样只要知道两个状态中同一个数字的位置，就可查出它们的相对距离，也就是该数字的偏移距离：
0 1 2 3 4 5 6 7 8
0 0 1 2 1 2 3 2 3 4
1 1 0 1 2 1 2 3 2 3
2 2 1 0 3 2 1 4 3 2
3 1 2 3 0 1 2 1 2 3
4 2 1 2 1 0 1 2 1 2
5 3 2 1 2 1 0 3 2 1
6 2 3 4 1 2 3 0 1 2
7 3 2 3 2 1 2 1 0 1
8 4 3 2 3 2 1 2 1 0
例如在一个状态中，数字8的位置是3，在另一状态中位置是7，那么从矩阵的3行7列可找到2，它就是8在两个状态中的偏移距离。
估价函数中的h就是全体数字偏移距离之和。
显然，要计算两个不同状态中同一数字的偏移距离，需要知道该数字在每个状态中的位置，这就要对数组p[9]进行扫描。由于状态发生变化，个数字的位置也要变化，所以每次计算h都沿线扫描数组，以确定每个数字在数组中的位置。为了简化计算，这里用一个数组存储状态中各个数字的位置，并让它在状态改变时随着变化，这样就不必在每次计算h时，再去扫描状态数组。
例如，某个状态中，数字5的位置是8，如果用数组r[9]存储位置，那么就有r[5]=8。
现在用数组r[9]存储当前状态的数字位置，而用s[9]存储目标状态的数字位置，那么当前状态数字i对目标状态的偏移距离就是矩阵中r[i]行s[i]列对应的值。
5.A*算法的类结构
A*算法的类声明如下：

[cpp] view plaincopyprint?

1. class TAstar:public TEight  
2. {  
3. public:  
4.     TAstar(){} //构造函数  
5.     TAstar(char *fname); //带参数构造函数   
6.     virtual void Search(); //A*搜索法  
7. private:  
8.     int f,g,h; //估价函数   
9.     int r[Num]; //存储状态中各个数字位置的辅助数组  
10.     static int s[Num]; //存储目标状态中各个数字位置的辅助数组  
11.     static int e[]； //存储各个数字相对距离的辅助数组   
12.         void Printl(TList L); //成员函数，输出搜索路径   
13.     int Expend(int i); //成员函数，A*算法的状态扩展函数   
14.     int Calcuf(); //成员函数，计算估价函数  
15.     void Sort(TList& L,int k); //成员函数，将新扩展结点按f从小到大  
16.     //顺序插入待扩展结点队列  
17.     int Repeat(TList &L); //成员函数，检查结点是否重复  
18. };  
19.   
20. int TAstar::s[Num],TAstar::e[Num*Num];  
21.   
22. TAstar::TAstar(char *fname):TEight(fname)  
23. {   
24.     for(int i=0;i<num;)  
25.     {  
26.         r[p[i]]=i; //存储初始状态个个数字的位置  
27.         s[q[i]]=i++; //存储目标状态个个数字的位置  
28.     }  
29.     ifstream fin;   
30.     fin.open(".\Eightd.txt",ios::in | ios::nocreate);//打开数据文件  
31.     if(!fin)  
32.     {  
33.         cout<<"不能打开数据文件!"<<endl;  
34.         return;  
35.     }  
36.     for(i=0;i<num*num;i++) 读入各个数字相对距离值=""   
37.         fin>>e[i];  
38.     fin.close();  
39.     f=g=h=0; //估价函数初始值   
40. }  
41.   
42. void TAstar::Printl(TList L)   
43. {  
44.     TAstar T=*this;  
45.     if(T.last==-1)  
46.         return;  
47.     else  
48.     {  
49.         T=L.GetData(T.last);  
50.         T.Printl(L);  
51.         T.Printf();  
52.     }  
53. }  
54.   
55. int TAstar::Expend(int i)   
56. {  
57.     if(Extend(i)) //结点可扩展  
58.     {  
59.         int temp=r[p[r[0]]]; //改变状态后数字位置变化，存储改变后的位置   
60.         r[p[r[0]]]=r[0];  
61.         r[0]=temp;  
62.         return 1;  
63.     }  
64.     return 0;  
65. }  
66.   
67.   
68. int TAstar::Calcuf()  
69. {  
70.     h=0;  
71.     for(int i=0;i<num;i++) 计算估价函数的h  
72.         h+=e[Num*r[i]+s[i]];  
73.     return ++g+h;   
74. }  
75.   
76. void TAstar::Sort(TList& L,int k)  
77. {  
78.     int n=L.Getlen();  
79.     for(int i=k+1;i<n;i++)  
80.     {  
81.         TAstar T=L.GetData(i);  
82.         if(this->f<=T.f)  
83.             break;  
84.     }  
85.     L.Insert(*this,i);  
86. }  
87.   
88. int TAstar::Repeat(TList &L)  
89. {  
90.     int n=L.Getlen();  
91.     for(int i=0;i<n;i++)  
92.         if(L.GetData(i)==*this)  
93.             break;  
94.     return i;  
95. }  
96.   
97. void TAstar::Search()  
98. {  
99.     TAstar T=*this; //初始结点  
100.     T.f=T.Calcuf(); //初始结点的估价函数  
101.     TList L; //建立队列  
102.     L.Append(T); //初始结点入队  
103.     int head=0,tail=0; //队列头和尾指针  
104.     while(head<=tail) //队列不空则循环  
105.     {  
106.         for(int i=0;i<4;i++) //空格可能移动方向  
107.         {  
108.             T=L.GetData(head); //去队列头结点  
109.             if(T.h==0) //是目标结点  
110.             {   
111.                 T.Printl(L);//输出搜索路径  
112.                 T.Print(); //输出目标状态  
113.                 return; //结束  
114.             }  
115.             if(T.Expend(i)) //若结点可扩展  
116.             {  
117.                 int k=T.Repeat(L); //返回与已扩展结点重复的序号 if(k<head) 如果是不能扩展的结点  
118.                 continue; //丢弃  
119.                 T.last=head; //不是不能扩展的结点，记录父结点  
120.                 T.f=T.Calcuf(); //计算f  
121.                 if(k<=tail) //新结点与可扩展结点重复  
122.                 {  
123.                     TAstar Temp=L.GetData(k);   
124.                     if(Temp.g>T.g) //比较两结点g值  
125.                         L.SetData(T,k); //保留g值小的  
126.                     continue;  
127.                 }   
128.                 T.Sort(L,head) ; //新结点插入可扩展结点队列 tail++; //队列尾指针后移   
129.             }  
130.         }  
131.         head++; //一个结点不能再扩展，队列头指针指向下一结点  
132.     }  
133. }  

七.测试程序
八.算法运行结果
1.BFS算法只能适用于到达目标结点步数较少的情况，如果步数超过15步，运行时间太长，实际上不再起作用。
2.对于随机生成的同一个可解状态，BFS算法最慢，DBFS算法较慢，A*算法较快。但在15步以内，DBFS算法与A*算法相差时间不大，超过15步后，随步数增加，A*算法的优势就逐渐明显，A*算法要比DBFS算法快5倍以上，并随步数增大而增大。到25步以上，DBFS同样因运行时间过长而失去价值。
3.一般来说，解答的移动步数每增加1，程序运行时间就要增加5倍以上。由于八数码问题本身的特点，需要检查的节点随步数增大呈指数形式增加，即使用A*算法，也难解决移动步数更多的问题。
九.问题可解性
八数码问题的一个状态实际上是0~9的一个排列，对于任意给定的初始状态和目标，不一定有解，也就是说从初始状态不一定能到达目标状态。因为排列有奇排列和偶排列两类，从奇排列不能转化成偶排列或相反。
如果一个数字0~8的随机排列871526340，用F(X)表示数字X前面比它小的数的个数，全部数字的F(X)之和为Y=∑(F(X))，如果Y为奇数则称原数字的排列是奇排列，如果Y为偶数则称原数字的排列是偶排列。
例如871526340这个排列的
Y=0+0+0+1+1+3+2+3+0=10
10是偶数，所以他偶排列。871625340
Y=0+0+0+1+1+2+2+3+0=9
9是奇数，所以他奇排列。
因此，可以在运行程序前检查初始状态和目标状态的窘是否相同，相同则问题可解，应当能搜索到路径。否则无解。

[cpp] view plaincopyprint?

1. #include   
2. #include   
3. /* \Author:Kenny wong  
4. * \Email:[e-mail]  
5. * \version:0.1  
6. * \Date:2006-11-26  
7. * \bref:八数码问题A*算法解法  
8. */   
9.   
10. class node   
11. {   
12. public:   
13.     int data[9];   
14.     int zeroPos;   
15.     int distanceToOri;   
16.     int distanceToDis;   
17.     node *parent;   
18. public:   
19.     node(){}   
20.     node(int *data, int zeroPos, node *parent, int disToOri, int disToDis)   
21.         :zeroPos(zeroPos),parent(parent),distanceToOri(disToOri),   
22.         distanceToDis(disToDis)   
23.     {   
24.         memcpy(this -> data, data,sizeof(int)*9);   
25.     }   
26.     void setParent(node *parent)   
27.     {   
28.         this -> parent = parent;   
29.     }   
30.     int getRight()   
31.     {   
32.         return distanceToDis + distanceToOri;   
33.     }   
34. };   
35. class AStar   
36. {   
37. private:   
38.     std::vector openList;   
39.     std::vector closeList;   
40.     int grid[9];//当前棋盘   
41.     int target[9];//目标棋盘   
42.     int zeroPos;   
43.     int finZeroPos;   
44.     typedef void (AStar::*p)(int*,int);   
45.     p direction[4];   
46. public:   
47.     AStar(int *grid,int *target,int pos,int finPos)   
48.     {   
49.         memcpy(this -> grid, grid, sizeof(int)*9);   
50.         memcpy(this -> target, target, sizeof(int)*9);   
51.         zeroPos = pos;   
52.         finZeroPos = finPos;   
53.         direction[0] = &AStar::up;   
54.         direction[1] = &AStar::down;   
55.         direction[2] = &AStar::left;   
56.         direction[3] = &AStar::right;   
57.     }   
58.     ~AStar(){}   
59.     void swap(int *a, int *b)   
60.     {   
61.         int temp = *a;   
62.         *a = *b;   
63.         *b = temp;   
64.     }   
65.     void up(int *input,int pos)   
66.     {   
67.         memcpy(grid,input,sizeof(int)*9);   
68.         zeroPos = pos;   
69.         swap(&grid[zeroPos], &grid[zeroPos - 3]);   
70.         zeroPos -= 3;   
71.     }   
72.     void down(int *input,int pos)   
73.     {   
74.         memcpy(grid,input,sizeof(int)*9);   
75.         zeroPos = pos;   
76.         swap(&grid[zeroPos], &grid[zeroPos + 3]);   
77.         zeroPos += 3;   
78.     }   
79.     void left(int *input,int pos)   
80.     {   
81.         memcpy(grid,input,sizeof(int)*9);   
82.         zeroPos = pos;   
83.         swap(&grid[zeroPos], &grid[zeroPos - 1]);   
84.         zeroPos -= 1;   
85.     }   
86.     void right(int *input,int pos)   
87.     {   
88.         memcpy(grid,input,sizeof(int)*9);   
89.         zeroPos = pos;   
90.         swap(&grid[zeroPos], &grid[zeroPos + 1]);   
91.         zeroPos += 1;   
92.     }   
93.     bool check(int i,int Pos)   
94.     {   
95.         if(Pos - 3 < 0)//第一行   
96.         {   
97.             switch(Pos)   
98.             {   
99.             case 0:   
100.                 if( i == 0 || i == 2)   
101.                     return false;   
102.                 break;   
103.             case 1:   
104.                 if( i == 0)   
105.                     return false;   
106.                 break;   
107.             case 2:   
108.                 if( i == 0 || i == 3)   
109.                     return false;   
110.                 break;   
111.             default:   
112.                 break;   
113.             }   
114.         }   
115.         if(Pos - 3 >= 3)//第三行   
116.         {   
117.             switch(Pos)   
118.             {   
119.             case 6:   
120.                 if( i == 1 || i == 2)   
121.                     return false;   
122.                 break;   
123.             case 7:   
124.                 if( i == 1)   
125.                     return false;   
126.                 break;   
127.             case 8:   
128.                 if( i == 1 || i == 3)   
129.                     return false;   
130.                 break;   
131.             default:   
132.                 break;   
133.             }   
134.         }   
135.         else   
136.         {   
137.             switch(Pos)   
138.             {   
139.             case 3:   
140.                 if( i == 2 )   
141.                     return false;   
142.                 break;   
143.             case 5:   
144.                 if( i == 3 )   
145.                     return false;   
146.                 break;   
147.             default:   
148.                 break;   
149.             }   
150.         }   
151.         return true;   
152.     }   
153.     bool complete()   
154.     {   
155.     }   
156.     node *findNext()//寻找下一个待处理结点   
157.     {   
158.         std::vector::iterator it = openList.begin();   
159.         node *result = NULL;   
160.         //printf("%d\n",openList.size());   
161.         if(openList.size() == 1)   
162.         {   
163.             result = *it;   
164.             addToClose(result);   
165.             openList.clear();   
166.         }   
167.         if(openList.size() > 1)   
168.         {   
169.             ++it;   
170.             std::vector::iterator end = openList.end();   
171.             std::vector::iterator minNode = openList.begin();   
172.             int minRight = (*minNode) -> getRight();   
173.             for( ; it != end; ++it)   
174.             {   
175.                 int cur = (*it) -> getRight();   
176.                 if(cur < minRight)   
177.                 {   
178.                     minRight = cur;   
179.                     minNode = it;   
180.                 }   
181.             }   
182.             result = *minNode;   
183.             addToClose(result);   
184.             openList.erase(minNode);   
185.         }   
186.         return result;   
187.     }   
188.     void addToClose(node *n)   
189.     {   
190.         closeList.push_back(n);   
191.     }   
192.     void addToOpen(node *n)   
193.     {   
194.         openList.push_back(n);   
195.     }   
196.     void processCur(node *cur)//处理当前结点   
197.     {   
198.         for(int i = 0 ; i < 4; ++i)   
199.         {   
200.             //检测是否可达   
201.             if(check(i,cur -> zeroPos))   
202.             {   
203.                 (this ->* direction[i])(cur -> data,cur -> zeroPos);   
204.                 node *tem1 = createNode();   
205.   
206.                 //检测是否在closeList中   
207.                 if(!inCloseList(tem1))   
208.                 {   
209.                     node *tem2 = NULL;   
210.                     //检测是否在openList中   
211.                     if( (tem2 = findinOpenList(tem1)) == NULL)   
212.                     {   
213.                         tem1 -> distanceToOri = cur -> distanceToOri + calDistance();   
214.                         tem1 -> distanceToDis = calDistanceToDis();   
215.                         tem1 -> parent = cur;   
216.                         openList.push_back(tem1);   
217.                     }   
218.                     else   
219.                     {   
220.                         int temToOri = cur -> distanceToOri + calDistance();//计算经过当前点到达那个点的距离   
221.                         if(temToOri < tem2 -> distanceToOri)//如果经过当前点到那个点更好   
222.                         {   
223.                             tem2 -> parent = cur;   
224.                             tem2 -> distanceToOri = temToOri;   
225.                             tem2 -> distanceToDis = calDistanceToDis();   
226.                         }   
227.                     }   
228.                 }   
229.             }   
230.             else   
231.             {   
232.                 continue;   
233.             }   
234.         }   
235.     }   
236.     int calDistanceToDis()   
237.     {   
238.         int result = 0;   
239.         for(int i = 0; i < 9; ++i)   
240.             if(grid[i] == target[i])   
241.                 ++result;   
242.         return result;   
243.     }   
244.     int calDistance()   
245.     {   
246.         return 1;   
247.     }   
248.     node *createNode()   
249.     {   
250.         return new node(grid,zeroPos,NULL,0,0);   
251.     }   
252.     bool inCloseList(node *_node)   
253.     {   
254.         std::vector::iterator it = closeList.begin();   
255.         std::vector::iterator end = closeList.end();   
256.         for( ; it != end; ++it)   
257.             if(memcmp(_node -> data, (*it) -> data,sizeof(int)*9) == 0 && _node -> zeroPos == (*it) -> zeroPos)   
258.                 return true;   
259.         return false;   
260.     }   
261.     node *findinOpenList(node *_node)   
262.     {   
263.         std::vector::iterator it = openList.begin();   
264.         std::vector::iterator end = openList.end();   
265.         for( ; it != end; ++it)   
266.             if(memcmp(_node -> data, (*it) -> data,sizeof(int)*9) == 0 && _node -> zeroPos == (*it) -> zeroPos)   
267.                 return *it;   
268.         return NULL;   
269.     }   
270.     void begin()   
271.     {   
272.         node *ori = new node(grid,zeroPos,NULL,0,0);   
273.         addToOpen(ori);   
274.         node *distination = NULL;   
275.         node *finish = new node(target,finZeroPos,NULL,0,0);   
276.         int first[9];   
277.         memcpy(first,grid,sizeof(int)*9);   
278.         while(openList.size() > 0)   
279.         {   
280.             if( (distination = findinOpenList(finish)) != NULL)   
281.             {   
282.                 break;   
283.             }   
284.             node *temp = findNext();   
285.             if(temp != NULL)   
286.                 processCur(temp);   
287.         }   
288.   
289.         std::vector path;   
290.         if(distination != NULL && distination -> parent != NULL)   
291.         {   
292.             node *temp = distination;   
293.             while(temp -> parent != NULL)   
294.             {   
295.                 path.push_back(temp -> data);   
296.                 temp = temp -> parent;   
297.             }   
298.         }   
299.         path.push_back(first);   
300.         std::vector::iterator it = path.end();   
301.         --it;   
302.         std::vector::iterator begin = path.begin();   
303.         for( ; it != begin; --it)   
304.         {   
305.             for(int i = 0; i < 9; ++i)   
306.             {   
307.                 printf(" %d",(*it)[i]);   
308.                 if( (i + 1) % 3 == 0)   
309.                     printf("\n");   
310.             }   
311.             printf("\n");   
312.         }   
313.         for(int i = 0; i < 9; ++i)   
314.         {   
315.             printf(" %d",(*it)[i]);   
316.             if( (i + 1) % 3 == 0)   
317.                 printf("\n");   
318.         }   
319.     }   
320. };   
321. int main()   
322. {   
323.     int ori[9] = {1,2,5,7,6,4,0,8,3};   
324.     int tar[9] = {7,1,5,8,2,4,6,0,3};   
325.     AStar A(ori,tar,6,7);   
326.     A.begin(); return 0;   
327. }