Python 最简代码实现万位高精度圆周率π的计算

利用Python的特殊功能,简单实现万位高精度圆周率计算。

啥也别说了,直接给代码:

# -*- coding: UTF-8 -*-
# calculating PI with a simple series
# Author: Idealguy,Shanghai,2018
# pi=2*(1+1/3+1/3*2/5+1/3*2/5+3/7+...)
#
n=10000+4
p=2*10**n
a=p//3; p+=a
i=2
while a>0:
    a=a*i//(i*2+1); i+=1
    p+=a
p//=10000
print ("PI=",p)

计算结果(on Python 3.6.3 )

>>> ==================== RESTART: D:\Python\math\Pi-Simp2.py ==================== 
PI= 
3141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117067
982148086513282306647093844609550582231725359408128481117450284102701938521105559644622948954930381
  • 3
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
A: 首先,需要明确一下什么是串行方法和并行计算。串行方法指在一个线程内按照顺序依次执行计算任务,而并行计算则是利用多个线程或多个计算单元同时进行计算,提高计算效率。 接下来是利用公式计算圆周率π的方法。常见的公式有莱布尼兹级数公式、马青公式、蒙特卡罗方法等。这里以莱布尼兹级数公式为例进行说明。 莱布尼兹级数公式:π/4 = 1 - 1/3 + 1/5 - 1/7 + 1/9 - ... 其中,公式右侧的每一项都可以作为一个独立的计算任务。串行方法的实现就是按照公式依次计算每一项的值,最终将它们相加得到π的值。代码如下: ```python def serial_pi(n): pi = 0.0 sign = 1.0 for i in range(n): pi += sign / (2 * i + 1) sign = -sign return pi * 4 ``` 其中,n为计算次数,次数越多计算精度越高但计算时间也会越长。 如果要进行并行计算,则可以将计算任务分配给多个线程或计算单元同时处理。比如,可以把公式右侧的每一项分给不同的线程处理,最终得到每个线程的计算结果后再将它们相加得到π的值。代码如下: ```python import threading pi = 0.0 lock = threading.Lock() def parallel_pi(start, end, step): global pi sign = 1.0 local_pi = 0.0 for i in range(start, end, step): local_pi += sign / (2 * i + 1) sign = -sign with lock: pi += local_pi def calc_pi(n, num_threads): global pi pi = 0.0 threads = [] for i in range(num_threads): start = i end = n step = num_threads t = threading.Thread(target=parallel_pi, args=(start, end, step)) threads.append(t) t.start() for t in threads: t.join() return pi * 4 ``` 其中,num_threads为线程数量,每个线程计算的范围是从start到end,步长为step。 以上就是单机并行计算圆周率的一个实现方法。需要注意的是,因为多线程之间的计算是并行进行的,所以需要采用线程安全的方式来保证计算结果的准确性,比如使用锁机制。同时,多线程计算的效率也受到计算单元的数量和负载均衡的影响,需要进行适当的优化才能发挥最大的性能优势。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值