动态模拟退火法解决超高阶魔方还原问题.PDFVIP

58 《高性能计算发展与应用》 2013年第二期 总第四十三期 动态模拟退火法解决超高阶魔方还原问题  陈 希 合肥微尺度物质科学国家实验室 中国科学技术大学物理系 合肥230026  李会民 中国科学技术大学超级运算中心 合肥230026  丁泽军 合肥微尺度物质科学国家实验室 中国科学技术大学物理系 合肥 摘要： 本论文提出了一种有较高效率的Monte Carlo算法——动态模拟退火算法解决任意超高阶魔 方 （Rubik’s Cube）还原问题。通过设计合适的能量函数、冷却过程和邻域搜索算法，移动序 列可以形成一条快速降低魔方无序度的路径，并用简洁高效的方法跳出能量局部极小点。与静 态模拟退火法中通过改变Boltzmann函数中的温度参数不同的是，该方法利用动态切换能量函数 跳出局部极小点。另外，本论文设计了一种以低阶魔方的计算用于高阶时的高并行度的并行算 法，使得超高阶魔方也可以在较短时间内被还原，其还原时间近似与阶数的平方成正比。典型 的计算时间为，笔记本电脑上求解101阶魔方用时996.6秒，用8个节点的计算机集群和并行程序 还原5001阶魔方用时1877秒。高阶魔方还原算法中涉及的原理，如分阶段法、恰当能量函数的 设计方式、如何找到与能量函数匹配的邻域搜索等，可能对其它的全局优化问题提供帮助，特 别是在避免陷入能量绘景中大量的局部极小点方面。 1. 引言 题中作为搜寻目标终点的系统基态是已知的而且是 自从1974年匈牙利的Ernő Rubik发明魔方（一 唯一的，而在其它全局优化问题中基态是未知的。 [1] 因此，问题归结为在无数条连接初态与基态的路径 种三维机械玩具）以来，寻找其最优解 、群论解 [2] 中找到一条最短的或较短的，其计算难度与其它类 法 和魔方在物理上的隐喻引起了人们的极大兴趣。 然而，直到今天，其研究还仅局限于低阶魔方，多 似的全局优化问题相当，只要可用的计算资源是有 数是从3到11阶。不管是理论上还是采用计算机还 限的。如果不能恰当地设计搜索方案，很容易在未 原，对高阶魔方和超高阶魔方的尝试很少。魔方问 达到求解目标之前就耗尽了计算资源。 题的主要困难在于，微观状态数A的数量级随阶数 魔方可以类比为一个六面的二维Ising模型，每 n的阶乘而增长，求解难度急剧攀升。例如，即使对 一小面的颜色可当作自旋矢量的一个方向，魔方最 有序的状态对应着Ising模型中零温下的铁磁相。魔 于n=5的低阶魔方，A=282,870,942,277,741,856,536,1 方求解转化为寻找一条通往能量最小的稳定态的路 80,333,107,150,328,293,127,731,985,672,134,721,536,0 [3] 径问题，只要把魔方的能量定义成是随无序度增加 00,000,000,000,000。所以，对到一定高阶的魔方， 如果让计算机从任意混乱状态开始，随机转动寻找 的。根据算法，接受的随机试探步都倾向于降低魔 初始态而完全没有还原方向的指引，所耗的时间将 方的能量。这种有导向的路径构