[已解决]如何在python里求pi

妖小圣

今天给大家带来用python求出一个圆形图案的pi，很赏心悦目，感觉还是不错滴，代码分享出来啦，还有源网页：

1. # {{{ http://code.activestate.com/recipes/578130/ (r5)
   
2. def pi(places=10):
   
3.   """Computes pi to given number of decimal places
   
4.     参数places表示要返回的pi的小数点后位数  
   
5.     方法：先整体扩大10**8（10的八次方）倍，然后计算完成后再缩小10的八次方倍
   
6.   """
   
7.    
   
8.   # 3 + 3*(1/24) + 3*(1/24)*(9/80) + 3*(1/24)*(9/80)*(25/168)
   
9.   # The numerators 1, 9, 25, ... are given by (2x + 1) ^ 2
   
10.   # The denominators 24, 80, 168 are given by (16x^2 -24x + 8)
    
11.   extra = 8
    
12.   one = 10 ** (places+extra)
    
13.   t, c, n, na, d, da = 3*one, 3*one, 1, 0, 0, 24
    
14.   #这里的n 和d 分别为每一项的分子与分母 ,na 和 da 分别为分子和分分母后一项比前一项增加的数值
    
15.   #这里的//可不是C++中的注释,而是除的意思
    
16.   while t > 1:
    
17.     n, na, d, da = n+na, na+8, d+da, da+32
    
18.     t = t * n // d
    
19.     c += t
    
20.   return c // (10 ** extra)
    
21. 
    
22. def picirc(radius, aspect_ratio=5):
    
23.   """Display the digit of pi in a circle of given radius
    
24.     radius:显示的半径
    
25.     aspect_ratio：调节显示的比率参数
    
26.   """
    
27.   #display_width为各行的显示长度
    
28.   display_width = int(radius * aspect_ratio + 10)
    
29.   pi_str = repr(pi(int(2 * radius ** 2 * aspect_ratio)))
    
30.   pos = 0
    
31.   #cols为每一行中要显示的数字个数
    
32.   for i in range(2 * radius):
    
33.     cols = int(0.5 + aspect_ratio * (radius**2 - (radius-(i+0.5))**2) ** 0.5)
    
34.     print(pi_str[pos:pos+cols].center(display_width)) #将产生的pi数值生成的文本列表中相应的位数取出来显示在当前行
    
35.     pos += cols
    
36. 
    
37. if __name__ == '__main__':
    
38.    
    
39.   picirc(16)
    
40. ## end of http://code.activestate.com/recipes/578130/ }}}

复制代码

最佳答案

月排行榜 / 总排行榜

jerryxjr1220

2017-4-8 23:51:40

我之前的小练习已经做过了
http://bbs.fishc.com/thread-80041-1-1.html

python2的代码：

1. def calcPi(n):
   
2.     w = 10**(n+10)
   
3.     x = int(16 * w / 5 - 4 * w / 239)
   
4.     for i in range(4,2*n,2):
   
5.         x1 = 16 * w / (i-1)
   
6.         x1 /= 5**(i-1)
   
7.         x2 = 4 * w / (i-1)
   
8.         x2 /= 239**(i-1)
   
9.         x += ((-1)**(i/2+1))*(x1-x2)
   
10.     return x // 10**10

复制代码

跳转到最佳答案楼层

～风介～

无图无真相~

jerryxjr1220

我之前的小练习已经做过了
http://bbs.fishc.com/thread-80041-1-1.html

python2的代码：

1. def calcPi(n):
   
2.     w = 10**(n+10)
   
3.     x = int(16 * w / 5 - 4 * w / 239)
   
4.     for i in range(4,2*n,2):
   
5.         x1 = 16 * w / (i-1)
   
6.         x1 /= 5**(i-1)
   
7.         x2 = 4 * w / (i-1)
   
8.         x2 /= 239**(i-1)
   
9.         x += ((-1)**(i/2+1))*(x1-x2)
   
10.     return x // 10**10

复制代码