我想很多人第一次学习递归的时候，老师或者书本上可能会举汉诺塔的例子。
但是今天，我们讨论的重点不是简单的汉诺塔算法，而是三柱汉诺塔的延伸。先来看看经典的三柱汉诺塔。
一、三柱汉诺塔（Hanoi_Three）：
我想大家对于三柱汉诺塔的理解以及算法的实现应该是很熟练了。
我在这里简单的过一遍三柱汉诺塔的算法思想：
有A、B、C三根柱子，A柱上有n个盘子，现在需要将A上所有的盘子转移到C上，请给出搬运次数最少的步骤。
算法思想：
1、将A上n-1个盘子以C为缓存，全部转移到 B 柱上。
2、将A上留下的第n个盘子，直接转移到 C 柱上。
3、将B上的n-1个盘子，以A为缓存，全部转移到 C 柱上。
很容易得到算法的递归方程为：T(n)=2*T(n-1)+1，不难算出步数是T(n)=2^n-1。
具体的代码如下：
void Move( char x, char y )
{
printf("%c --> %c\n", x, y); //打印路径，x-->y
}
void Hanoi_Three( int n, char a, char b, char c )
{
if( n <= 0 )
return ;
step++; //步骤+1
if( n == 1 )
{
Move( a, c ); //将a柱上的一个盘子直接移动到c柱
return ;
}
else
{
Hanoi_Three( n-1, a, c, b); //将a上n-1个盘子以c为缓存，全部转移到 b 柱上
Move( a, c ); //将A上留下的第n个盘子，直接转移到 C 柱上
Hanoi_Three( n-1, b, a, c); //将b上的n-1个盘子，以a为缓存，全部转移到 c 柱上
}
}
二、四柱汉诺塔（Hanoi_Four）：
当柱子为四根时，对于只是将盘子全部转移到另一根柱子上这个目的来说，是大大降低了难度，而且算法的复杂度也大大降低了。但是，这个时候，如果要你找到一个最优的、步骤最少的实现方法，可以说难度是提升了一个数量级。
有些人可能会质疑，为什么，我用三柱汉诺塔的思想不是很优化了吗？
别急，且让我慢慢向你道来。
先来看看这种‘看上去很合理’的解法：
假设，A，B，C，D，分别为：源位置，缓存，缓存，目的位置。
因为三柱的时候，我们是将A的前n-1个盘子放到B上缓存，然后将第n个盘子放到C柱上。
现在的情况好很多，有两个可以缓存的柱子，因此，看上去移动起来更加方便，原来B上需要缓存的n-1个盘子，现在可以只是n-2个盘子，而将第n-1个盘子放到C上缓存。
具体的流程如下（非最优解法）：
1、从A借助C、D将 n-2个盘子移动到B上。
2、将第n-1个盘子移动到C上。
3、将第n个盘子移动到D上。
4、将第n-1个盘子移动到D上。
5、从B借助A、C将 n-2个盘子全部移动到D上。
看上去，非常完美，笔者也一度觉得这个思想没有破绽，甚至我还自以为找到了k根柱子汉诺塔的通用方法（想当然的将B柱上缓存的数量从n-2个，改为n-（k-2）个盘子）。直到我看了这篇文章：多柱汉诺塔最优算法设计探究。
虽然我们想到让盘子尽量不发生重叠来保证步数的最少，但是这并不能绝对保证。或许在盘子较少的情况下是可行的，但是盘子增多时，那些多余的只有一个盘子的柱子是可以加以利用的（可能的优化在这里）。虽然这么做加多了每次的移动步数，但是却从另一个侧面减少了递归的数量，因此我们需要从这里边找一个平衡点。
下面我们来看看，1941年，美国的J. S. Frame，给出的四柱汉诺塔的算法思想，也叫Frame算法：
1、用4柱汉诺塔算法把A柱上部分的n- r个碟子通过C柱和D柱移到B柱上【F( n- r )步】。
2、用3柱汉诺塔经典算法把A柱上剩余的r个碟子通过C柱移到D柱上【2^r-1步】（参照上述三柱时的情况）。
3、用4柱汉诺塔算法把B柱上的n-r个碟子通过A柱和C柱移到D柱上【F(n-r)步】。
4、依据上边规则求出所有r（1≤r≤n）情况下步数f(n)，取最小值得最终解。
因此Frame算法的递归方程如下：
F(n)=min(2*F(n-r)+2^r-1)，（1≤r≤n）。
大家有没有发现，其实，这个算法思想跟我们之前认为合理的算法基本一致，差别只是在于他将我们的n-2个碟子缓存到B上，改为了将n- r个碟子转移到B柱上。
差别即使核心，这个算法的核心，就是计算n个盘子的情况下，r为何值时，能够使得算法最优。
找到了核心，我们现在的任务就明确了，就是对r值的计算。
这里给出了一个较笨的方法--枚举（不知道各位有没有其他方法）。就是将一定范围内的n与r的所有取值带入，得到满足F（n）为最小值的r的值，记为K[n] = r；
具体的代码如下：
得到各个n对于的r之后，算法将变的非常简单，具体实现如下：
#include
#include
#define Max 100
#define INT_MAX 0xfffffffffffffff
int K[Max+1] = {0};
int step = 0;
void Hanoi_Four( int n, char a, char b, char c, char d );
void Hanoi_Three( int n, char a, char b, char c );
void Move( char x, char y );
void Move( char x, char y )
{
printf("%c --> %c\n", x, y); //打印路径，x-->y
}
void Hanoi_Three( int n, char a, char b, char c )
{
if( n <= 0 )
return ;
step++; //步骤+1
if( n == 1 )
{
Move( a, c ); //将a柱上的一个盘子直接移动到c柱
return ;
}
else
{
Hanoi_Three( n-1, a, c, b); //将a上n-1个盘子以c为缓存，全部转移到 b 柱上
Move( a, c ); //将A上留下的第n个盘子，直接转移到 C 柱上
Hanoi_Three( n-1, b, a, c); //将b上的n-1个盘子，以a为缓存，全部转移到 c 柱上
}
}
void Hanoi_Four( int n, char a, char b, char c, char d )
{
if( n <= 0 )
return ;
if( n == 1 )
{
step++;
Move( a, d );
return ;
}
else
{
int kn = K[n];
//printf("kn = %d\n", K[n]);
Hanoi_Four( n-kn, a, c, d, b ); //用4柱汉诺塔算法把A柱上部分的n- kn个碟子通过C柱和D柱移到B柱上
Hanoi_Three( kn, a, c, d ); //用3柱汉诺塔经典算法把A柱上剩余的kn个碟子通过C柱移到D柱上。
Hanoi_Four( n-kn, b, a, c, d ); //用4柱汉诺塔算法把B柱上的n-r个碟子通过A柱和C柱移到D柱上
}
}
void Init_K(void )
{
int i, k;
__int64 temp;
__int64 m[Max+1] = {0};
for( i = 1; i <= Max; i++ )
{
m[i] = INT_MAX;
for( k = 1; k <= i; k++ )
{
temp = 2*m[i-k] + (__int64)pow(2,k) - 1;
if( temp < m[i] )
{
m[i] = temp;
K[i] = k;
//printf("K[%d] = %d, m[i] = %d\n", i, k, temp );
}
}
}
}
int main()
{
int n;
Init_K();
printf("Please enter the number of the Plates: \n");
while( scanf("%d", &n) != EOF )
{
step = 0;
Hanoi_Four( n, 'A', 'B', 'C', 'D' );
//Hanoi_Three( n, 'A', 'B', 'C' );
printf("**************************\nTotal Step: %d\n", step );
printf("Please enter the number of the Plates: \n");
}
return 0;
}
到这里，四柱汉诺塔的算法基本讲完了。
有兴趣的同学，可以继续归纳多柱汉诺塔的实现方法，欢迎交流指导！
很多时候，看似合理的背后，其实是一种思维定势。。。
--------------------------by-----------cyh24----------2012年10月5日22:54:17
