【算法专题】汉诺塔问题

汉诺塔问题

有三根相邻的柱子，标号为A,B,C，A柱子上从下到上按金字塔状叠放着n个不同大小的圆盘，要把所有盘子一个一个移动到柱子C上，并且每次移动同一根柱子上都不能出现大盘子在小盘子上方，请问至少需要多少次移动以及如何移动？

一. 基本汉诺塔问题

如何移动？

• 这是一个递归问题。

• 假设移动有n个圆盘，最终从A移动到C： 可以先将上面的n-1个圆盘移动到B，然后把最下面的圆盘移动到C，最后再把n-1根圆盘移动到C即可。

• 具体代码如下：

/**
 * 汉诺塔求具体方案
 */
public class Main1 {

    // 将 a 上的圆盘借助 b 最终移动到 c 上
    private void hanoi(int n, char a, char b, char c) {

        if (n == 1) System.out.println(a + " -> " + c);
        else {
            hanoi(n - 1, a, c, b);  // n-1个圆盘从a移动到b上
            System.out.println(a + " -> " + c);  // 最下面的圆盘从a移动到c上
            hanoi(n - 1, b, a, c);  // n-1个圆盘从b移动到c上
        }
    }

    public static void main(String[] args) {

        Main1 t = new Main1();
        t.hanoi(3, 'A', 'B', 'C');
    }
}

• 运行结果：

A -> C
A -> B
C -> B
A -> C
B -> A
B -> C
A -> C

移动次数？

• 结论：n个圆盘的汉诺塔需要的移动次数为$2^n-1$。分析如下：
• 先将上面的n-1个圆盘移动到B，然后把最下面的圆盘移动到C，最后再把n-1根圆盘移动到C即可，因此有递推式：

$$d[i]=d[i-1]+1+d[i-1]=2\times d[i-1]+1$$

• 具体代码如下：

/**
 * 求解汉诺塔移动次数
 */
public class Main2 {

    public static final int N = 15;  // 数据范围
    public static int[] d = new int[N];  // 存储最小移动次数

    private static void init() {
        d[1] = 1;
        for (int i = 2; i < N; i++)
            d[i] = 2 * d[i - 1] + 1;
    }

    public static void main(String[] args) {

        init();

        System.out.println(d[3]);  // 三个圆盘需要的最小移动次数
    }
}

• 运行结果：

7

二. 拓展问题

AcWing 96. 奇怪的汉诺塔

问题描述

• 问题链接：AcWing 96. 奇怪的汉诺塔

分析

• 如果是三个柱子，最终从A移动到C: 可以先将上面的n-1个圆盘移动到B，然后把最下面的圆盘移动到C，最后再把n-1根圆盘移动到C即可，因此有递推式：$d[i]=d[i-1]+1+d[i-1]=2\times d[i-1]+1$。

• 如果是四个柱子，最终从A移动到D:（1）选择将j个圆盘移动到B(四根柱子的状态下)，（2）将剩余的i-j个圆盘移动到D(三柱子状态下)，（3）将j个圆盘移动到D(四根柱子的状态下)。因此有递推式：$f[i]=f[j]+d[i-j]+f[j]=f[j]\times 2+d[i-j]$。

代码

• C++

#include <iostream>
#include <cstring>

using namespace std;

const int N = 15;

int d[N];  // 三根柱子对应方案数
int f[N];  // 四根柱子对应方案数

int main() {

    d[1] = 1;
    for (int i = 2; i <= 12; i++) d[i] = d[i - 1] * 2 + 1;

    memset(f, 0x3f, sizeof f);
    f[0] = 0;
    for (int i = 1; i <= 12; i++)
        for (int j = 0; j < i; j++)
            f[i] = min(f[i], f[j] * 2 + d[i - j]);

    for (int i = 1; i <= 12; i++) cout << f[i] << endl;

    return 0;
}