CC BY 4.0 （除特别声明或转载文章外）
如果这篇博客帮助到你，可以请我喝一杯咖啡~
IDA* Algorithm
Description
Iterative deepening A* (IDA*) was first described by Richard Korf in 1985, which is a graph traversal and path search algorithm that can find the shortest path between a designated start node and any member of a set of goal nodes in a weighted graph.
It is a variant of iterative deepening depth-first search that borrows the idea to use a heuristic function to evaluate the remaining cost to get to the goal from the A* search algorithm.
Since it is a depth-first search algorithm, its memory usage is lower than in A*, but unlike ordinary iterative deepening search, it concentrates on exploring the most promising nodes and thus does not go to the same depth everywhere in the search tree.
Iterative-deepening-A* works as follows: at each iteration, perform a depth-first search, cutting off a branch when its total cost f(n) = g(n)+h(n) exceeds a given threshold. This threshold starts at the estimate of the cost at the initial state, and increases for each iteration of the algorithm. At each iteration, the threshold used for the next iteration is the minimum cost of all values that exceeded the current threshold.
Tasks
- Please solve 15-Puzzle problem by using IDA* (Python or C++). You can use one of the two commonly used heuristic functions: h1 = the number of misplaced tiles. h2 = the sum of the distances of the tiles from their goal positions.
- Here are 4 test cases for you to verify your algorithm correctness. You can also play this game (15puzzle.zip) for more information.
Codes
代码做了如下剪枝：
- 使用所有元素到目标位置的曼哈顿距离之和作为预估函数
- 通过逆序对判断是否有解
- 针对路径判重或判环实测代价比较大，这里仅判断了回路
- 在算 N 数码的逆序数时，不把 0 算入在内；
- 当 N 为奇数时，当两个 N 数码的逆序数奇偶性相同时，可以互达，否则不行；
- 当 N 为偶数时，当两个 N 数码的奇偶性相同的话，那么两个 N 数码中的 0 所在行的差值 k，k 也必须是偶数时，才能互达；
- 当两个 N 数码的奇偶性不同时，那么两个 N 数码中的 0 所在行的差值 k，k 也必须是奇数时，才能互达；
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
#include <bits/stdc++.h>
#define DIST(x, y) (abs((a[x][y] - 1) / 4 - (x)) + abs((a[x][y] - 1) % 4 - (y)))
using namespace std;
int t, a[4][4], ans[64];
int dfs(int cur, int h, int x, int y, int dep)
{
if (cur + h > dep)
return 0;
if (!h)
{
for (int i = 0; i < cur; ++i)
printf("%c", "RLDU"[ans[i]]);
return 1;
}
for (int i = 0, dx[] = {0, 0, 1, -1}, dy[] = {1, -1, 0, 0}, tx, ty; i < 4; ++i)
if (0 <= min(tx = x + dx[i], ty = y + dy[i]) && max(tx, ty) < 4)
if (!cur || i != (ans[cur - 1] ^ 1))
{
ans[cur] = i;
int t = DIST(tx, ty);
swap(a[x][y], a[tx][ty]);
if (dfs(cur + 1, h - t + DIST(x, y), tx, ty, dep))
return 1;
swap(a[x][y], a[tx][ty]);
}
return 0;
}
int main()
{
int t;
for (scanf("%d", &t); t--; printf("\n"))
{
int h = 0, inv = 0, sx, sy;
for (int x = 0; x < 4; ++x)
for (int y = 0; y < 4; ++y)
{
scanf("%d", &a[x][y]);
if (a[x][y])
{
h += DIST(x, y);
for (int i = x << 2 | y, j = 0; j < i; ++j)
inv += a[j / 4][j % 4] > a[x][y];
}
else
sx = x, sy = y;
}
if ((inv + sx) & 1)
for (int dep = h; !dfs(0, h, sx, sy, dep);)
++dep;
else
printf("This puzzle is not solvable.");
}
}
Results
本地测试
input
使用了 TA 演示中的数据。
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
4
11 3 1 7
4 6 8 2
15 9 10 13
14 12 5 0
14 10 6 0
4 9 1 8
2 3 5 11
12 13 7 15
0 5 15 14
7 9 6 13
1 2 12 10
8 11 4 3
6 10 3 15
14 8 7 11
5 1 0 2
13 12 9 4
output
1
2
3
4
LLURRDLURULLURDLLURDDRRULLDLURRDLLDRRULLDRUURULDRRULDRDD
LLDLURDRULDDLUURDDRRUULLDDDLURDRULDLUURRRDLLURDDR
DRDLURULDRDDLUURRDRDLLURURDDLUUURDDDLUUULLDDRDRUUURDDDLULLDRRR
DLLURRURDDLLURURULDRDLULLDDRRULUULDDDRUUURRDLDRD
在线测试
WuK’s solution for [UVA-10181]，运行时间 150ms。