请选择 进入手机版 | 继续访问电脑版

技术控

    今日:34| 主题:61586
收藏本版 (1)
最新软件应用技术尽在掌握

[其他] Python开启尾递归优化!

[复制链接]
哦吼啦~] 发表于 2016-11-30 12:39:34
189 2
一般递归与尾递归

  一般递归

  1. def normal_recursion(n):
  2.     if n == 1:
  3.         return 1
  4.     else:
  5.         return n + normal_recursion(n-1)
复制代码
执行:
  1. normal_recursion(5)
  2. 5 + normal_recursion(4)
  3. 5 + 4 + normal_recursion(3)
  4. 5 + 4 + 3 + normal_recursion(2)
  5. 5 + 4 + 3 + 2 + normal_recursion(1)
  6. 5 + 4 + 3 + 3
  7. 5 + 4 + 6
  8. 5 + 10
  9. 15
复制代码
可以看到, 一般递归, 每一级递归都需要调用函数, 会创建新的栈,
  随着递归深度的增加, 创建的栈越来越多, 造成爆栈:boom:
  尾递归

   尾递归基于函数的尾调用 , 每一级调用直接返回函数的返回值更新调用栈,而不用创建新的调用栈, 类似迭代的实现, 时间和空间上均优化了一般递归!
  1. def tail_recursion(n, total=0):
  2.     if n == 0:
  3.         return total
  4.     else:
  5.         return tail_recursion(n-1, total+n)
复制代码
执行:
  1. tail_recursion(5)
  2. tail_recursion(4, 5)
  3. tail_recursion(3, 9)
  4. tail_recursion(2, 12)
  5. tail_recursion(1, 14)
  6. tail_recursion(0, 15)
  7. 15
复制代码
可以看到, 每一级递归的函数调用变成"线性"的形式.
  深入理解尾递归

  呃, 所以呢? 是不是感觉还不够过瘾... 谁说尾递归调用就不用创建新的栈呢?
  还是让我们去底层一探究竟吧
  1. int tail_recursion(int n, int total) {
  2.     if (n == 0) {
  3.         return total;
  4.     }
  5.     else {
  6.         return tail_recursion(n-1, total+n);
  7.     }
  8. }
  9. int main(void) {
  10.     int total = 0, n = 4;
  11.     tail_recursion(n, total);
  12.     return 0;
  13. }
复制代码
反汇编

  
       
  • $ gcc -S tail_recursion.c -o normal_recursion.S
       
  •   $ gcc -S -O2 tail_recursion.c -o tail_recursion.S gcc开启尾递归优化
      
  对比反汇编代码如下(AT&T语法)
   

Python开启尾递归优化!

Python开启尾递归优化!-1-技术控-return,空间

  可以看到, 开启尾递归优化前, 使用call调用函数, 创建了新的调用栈(LBB0_3);
  而开启尾递归优化后, 就没有新的调用栈生成了, 而是直接pop
   bp指向的 _tail_recursion 函数的地址(pushq %rbp)然后返回,
   

Python开启尾递归优化!

Python开启尾递归优化!-2-技术控-return,空间

  仍旧用的是同一个调用栈!
  存在的问题

  虽然尾递归优化很好, 但python 不支持尾递归,递归深度超过1000时会报错
  1. RuntimeError: maximum recursion depth exceeded
复制代码
一个牛人想出的解决办法

  实现一个 tail_call_optimized 装饰器

  1. #!/usr/bin/env python2.4
  2. # This program shows off a python decorator(
  3. # which implements tail call optimization. It
  4. # does this by throwing an exception if it is
  5. # it's own grandparent, and catching such
  6. # exceptions to recall the stack.
  7. import sys
  8. class TailRecurseException:
  9.     def __init__(self, args, kwargs):
  10.         self.args = args
  11.         self.kwargs = kwargs
  12. def tail_call_optimized(g):
  13.     """
  14.     This function decorates a function with tail call
  15.     optimization. It does this by throwing an exception
  16.     if it is it's own grandparent, and catching such
  17.     exceptions to fake the tail call optimization.
  18.     This function fails if the decorated
  19.     function recurses in a non-tail context.
  20.     """
  21.     def func(*args, **kwargs):
  22.         f = sys._getframe()
  23.         # 为什么是grandparent, 函数默认的第一层递归是父调用,
  24.         # 对于尾递归, 不希望产生新的函数调用(即:祖父调用),
  25.         # 所以这里抛出异常, 拿到参数, 退出被修饰函数的递归调用栈!(后面有动图分析)
  26.         if f.f_back and f.f_back.f_back \
  27.             and f.f_back.f_back.f_code == f.f_code:
  28.             # 抛出异常
  29.             raise TailRecurseException(args, kwargs)
  30.         else:
  31.             while 1:
  32.                 try:
  33.                     return g(*args, **kwargs)
  34.                 except TailRecurseException, e:
  35.                     # 捕获异常, 拿到参数, 退出被修饰函数的递归调用栈
  36.                     args = e.args
  37.                     kwargs = e.kwargs
  38.     func.__doc__ = g.__doc__
  39.     return func
  40. @tail_call_optimized
  41. def factorial(n, acc=1):
  42.     "calculate a factorial"
  43.     if n == 0:
  44.         return acc
  45.     return factorial(n-1, n*acc)
  46. print factorial(10000)
复制代码
  为了更清晰的展示开启尾递归优化前、后调用栈的变化和tail_call_optimized装饰器抛异常退出递归调用栈的作用, 我这里利用 pudb调试工具 做了动图 <br/>
  开启尾递归优化前的调用栈
   

Python开启尾递归优化!

Python开启尾递归优化!-3-技术控-return,空间

  开启尾递归优化后(tail_call_optimized装饰器)的调用栈
   

Python开启尾递归优化!

Python开启尾递归优化!-4-技术控-return,空间

  通过pudb右边栏的stack, 可以很清晰的看到调用栈的变化.
   因为尾递归没有调用栈的嵌套, 所以Python也不会报 RuntimeError: maximum recursion depth exceeded 错误了!
  这里解释一下 sys._getframe() 函数:
  1. sys._getframe([depth]):
  2. Return a frame object from the call stack.
  3. If optional integer depth is given, return the frame object that many calls below the top of the stack.
  4. If that is deeper than the call stack, ValueEfror is raised. The default for depth is zero,
  5. returning the frame at the top of the call stack.
  6. 即返回depth深度调用的栈帧对象.
  7. import sys
  8. def get_cur_info():
  9.     print sys._getframe().f_code.co_filename  # 当前文件名
  10.     print sys._getframe().f_code.co_name  # 当前函数名
  11.     print sys._getframe().f_lineno # 当前行号
  12.     print sys._getframe().f_back # 调用者的帧
复制代码
  更多关于 sys._getframe 请看 Frame Hacks
liu2196 发表于 2016-11-30 13:03:27
看来楼下的有话说!
回复 支持 反对

使用道具 举报

南生 发表于 2016-11-30 17:29:57
哦吼啦~]的内心很强大!
回复 支持 反对

使用道具 举报

我要投稿

推荐阅读


回页顶回复上一篇下一篇回列表
手机版/c.CoLaBug.com ( 粤ICP备05003221号 | 文网文[2010]257号 | 粤公网安备 44010402000842号 )

© 2001-2017 Comsenz Inc.

返回顶部 返回列表