python中的递归问题,求圆周率

python中的递归问题,求圆周率

<span>python中的递归问题,求圆周率</span>

以上面一个公式为例:

import numpy as np
def getPi(n):
    if n == 0:
        return np.power(-1,n)*(1.0/(2*n+1))
    else:
        return np.power(-1,n)*(1.0/(2*n+1))+getPi(n-1)
    
print 4*getPi(100)

  

 

可以通过上面一个递归实现。

参考

特点:

①递归就是在过程或者函数里调用自身。

②在使用递归策略时,必须有一个明确的递归条件,称为递归出口。

③递归算法解题通常显得很简洁,但递归算法解题的效率较低。所以一般不倡导使用递归算法设计程序。

④在递归调用的过程当中系统的每一层的返回点、局部变量等开辟了栈来存储。递归函数次数过多容易造成栈溢出等。

   所以一般不倡导用递归算法设计程序。

要求:

递归算法所体现的”重复”一般有三个条件:

①每次在调用规模上都有所缩小(通常是减半)。

②相邻两次重复之间有紧密的联系,前一次要为后一次做准备(通常前一次的输出就作为后一次的输入)。

③在问题的规模极小时必须用直接接触解答而不再进行递归调用,因而每次递归调用都是有条件的(以规模未达到直接解答的大小为条件),

   无条件的递归调用将会成为死循环而不能正常结束。

 

每当你调用一个函数,在这个函数执行前都会将之前的代码地址(也就是调用点)入栈,等被调用的函数执行完将地址出栈,程序根据这个数据返回调用点。
若递归调用次数太多,就会只入栈不出栈,于是堆栈就被压爆了,此为栈溢出。

 

python中的解决办法:

1、人为设置递归深度

import sys
sys.setrecursionlimit(1000000) #括号中的值为递归深度

  事实上并不能完全解决,太多还是会程序崩溃的。

 

2、所谓的尾递归

import numpy as np

def getPi(n,m):
    if n == 0:
        return m+np.power(-1,n)*(1.0/(2*n+1))
    else:
        return getPi(n-1,m+np.power(-1,n)*(1.0/(2*n+1)))
    
print 4*getPi(100,0)

  尾递归参考:https://www.cnblogs.com/hello–the-world/archive/2012/07/19/2599003.html

尾递归的写法就是将操作的值作为参数传递,事实上,python并不支持尾递归的优化!而且对递归的次数有限制,当递归深度超过1000时,会抛出异常。

故对于继续研究突破递归次数的话,虽然有高手找到解决办法,并没有太大意义。

 

3、利用迭代的方式,而不是使用递归(譬如,reduce)

import numpy as np

s = reduce(lambda x,n:x+np.power(-1,n)*(1.0/(2*n+1)),[np.power(-1,0)*(1.0/(2*0+1))]+range(1,100000))
print 4*s

s = reduce(lambda x,n:x+np.power(-1,n)*(1.0/(2*n+1)) if x>0 else np.power(-1,x)*(1.0/(2*x+1))+np.power(-1,n)*(1.0/(2*n+1)),\
          range(0,100000))
print 4*s

  可以看到,利用reduce函数是处理这类问题的比较好的办法。用一句话总结:

                                                  普通程序员用迭代,天才程序员用递归!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/119477.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号