可能出现故障? BBP python算法

时间:2017-02-24 00:18:17

标签: python-2.7

出于某种原因,我的在线和桌面编译环境都认为一切都是无效的语法突然!这是我的代码:

def sumn(n): #summation for sigma
return ((n + 1) * n) / 2
ipt = raw_input('How In Depth Would You Like To Go? ')
ipt = int(ipt)
pi = sumn(ipt) * ((4 / (8 * ipt + 1)) - (2 / (8 * ipt + 4)) - (1 / (8 * ipt + 5)) - (1 / (8 * ipt + 6)) * (1 / (16 ^ ipt))
print pi

2 个答案:

答案 0 :(得分:0)

  1. 您对函数sumn(n)的缩进无效。
  2. 您的行语法ipt = int(raw_input('How In Depth Would You Like To Go? '))不正确。
  3. 供电使用**代替^
  4. 我的代码使用pi的正确缩进和语法:

    >>> def sumn(n): #summation for sigma
            return ((n + 1) * n) / 2
    >>> ipt = float(raw_input('How In Depth Would You Like To Go? '))
    How In Depth Would You Like To Go? 4
    >>> pi = sumn(ipt) * ((4 / (8 * ipt + 1)) - (2 / (8 * ipt + 4)) - (1 / (8 * ipt + 5)) - (1 / (8 * ipt + 6)) * (1 / (16 ** ipt)))
    >>> print pi
    0.386291370825
    

答案 1 :(得分:0)

Yayy,我终于解决了!

from decimal import * 
def sumn(n): #summation for sigma
    return ((n + 1) * n) / 2
ipt = int(raw_input('How In Depth Would You Like To Go? '))
getcontext().prec = ipt


def factorial(n):
    if n<1:
        return 1
    else:
        return n * factorial(n-1)
def BBP(n):
    pi = Decimal(0)
    k = 0
    while k < n:
        pi += (Decimal(1) / (16 ** k)) * ((Decimal(4) / (8 * k + 1)) - (Decimal(2) / (8 * k + 4)) - (Decimal(1) / (8 * k + 5)) - (Decimal(1) / (8 * k + 6)))
        k += 1
    return pi   
print BBP(ipt)

离。

How In Depth Would You Like To Go? 1000
3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117067982148086513282306647093844609550582231725359408128481117450284102701938521105559644622948954930381964428810975665933446128475648233786783165271201909145648566923460348610454326648213393607260249141273724587006606315588174881520920962829254091715364367892590360011330530548820466521384146951941511609433057270365759591953092186117381932611793105118548074462379962749567351885752724891227938183011949129833673362440656643086021394946395224737190702179860943702770539217176293176752384674818467669405132000568127145263560827785771342757789609173637178721468440901224953430146549585371050792279689258923542019956112129021960864034418159813629774771309960518707211349999998372978049951059731732816096318595024459455346908302642522308253344685035261931188171010003137838752886587533208381420617177669147303598253490428755468731159562863882353787593751957781857780532171226806613001927876611195909216420193

最后两个值通常是错误的,我仍在尝试截断它们。然后我所要做的就是在ipt上加2,这将花费我5秒钟的时间☺