人工智能经典问题,八数码问题求解,多种搜索算法大全,BFS,DFS,UCS,A*,贪婪算法,迭代加深IDS,有界深度搜索,C语言版,保证看懂,分析到位,注释详细,没有bug
人工智能经典问题,八数码问题求解,多种搜索算法大全,BFS,DFS,UCS,A*,贪婪算法,迭代加深IDS,有界深度搜索,C语言版,保证看懂,分析到位,注释详细,没有bug
·
目录
废话不多说,直接上代码
想看详细分析的可以在我的博客里找,每个算法有对应的分析。
1、BFS(广度优先搜索法)
#include <stdio.h> #include <stdlib.h> #include <time.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 回收内存 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 tail_insert(open, p); // BFS用尾插法,DFS用头插法 } // 广度优先搜索法 void BFS() { printf("start BFS\n"); node *tmp; int x, y, k, cmp[N][N]; add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return ; } // 从open表中取出首元节点,插入到close表的末尾 tmp = open->next; Remove(tmp); tail_insert(close, tmp); printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); showWin(); return; } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); BFS(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("BFS finish\n"); return 0; }
2、DFS(深度优先搜索法)版本一
#include <stdio.h> #include <stdlib.h> #include <time.h> #define N 3 // 阶数 #define Up 1 // 上移 #define Down -1 // 下移 #define Left 2 // 左移 #define Right -2 // 右移 // 定义一个结构体来表示棋盘状态 typedef struct node { int x; // 父节点的x坐标 int y; // 父节点的y坐标 int op; // 记录空格操作 struct node *prev; // 链表中的前指针 struct node *next; // 链表中后指针 }node; node *open; node *close; static int n = 0; int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->x = -1; p->y = -1; p->op = 0; p->prev = p; p->next = p; return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 检查是否到达目的状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到不相同的立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印成功搜索出的路径 void showWin() { node *p; int i, j, k, tmp[N][N]; k = 1; copy(src, tmp); printf("the solution as follow:\n"); // 从初始状态src经过路径上是空格(上下左右)操作,一步一步打印出路径 for(p = close->next; p != close; p = p->next, k++) { switch(p->op) { case Up: up(tmp, p->x, p->y); break; case Down: down(tmp, p->x, p->y); break; case Left: left(tmp, p->x, p->y); break; case Right: right(tmp, p->x, p->y); break; } printf("step%d:\n", k); display(tmp); printf("\n"); } } // 检查是否有重复的状态 int checkSame(int cmp[N][N]) { node *p; int i, j, tmp[N][N]; copy(src, tmp); // 复制当前状态 for(p = close->next; p != close; p = p->next) { switch(p->op) { case Up: up(tmp, p->x, p->y); break; case Down: down(tmp, p->x, p->y); break; case Left: left(tmp, p->x, p->y); break; case Right: right(tmp, p->x, p->y); break; } // 检测从初始状态沿路径节点的操作推出的状态与待插入到open表的状态比较 // 如果所有的路径状态点与当前状态不相同就可以插入到open表中 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 如果有不相同的位置点,则修改i值 // 然后直接退出本次比较,进行下一个状态比较 if(tmp[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 读取当前状态 void getCurrent(node *p) { if(!p || p->next == p) { return; } // 根据当前路径节点的空格操作来推出当前状态 switch(p->op) { case Up: up(cur, p->x, p->y); break; case Down: down(cur, p->x, p->y); break; case Left: left(cur, p->x, p->y); break; case Right: right(cur, p->x, p->y); break; } } // 删除close表的末尾的节点,同时把cur回溯到上一个状态 void getBack(node *tmp, int x, int y) { // 循环删除close表末尾(tmp)中不符合的节点 // 因为可能存在tmp的父节点只有一个后继节点,即tmp // 此时tmp的父节点也需要删除 while(1) { tmp = close->prev; // tmp指向close表末尾 tmp->op = -tmp->op; // 空格操作取反 tmp->x = x; // 记录空格x坐标 tmp->y = y; // 记录空格y坐标 getCurrent(tmp); // 修改当前状态 Remove(tmp); // 弹出tmp节点 free(tmp); // 释放tmp节点 findSpace(cur, &x, &y); // 寻找当前节点的空格坐标 // 如果当前状态cur是由close表中的最末节点推算来的 // 且open表的首元节点,即下一个将插入close表的节点 // 那么close->prev就是open->next的父节点 // 即open->next中的坐标(x, y)等于cur的空格坐标 if(x == open->next->x && y == open->next->y) { break; } } } // 把当前节点(tmp) 的后继节点插入到open表中 void add(int x, int y, int op) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } p->op = op; // 修改空格操作 p->x = x; // 修改空格x坐标 p->y = y; // 修改空格y坐标 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 深度优先搜索法 void DFS() { printf("start DFS\n"); int x, y, k, cmp[N][N]; node *p, *tmp; copy(src, cur); findSpace(cur, &x, &y); add(x, y, 0); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return ; } // 从open表中取出首元节点,插入到close表的末尾 tmp = open->next; Remove(tmp); tail_insert(close, tmp); printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 getCurrent(tmp); printf("current->op: %d\n", tmp->op); printf("current:\n"); display(cur); if(checkWin()) { printf("success\n\n"); showWin(); return; } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(x, y, Up); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(x, y, Down); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(x, y, Left); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(x, y, Right); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,则应该从close表中删除该节点 if(k == 0) { // 删除节点的同时把cur回溯到上一个状态 getBack(tmp, x, y); printf("del done\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); DFS(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
3、DFS(深度优先搜索法)版本二
#include <stdio.h> #include <stdlib.h> #include <time.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 回收内存 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 深度优先搜索法 void DFS() { printf("start BFS\n"); node *tmp; int x, y, k, cmp[N][N]; add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return ; } // 从open表中取出首元节点,插入到close表的末尾 tmp = open->next; Remove(tmp); tail_insert(close, tmp); printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); printf("solution as follow:\n"); showWin(); return; } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); DFS(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
4、IDS(有界深度优先搜索法)
#include <stdio.h> #include <stdlib.h> #include <time.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 内存回收 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 重新搜索 void restart(int *depth) { node *p, *tmp; tmp = close->prev; // 使tmp指向close表末尾 Remove(tmp); // 弹出tmp节点 free(tmp); // 释放tmp节点 (*depth)--; // 搜索深度减一 // 如果下一个在open表中的第一个待搜索节点(open->next) // 其父节点(open->next->father)刚好是close表中的末尾节点 // 即close表中的末尾节点是删除了一个超过深度的节点(tmp)后节点 // 那么可以断定open->next与tmp是兄弟节点 if(open->next->father == close->prev) { return; } // 如果删除了close表末尾的一个节点后,close->prev与open->next不是兄弟节点 // 那么说明刚删除了tmp后,close表中的末尾节点只有一个后继节点 // close->prev只有一个后继节点tmp,则需要循环删除close的末尾节点 // 直到open->next与close->prev是子父节点为止,同时搜索深度不断减一 // 这里采用尾删法,与尾插法相似 for(tmp = close->prev; tmp != close; tmp = close->prev) { // 找到父节点,立马返回,进行下次搜索 if(open->next->father == tmp) { return; } Remove(tmp); // 弹出tmp节点 free(tmp); // 释放tmp节点 (*depth)--; // 搜索深度减一 } } // 有界深度优先搜索法 void DFS() { printf("start DFS\n"); node *tmp; int x, y, k, cmp[N][N]; int depth = 0; // 计算搜索的深度 int boundary = 5; // 设定深度不能超过5,这里解的深度为5 add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return ; } // 从open表中取出首元节点,插入到close表的末尾 tmp = open->next; Remove(tmp); tail_insert(close, tmp); depth++; // 深度加一 printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); printf("solution as follow:\n"); showWin(); return; } // 如果超过最深限度,则重新搜索 if(depth >= boundary) { printf("cross the boundary\n"); printf("restart\n\n"); restart(&depth); // 修改close表和open表 continue; // 重新搜索 } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 depth--; // 搜索深度减一 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); DFS(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
5、迭代加深优先搜索法
#include <stdio.h> #include <stdlib.h> #include <time.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 内存回收 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 重新搜索 void restart(int *depth) { node *p, *tmp; tmp = close->prev; // 使tmp指向close表末尾 Remove(tmp); // 弹出tmp节点 free(tmp); // 释放tmp节点 (*depth)--; // 搜索深度减一 // 如果下一个在open表中的第一个待搜索节点(open->next) // 其父节点(open->next->father)刚好是close表中的末尾节点 // 即close表中的末尾节点是删除了一个超过深度的节点(tmp)后节点 // 那么可以断定open->next与tmp是兄弟节点 if(open->next->father == close->prev) { return; } // 如果删除了close表末尾的一个节点后,close->prev与open->next不是兄弟节点 // 那么说明刚删除了tmp后,close表中的末尾节点只有一个后继节点 // close->prev只有一个后继节点tmp,则需要循环删除close的末尾节点 // 直到open->next与close->prev是子父节点为止,同时搜索深度不断减一 // 这里采用尾删法,与尾插法相似 for(tmp = close->prev; tmp != close; tmp = close->prev) { // 找到父节点,立马返回,进行下次搜索 if(open->next->father == tmp) { return; } Remove(tmp); // 弹出tmp节点 free(tmp); // 释放tmp节点 (*depth)--; // 搜索深度减一 } } // 迭代加深搜索法 void DFS() { printf("start BFS\n"); node *tmp; int x, y, k, cmp[N][N]; int depth = 0; // 计算搜索的深度 int boundary = 1; // 设定初始深度不能超过100,这里解的深度为405 add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); printf("the boundary is %d\n", boundary); printf("improve boundary\n"); boundary += 1; // 迭代加深,深度加1 clearList(close); // 清空close表 add(NULL, src); // 把初始状态加入到open表中 depth = 0; // 重新定义搜索深度 sleep(3); if(boundary > 1000) // 如果深度太深就直接结束 { printf("boundary can not improve any more\n"); return; } } // 从open表中取出首元节点,插入到close表的末尾 tmp = open->next; Remove(tmp); tail_insert(close, tmp); depth++; // 深度加一 printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); printf("solution as follow:\n"); showWin(); printf("boundary is %d\n", boundary); return; } // 如果超过最深限度,则重新搜索 if(depth >= boundary) { printf("cross the boundary\n"); printf("restart\n\n"); restart(&depth); // 修改close表和open表 continue; // 重新搜索 } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 depth--; // 搜索深度减一 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); DFS(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
6、UCS(等代价搜索法)
// 备注:采用曼哈顿距离计算初始状态与当前状态的水平距离和垂直距离之和 #include <stdio.h> #include <stdlib.h> #include <time.h> #include <math.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int cost; // 从初始状态到本状态的代价 int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; p->cost = 0; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 采用曼哈顿距离,计算从初始状态到本状态的代价 // 曼哈顿距离指,与目标位置的水平距离和垂直距离之和 // 书上给出的代价是从起始节点S到任一节点i的路径代价g(i) // 修改代价时,采用g(j) = g(i) + cost(i, j); // cost(i, j)指从当前i节点到它的后继节点j的代价 // 但是这里因为每次空格操作,只会改变2个位置 // 所以cost(i, j)都相同,因此可以省略,直接计算g(j) void compute_cost(node *p) { int i, j, m, n; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 如果有对应位置不匹配的 // 则重头搜索出与src[i][j]匹配的p->data[m][n] if(p->data[i][j] != src[i][j]) { for(m = 0; m < N; m++) { for(n = 0; n < N; n++) { // 搜索到与src[i][j]匹配的p->data[m][n] // 计算曼哈顿距离后立马结束搜索 if(p->data[m][n] == src[i][j]) { // 累加每个位置的曼哈顿距离 p->cost += abs(m-i) + abs(n-j); m = N; // 修改m值,结束搜索 break; } } } } } } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 内存回收 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 compute_cost(p); // 修改从初始状态到本状态的代价 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 从表中选择一个节点,使其代价p->cost最小 // 通常是从open表中选择,如果有多个点,则随便选一个 node *min_cost(node *head) { // 空表直接退出 if(!head || head->next == head) { return NULL; } node *p, *min; min = head->next; for(p = min->next; p != head; p = p->next) { // 遇到更小的,就更新min if(p->cost < min->cost) { min = p; } } return min; } // 等代价搜索法 void UCS() { printf("start UCS\n"); node *tmp; int x, y, k, cmp[N][N]; add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return; } // 从open表中选择一个节点,使其代价p->cost最小 tmp = min_cost(open); // 如果tmp为空,则说明open表为空 // 但是一般不会出现,因为前面已经对open表进行判空处理 // 这里写的原因就是特意为了增强算法的健壮性 if(!tmp) { printf("fail\n"); return; } Remove(tmp); tail_insert(close, tmp); printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); printf("solution as follow:\n"); showWin(); return; } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); UCS(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
7、贪婪算法
// 备注:采用曼哈顿距离计算当前状态与目标状态的水平距离和垂直距离之和 #include <stdio.h> #include <stdlib.h> #include <time.h> #include <math.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int cost; // 从当前状态到目标状态的代价 int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; p->cost = 0; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 采用曼哈顿距离,计算从当前状态到目标状态的代价 // 曼哈顿距离指,与目标位置的水平距离和垂直距离之和 // 书上给出的代价是从起始节点S到任一节点i的路径代价g(i) // 修改代价时,采用g(j) = g(i) + cost(i, j); // cost(i, j)指从当前i节点到它的后继节点j的代价 // 但是这里因为每次空格操作,只会改变2个位置 // 所以cost(i, j)都相同,因此可以省略,直接计算g(j) void compute_cost(node *p) { int i, j, m, n; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 如果有对应位置不匹配的 // 则重头搜索出与dest[i][j]匹配的p->data[m][n] if(p->data[i][j] != dest[i][j]) { for(m = 0; m < N; m++) { for(n = 0; n < N; n++) { // 搜索到与dest[i][j]匹配的p->data[m][n] // 计算曼哈顿距离后立马结束搜索 if(p->data[m][n] == dest[i][j]) { // 累加每个位置的曼哈顿距离 p->cost += abs(m-i) + abs(n-j); m = N; // 修改m值,结束搜索 break; } } } } } } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 内存回收 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 compute_cost(p); // 修改从当前状态到目标状态的代价 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 从表中选择一个节点,使其代价p->cost最小 // 通常是从open表中选择,如果有多个点,则随便选一个 node *min_cost(node *head) { // 空表直接退出 if(!head || head->next == head) { return NULL; } node *p, *min; min = head->next; for(p = min->next; p != head; p = p->next) { if(p->cost < min->cost) { min = p; } } return min; } // 贪婪算法 void Greedy() { printf("start Greedy\n"); node *tmp; int x, y, k, cmp[N][N]; add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return; } // 从open表中选择一个节点,使其代价p->cost最小 tmp = min_cost(open); // 如果tmp为空,则说明open表为空 // 但是一般不会出现,因为前面已经对open表进行判空处理 // 这里写的原因就是特意为了增强算法的健壮性 if(!tmp) { printf("fail\n"); return; } Remove(tmp); tail_insert(close, tmp); printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); printf("solution as follow:\n"); showWin(); return; } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); Greedy(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
8、启发式搜索法(A*算法)
// 备注:曼哈顿距离等于当前状态与目标状态的水平距离和垂直距离之和 #include <stdio.h> #include <stdlib.h> #include <time.h> #include <math.h> #define N 3 // 阶数,可以改为更高阶 // 定义一个结构体来表示棋盘状态 typedef struct node { int cost; // 从初始状态到本状态的代价 int data[N][N]; // 存放棋盘状态 struct node *prev; // 链表中的前指针 struct node *next; // 链表中的后指针 struct node *father;// 搜索树的父节点 }node; node *open; // open表,存放未拓展的节点 node *close; // close表,存放已经拓展的节点 static int n = 0; // 用来记录总共搜索的次数 int src[N][N] = {2, 8, 3, 1, 0, 4, 7, 6, 5}; // 初始状态 int cur[N][N]; // 当前状态 int dest[N][N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; // 目标状态 // 生成一个节点 node *initList() { int i, j; node *p = malloc(sizeof(node)); if(!p) { printf("malloc fail\n"); return NULL; } p->prev = p; p->next = p; p->father = NULL; p->cost = 0; // 初始化棋盘状态 for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { p->data[i][j] = -1; } } return p; } // 头插法,把节点插入到链表头部 void head_insert(node *head, node *p) { p->next = head->next; p->prev = head; head->next->prev = p; head->next = p; } // 尾插法,把节点插入到链表尾部 void tail_insert(node *head, node *p) { p->prev = head->prev; head->prev->next = p; p->next = head; head->prev = p; } // 弹出节点 void Remove(node *p) { p->prev->next = p->next; p->next->prev = p->prev; p->next = p; p->prev = p; } // 清空链表 void clearList(node *head) { // 空表不需要再清空 if(!head || head->next == head) { printf("list is null\n"); return; } // 头删法,跟头插法相似,循环删除首元节点 for(node *p = head->next; p != head; p = head->next) { Remove(p); // 弹出节点p free(p); // 释放节点p的空间 } } // 初始化棋盘 void init() { int i, j, k; int a[N*N] = {1, 2, 3, 8, 0, 4, 7, 6, 5}; int flag[N*N] = {0}; // 用来标记0~8中哪个数字已经出现过 srand(time(0)); for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 尝试生成一个0~8间的数字 k = rand()%(N*N); while(1) { // 如果该数字未出现 if(flag[k] == 0) { flag[k] = 1; // 修改标志位 src[i][j] = a[k]; // 初始化对应src的位置 break; } // 重新生成数字 k = rand()%(N*N); } } } } // 复制棋盘状态 void copy(int source[N][N], int dest1[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { dest1[i][j] = source[i][j]; } } } // 打印棋盘状态 void display(int s[N][N]) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 把0当做空格输出 if(s[i][j] > 0) { printf("%d ", s[i][j]); } else if(s[i][j] == 0) { printf(" "); } } printf("\n"); } } // 寻找出棋盘中的空格键 void findSpace(int s[N][N], int *x, int *y) { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(s[i][j] == 0) // 找到空格 { *x = i; *y = j; } } } } // 空格上移 void up(int s[N][N], int x, int y) { // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { s[x][y] = s[x-1][y]; s[x-1][y] = 0; } } // 空格下移 void down(int s[N][N], int x, int y) { // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { s[x][y] = s[x+1][y]; s[x+1][y] = 0; } } // 空格左移 void left(int s[N][N], int x, int y) { // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { s[x][y] = s[x][y-1]; s[x][y-1] = 0; } } // 空格右移 void right(int s[N][N], int x, int y) { // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { s[x][y] = s[x][y+1]; s[x][y+1] = 0; } } // 采用曼哈顿距离,计算从某个状态cmp到当前状态的代价 // 曼哈顿距离指,与目标位置的水平距离和垂直距离之和 // 书上给出的代价是从起始节点S到任一节点i的路径代价g(i) // 修改代价时,采用g(j) = g(i) + cost(i, j); // cost(i, j)指从当前i节点到它的后继节点j的代价 // 但是这里因为每次空格操作,只会改变2个位置 // 所以cost(i, j)都相同,因此可以省略,直接计算g(j) void compute_cost(node *p, int cmp[N][N]) { int i, j, m, n; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 如果有对应位置不匹配的 // 则重头搜索出与cmp[i][j]匹配的p->data[m][n] if(p->data[i][j] != cmp[i][j]) { for(m = 0; m < N; m++) { for(n = 0; n < N; n++) { // 搜索到与cmp[i][j]匹配的p->data[m][n] // 计算曼哈顿距离后立马结束搜索 if(p->data[m][n] == cmp[i][j]) { // 累加每个位置的曼哈顿距离 p->cost += abs(m-i) + abs(n-j); m = N; // 修改m值,结束搜索 break; } } } } } } } // 检查close表中是否有重复的状态 int checkSame(int cmp[N][N]) { int i, j; node *p; for(p = close->next; p != close; p = p->next) { for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { if(p->data[i][j] != cmp[i][j]) { i = N + 1; break; } } } if(i == N && j == N) { return 1; } } return 0; } // 检查当前棋盘状态是否为目标状态 int checkWin() { int i, j; for(i = 0; i < N; i++) { for(j = 0; j < N; j++) { // 检测到有不同的就立马返回,节省时间 if(cur[i][j] != dest[i][j]) { return 0; } } } return 1; } // 打印从初始状态到目标状态的路径 void showWin() { int i = 1; node *p, *win; // 创建一个win表,用来存放路径 win = initList(); if(!win) { printf("win malloc fail\n"); return; } // 由于节点中只有father(父节点),而打印路径是从父节点开始的 // 所以必须先从close表选择出正确的路径,因为close中可能存放不是正确路径的节点 // 可以唯一确定的是 close->prev 一定是目的状态,由此循环往上寻找其父节点 for(p = close->prev; p && p != close; p = p->father) { Remove(p); // 从close表中弹出p节点 head_insert(win, p); // 把p节点用头插法插入到win表中 } printf("the solution as follow:\n"); // 经过头插法创建的win表,一定是按从初始状态到目标状态的排序的 for(p = win->next; p != win; p = p->next, i++) { printf("step %d:\n", i); display(p->data); printf("\n"); } clearList(win); // 清空win表 free(win); // 内存回收 } // 把当前节点(tmp) 的后继节点插入到open表中 void add(node *tmp, int cmp[N][N]) { node *p = initList(); if(!p) { printf("malloc fail\n"); return; } copy(cmp, p->data); // 复制棋盘状态 p->father = tmp; // 修改父节点指针 compute_cost(p, src); // 修改从初始状态到当前状态的代价 compute_cost(p, dest); // 修改从当前状态到目标状态的代价 head_insert(open, p); // BFS用尾插法,DFS用头插法 } // 从表中选择一个节点,使其代价p->cost最小 // 通常是从open表中选择,如果有多个点,则随便选一个 node *min_cost(node *head) { // 空表直接退出 if(!head || head->next == head) { return NULL; } node *p, *min; min = head->next; for(p = min->next; p != head; p = p->next) { if(p->cost < min->cost) { min = p; } } return min; } // 等代价搜索法 void Astar() { printf("start Astar\n"); node *tmp; int x, y, k, cmp[N][N]; add(NULL, src); // 把初始状态加入到open表中 while(1) { // 如果open表为空,就可以退出了 if(open->next == open) { printf("fail\n"); return; } // 从open表中选择一个节点,使其代价p->cost最小 tmp = min_cost(open); // 如果tmp为空,则说明open表为空 // 但是一般不会出现,因为前面已经对open表进行判空处理 // 这里写的原因就是特意为了增强算法的健壮性 if(!tmp) { printf("fail\n"); return; } Remove(tmp); tail_insert(close, tmp); printf("now go :%d\n", ++n); // 打印当前搜索出的棋盘状态 copy(tmp->data, cur); printf("current:\n"); display(cur); // 检测当前棋盘状态是否为目标状态 if(checkWin()) { printf("success\n\n"); printf("solution as follow:\n"); showWin(); return; } k = 0; // 记录当前节点的后继节点个数 findSpace(cur, &x, &y); // 寻找当前棋盘状态的空格坐标 printf("try up\n"); // 空格必须是处于第2行或者第3行,才能上移 if(x >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 up(cmp, x, y); // cmp尝试向上移 // 判断cmp上移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can up\n"); } } printf("try down\n"); // 空格必须是处于第1行或者第2行,才能下移 if(x <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 down(cmp, x, y); // cmp尝试向下移 // 判断cmp下移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can down\n"); } } printf("try left\n"); // 空格必须是处于第2列或者第3列,才能左移 if(y >= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 left(cmp, x, y); // cmp尝试向左移 // 判断cmp左移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can left\n"); } } printf("try right\n"); // 空格必须是处于第1列或者第2列,才能右移 if(y <= 1) { copy(cur, cmp); // 先复制当前状态cur到cmp中 right(cmp, x, y); // cmp尝试向右移 // 判断cmp右移后close表中是否有重复状态 // 没有重复才可以把cmp状态加入到close表中 if(checkSame(cmp) == 0) { k++; // 当前节点的后继节点数加一 add(tmp, cmp); // 把后继节点加入到open表中 printf("can right\n"); } } // 如果当前状态没有后继节点,应该从close表中删除该状态 if(k == 0) { Remove(tmp); // 弹出当前状态 free(tmp); // 释放当前状态 printf("del done\n\n"); } } } int main() { open = initList(); // 初始化open表 close = initList(); // 初始化close表 if(!open || !close) { printf("初始化状态失败\n"); return -1; } // 测试阶段建议直接定义src, init()后的src很大程度上无解 // init(); printf("src:\n"); display(src); Astar(); // 清空open表和close表并回收内存 clearList(open); clearList(close); free(open); free(close); printf("DFS finish\n"); return 0; }
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
已为社区贡献4条内容
所有评论(0)