人工智能课程的作业,正好把之前写过的八数码问题重新再学习,用C语言写一遍:
可以看看这篇文章:
一、使用的方法:暴力BFS+康托(逆)展开+hash判重
直接暴力bfs所有的询问,用康托展开计算全排列数值从而得到哈希值判重。
C语言学艺不精,没敢用指针和结构体,只用了简单的数组来实现链表,hash,队列等逻辑....QAQ....
虽然只用了数组,不过速度还是蛮快的,在我的电脑上基本上是秒出结果~。
#include<stdio.h>
int BFS(int origin[],int n,int target []);
int KT( int s[]);
int InsertHash(int h,int sum);
int Expand(int n,int Target);
//主函数,输入初始状态和目标状态
int main()
{
int origin[9]={1,2,3,4,5,6,7,8,0},target[9]={0,2,3,1,5,6,4,7,8};
/* printf("请输入初始状态:例如:\n1 2 3\n4 5 6\n7 8 0\n------*请输入:*------\n");
for(int i=0;i<9;i++)
{
scanf("%d",&origin[i]);
}
printf("请输入目标状态:例如:\n1 2 3\n4 5 6\n7 8 0\n------*请输入:*------\n");
for(i=0;i<9;i++)
{
scanf("%d",&target[i]);
}
*/
BFS(origin,8,target);
return 0;
}
//BFS函数,接收origin和target并开始广度优先搜索,使用hash判重.
//用Hash实现判重
const int MAXHASHSIZE=1200;//362880
int queue[MAXHASHSIZE],hash[MAXHASHSIZE],loc0[MAXHASHSIZE];//loc0[]记录0的位置
int num,num2;//num记录队尾,num2记录队头
//0所在每个位置的可扩展的位置
int loc[9][5];
int Ceng[MAXHASHSIZE],C,C1;//记录层数
int BFS(int origin[],int n,int target [])//n记录0的初始位置
{
//初始化
int h=0,sum=0; num=num2=0;C=C1=0;Ceng[C]=1;
int Target=0;
Target=KT(target);
//初始化0所在每个位置的可扩展的位置
loc[0][0]=1;loc[1][0]=0;loc[2][0]=1;loc[3][0]=0;loc[4][0]=1;loc[5][0]=2;loc[6][0]=3;loc[7][0]=4;loc[8][0]=5;
loc[0][1]=3;loc[1][1]=2;loc[2][1]=5;loc[3][1]=4;loc[4][1]=3;loc[5][1]=4;loc[6][1]=7;loc[7][1]=6;loc[8][1]=7;
loc[1][2]=4; loc[3][2]=6;loc[4][2]=5;loc[5][2]=8; loc[7][2]=8;
loc[4][3]=7;
//首元素放入队列
sum=KT(origin);
h=sum%MAXHASHSIZE;//得到哈希值
queue[num]=sum;
InsertHash(h,sum);
printf("%d",n);
printf("\n");
C1=C+1;
//广度搜索
loc0[num]=n;
Expand(n,Target);
if(InsertHash(h,sum))
{
num++;
queue[num]=sum;
}
return 0;
}
int Expand(int n,int Target)
{
///逆康托展开
int fac[] = {1,1,2,6,24,120,720,5040,40320};
int i, j, t, vst[8]={0},s[9];
int k=queue[num2];
k--;
for (i=0; i<9; i++)
{
t = k/fac[9-i-1];
k %= fac[9-i-1];
for (j=0; j<9; j++)
if (!vst[j])
{
if(t==0)break;
t--;
}
s[i] = j;
vst[j] = 1;
}
//显示最终结果
if(queue[num2]==Target)
{
printf("\n");
for(i=0;i<9;i++)
{
printf("%d",s[i]);
if(i%3==2)
printf("\n");
}
return 1;
}
int temp[9];
//队头扩展
int next=0;
while(loc[n][next]!=NULL || (n==1&&next==0&&loc[n][next]==0) || (n==3&&next==0&&loc[n][next]==0))
{
printf("%d ",loc[n][next]);
//广度扩展
for(i=0;i<9;i++)temp[i]=s[i];
temp[n]=temp[loc[n][next]];
temp[loc[n][next]]=0;
//判断是否已存在
int sum=KT(temp);
int h=sum%MAXHASHSIZE;//得到哈希值
if(InsertHash(h,sum))
{
//插入队列
num++;
queue[num]=sum;
loc0[num]=loc[n][next];
Ceng[C1]++;
}
next++;
}
printf("-");
Ceng[C]--;
if(Ceng[C]<=0)
{
printf("\n");
C++;
C1=C+1;
}
num2++;//读取下一个
n=loc0[num2];
if(Expand(n,Target))
return 1;
if(num2==362879)
return 0;
}
/* 康托函数 */
int KT( int s[])
{
int fac[] = {1,1,2,6,24,120,720,5040,40320};
int i, j, t, sum;
sum = 0;
for (i=0; i<9; i++)
{
t = 0;
for (j=i+1; j<9; j++)
if (s[j] < s[i])
t++;
sum += t*fac[9-i-1];
}
return sum+1;
}
//插入hash表函数
int InsertHash(int h,int sum)
{
if(hash[h]==sum)
return 0;
else
{
if(hash[h]==NULL)
{
hash[h]=sum;
return 1;
}
h++;
if(InsertHash(h,sum))
return 1;
else
return 0;
}
}