c++阿克曼函数详解
不爱吃饭的蓝胖子要开始整活了大家好我是蓝胖子好久不见倍感思念今天带来的是--C阿克曼函数~~希望你能看到最后有惊喜哈正片开始——————————————————————————————————————————标签C、递归、算法、数学、数据结构入门适合人群有 C 基础了解函数与递归概念的初学者前言你有没有想过一个看起来只有三行的函数竟然能让你的电脑算到宇宙热寂都得不出结果今天我们来认识一个计算机科学史上鼎鼎大名的函数——阿克曼函数Ackermann Function。它不仅是递归的极致体现更是可计算性理论中的重要案例。本文将带你从零认识它并用 C 一步步实现。一、什么是阿克曼函数阿克曼函数由德国数学家Wilhelm Ackermann于 1928 年提出其定义如下A(m,n){n1if m0A(m−1,1)if m0 and n0A(m−1,A(m,n−1))if m0 and n0A(m,n)⎩⎨⎧n1A(m−1,1)A(m−1,A(m,n−1))if m0if m0 and n0if m0 and n0其中 mm 和 nn 均为非负整数。乍一看很简单但它有一个令人震惊的特性增长速度极其恐怖。二、手算几个值感受一下我们来手动推算几个小值建立直觉A(0, n)A(0,n)n1A(0,n)n1直接加 1没什么特别的。A(1, n)A(1,0)A(0,1)2A(1,0)A(0,1)2 A(1,1)A(0,A(1,0))A(0,2)3A(1,1)A(0,A(1,0))A(0,2)3 A(1,n)n2A(1,n)n2相当于加 2。A(2, n)A(2,0)A(1,1)3A(2,0)A(1,1)3 A(2,1)A(1,A(2,0))A(1,3)5A(2,1)A(1,A(2,0))A(1,3)5 A(2,2)A(1,A(2,1))A(1,5)7A(2,2)A(1,A(2,1))A(1,5)7 A(2,n)2n3A(2,n)2n3相当于乘 2 再加 3。A(3, n)A(3,n)2n3−3A(3,n)2n3−3增长变成了指数级A(4, n)A(4,0)13,A(4,1)65533,A(4,2)265536−3A(4,0)13,A(4,1)65533,A(4,2)265536−3A(4,2)A(4,2) 是一个位数比宇宙中原子数还多的天文数字小结mm 每增加 1函数增长模式就从加法→乘法→指数→幂塔跃升一个层次。三、C 递归实现根据定义递归实现非常直接cpp复制运行结果A(0, 0) 1 A(1, 1) 3 A(2, 3) 9 A(3, 3) 61 A(3, 4) 125⚠️警告A(4, 0) 13还能算出来但从A(4, 1)开始递归层数爆炸会导致栈溢出Stack Overflow四、栈溢出问题与解决方案为什么会栈溢出每次递归调用都会在调用栈上压入一个栈帧存放局部变量、返回地址等。阿克曼函数的递归深度是超指数级的哪怕是A(3, 10)就需要数百万次递归调用远超系统默认栈大小通常 1MB ~ 8MB。方案一用栈模拟递归非递归版本我们可以用**显式栈std::stack**来模拟递归调用避免系统栈溢出cpp复制这样堆内存可以动态增长比系统栈大得多可以处理更大的输入。方案二记忆化Memoization优化对于重复子问题我们可以用map缓存已计算的结果cpp复制五、可计算的值速查表m \ n01234012345123456235791135132961125413655332^65536 - 3……从表中可以清晰看到增长速度的层次感。六、阿克曼函数有什么用你可能会问这个函数增长这么快实际中有什么用呢可计算性理论阿克曼函数是可计算但不是原始递归函数的经典例子。它证明了可计算比原始递归范围更广是计算理论中的里程碑。算法复杂度分析阿克曼函数的反函数 α(n)增长极其缓慢实际中几乎等于常数出现在**并查集Union-Find**的时间复杂度分析中。你在学并查集时看到 O(α(n)) 的复杂度说的就是它编译器测试阿克曼函数常被用来测试编译器对递归的优化能力如尾递归优化等。压力测试在硬件或语言运行时的压力测试中阿克曼函数是给系统一个极限挑战的常用工具。七、完整代码汇总#include iostream #include stack #include map using namespace std; // 方法一纯递归仅适合小值 long long ackermann_recursive(long long m, long long n) { if (m 0) return n 1; if (n 0) return ackermann_recursive(m - 1, 1); return ackermann_recursive(m - 1, ackermann_recursive(m, n - 1)); } // 方法二显式栈模拟推荐 long long ackermann_iterative(long long m, long long n) { stacklong long stk; stk.push(m); while (!stk.empty()) { m stk.top(); stk.pop(); if (m 0) { n; } else if (n 0) { stk.push(m - 1); n 1; } else { stk.push(m - 1); stk.push(m); n--; } } return n; } // 方法三记忆化递归 mappairlong long, long long, long long memo; long long ackermann_memo(long long m, long long n) { auto key make_pair(m, n); if (memo.count(key)) return memo[key]; long long res; if (m 0) res n 1; else if (n 0) res ackermann_memo(m - 1, 1); else res ackermann_memo(m - 1, ackermann_memo(m, n - 1)); return memo[key] res; } int main() { cout 递归版 endl; cout A(3, 4) ackermann_recursive(3, 4) endl; // 125 cout \n 显式栈版 endl; cout A(3, 6) ackermann_iterative(3, 6) endl; // 509 cout \n 记忆化版 endl; cout A(3, 7) ackermann_memo(3, 7) endl; // 1021 return 0; } //编写不易尊重原创八、总结特性说明定义方式双参数递归定义增长速度超过任意原始递归函数极其恐怖实现挑战深度递归导致栈溢出解决方案显式栈模拟 / 记忆化实际应用可计算性理论、并查集复杂度分析、编译器测试阿克曼函数是简单规则产生极端复杂行为的完美示范。它提醒我们递归的力量不可小觑写代码时一定要估算递归深度否则一个看似无害的函数调用可能让程序瞬间崩溃。如果你觉得这篇文章有帮助欢迎点赞收藏有问题欢迎在评论区交流 本文花费了蓝胖子的大量心血也请大家尊重原创不要转载谢谢
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2472192.html
如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!