Python中的1000位pi
Python中的1000位pi
提问于 2012-01-25 22:59:48
我一直在思考这个问题，但我想不出来。也许你能帮我。问题是我的代码不能用Python编码语言输出1000位pi。
下面是我的代码：
def make_pi():
q, r, t, k, m, x = 1, 0, 1, 1, 3, 3
while True:
if 4 * q + r - t < m * t:
yield m
q, r, t, k, m, x = (10*q, 10*(r-m*t), t, k, (10*(3*q+r))//t - 10*m, x)
else:
q, r, t, k, m, x = (q*k, (2*q+r)*x, t*x, k+1, (q*(7*k+2)+r*x)//(t*x), x+2)
digits = make_pi()
pi_list = []
my_array = []
for i in range(1000):
my_array.append(str("hello, I'm an element in an array \n" ))
big_string = "".join(my_array)
print "here is a big string:\n %s" % big_string 我知道这段代码可以修复，但我不确定该修复什么……这里的print语句是一个很大的字符串，而my_array.append(str("hello, im an element in an array \n))现在只是一个填充符。我知道所有的代码是如何工作的，但就像我之前说过的，我不能让它输出这些代码。
回答 11
Stack Overflow用户
回答已采纳
发布于 2012-01-25 23:22:39
运行此命令
def make_pi():
q, r, t, k, m, x = 1, 0, 1, 1, 3, 3
for j in range(1000):
if 4 * q + r - t < m * t:
yield m
q, r, t, k, m, x = 10*q, 10*(r-m*t), t, k, (10*(3*q+r))//t - 10*m, x
else:
q, r, t, k, m, x = q*k, (2*q+r)*x, t*x, k+1, (q*(7*k+2)+r*x)//(t*x), x+2
my_array = []
for i in make_pi():
my_array.append(str(i))
my_array = my_array[:1] + ['.'] + my_array[1:]
big_string = "".join(my_array)
print "here is a big string:\n %s" % big_string 并从此处阅读有关yield operator的信息：What does the "yield" keyword do?
这就是答案：
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211055596446229489549303819644288109756659334461284756482337Stack Overflow用户
发布于 2012-11-10 06:30:53
Stack Overflow用户
发布于 2019-03-03 17:25:07
接受的答案是不正确的，如注释中所述。
OP的代码似乎是基于从here复制的Spigot's algorithm实现。
要根据OP的问题修复代码(尽管我重命名了变量和函数以与原始源代码中的变量和函数匹配)，一种解决方案可能是：
#!/usr/bin/env python
DIGITS = 1000
def pi_digits(x):
"""Generate x digits of Pi."""
q,r,t,k,n,l = 1,0,1,1,3,3
while x >= 0:
if 4*q+r-t < x*t:
yield n
x -= 1
q,r,t,k,n,l = 10*q, 10*(r-n*t), t, k, (10*(3*q + r))/t-10*n, l
else:
q,r,t,k,n,l = q*k, (2*q+r)*l, t*l, k+1, (q*(7*k+2)+r*l)/(t*l), l+2
digits = [str(n) for n in list(pi_digits(DIGITS))]
print("%s.%s\n" % (digits.pop(0), "".join(digits)))此外，这里有一个更快的*实现，显然也是基于Spigot的算法：
#!/usr/bin/env python
DIGITS = 1000
def pi_digits(x):
"""Generate x digits of Pi."""
k,a,b,a1,b1 = 2,4,1,12,4
while x > 0:
p,q,k = k * k, 2 * k + 1, k + 1
a,b,a1,b1 = a1, b1, p*a + q*a1, p*b + q*b1
d,d1 = a/b, a1/b1
while d == d1 and x > 0:
yield int(d)
x -= 1
a,a1 = 10*(a % b), 10*(a1 % b1)
d,d1 = a/b, a1/b1
digits = [str(n) for n in list(pi_digits(DIGITS))]
print("%s.%s\n" % (digits.pop(0), "".join(digits)))我在this在线Pi数字生成器上对这两个工具进行了多次测试。
*基于测试10,000个数字，其中我得到了大约7秒，而不是大约1秒。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/9004789
复制相关文章
相似问题
腾讯云开发者
