python求圆周率马青公式_Python 实现丘德诺夫斯基（Chudnovsky）法計算高精度圓周率...

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

源码：(需要将 @ 替换成 ASCII空格)

# -*- coding: UTF-8 -*-

# 丘德诺夫斯基法計算高精度圓周率程序

# Calculating PI with Chudnovsky-Series

# Author: Idealguy，2018, Shanghai

import time

# In following functions, High-Prec Nums are both amplified 10**n

# pre-defined: Base=10**n

def Sqrt10005(): # Sqrt(10005L) by Imitate-Manual Method

@@@@n1=0

@@@@c=10002499687 #100.02499687

@@@@mc=8; m=mc

@@@@f1=10**mc

@@@@f2=f1*f1

@@@@a=10005*f2-c*c

@@@@while mc

@@@@@@@@a*=f2

@@@@@@@@b=c*2*f1

@@@@@@@@d=a//b

@@@@@@@@c=c*f1+d

@@@@@@@@a-=d*(b+d)

@@@@@@@@mc+=m

@@@@@@@@if mc*2>n: m=n-mc

@@@@@@@@else: m=mc

@@@@@@@@f1=10**m

@@@@@@@@f2=f1*f1

@@@@@@@@n1+=1

@@@@return c

# Main Program

print ("Chudnovsky法計算高精度圓周率程序")

while 1:

@@@@n=int(input('計算位數[1..50000],0:退出：'))

@@@@if n<=10: break

@@@@n+=2

@@@@base=10**n

@@@@t=time.clock()

@@@@# Start Calculating

@@@@A=13591409*base; B=A

@@@@c3=13591409

@@@@i=1

@@@@while abs(A)>5:

@@@@@@@@c1=((108-72*i)*i-46)*i+5

@@@@@@@@c2=10939058860032000*i**3

@@@@@@@@c4=c3; c3+=545140134

@@@@@@@@i+=1

@@@@@@@@A=A*c1*c3//(c2*c4) # Must in form: A=A*...

@@@@@@@@B+=A

@@@@p=426880*base*Sqrt10005()//B//100

@@@@# Post access

@@@@print ("用時= %8.3f 秒" % (time.clock()-t))

@@@@s=input('是否显示结果(Y/N):')

@@@@if (s=='Y')|(s=="y"): print ("PI="+str(p))

# end

python求圆周率马青公式_Python 实现丘德诺夫斯基（Chudnovsky）法計算高精度圓周率...相关推荐

Python 实现丘德诺夫斯基（Chudnovsky）法計算高精度圓周率
[原创] 本文利用Python自带的高精度整型计算功能,实现丘德诺夫斯基(Chudnovsky)法計算高精度圓周率.计算速度比马青公式要快. 经过艰苦的调试,实现该计算,10万位的计算在秒级完成.以下 ...
Python:实现chudnovsky algorithm楚德诺夫斯基算法(附完整源码)
Python:实现chudnovsky algorithm楚德诺夫斯基算法 from decimal import Decimal, getcontext from math import ceil, ...
python求伽马分布模型参数_Python实现的各种常见分布算法示例
本文实例讲述了Python实现的各种常见分布算法.分享给大家供大家参考,具体如下: #-*- encoding:utf-8 -*- import numpy as np from scipy impo ...
Python利用马青公式计算圆周率Π并写入文件
一.什么是马青公式马青公式由英国天文学教授约翰·马青(John Machin ,1686 –1751)于1706年发现,他利用这个公式计算到了100位的圆周率. 马青公式每计算一项可 ...
【Python、数学】计算任意位数的圆周率π（马青公式）
1. 公式准备计算准确圆周率的马青公式: 对反正切进行级数展开: 就可以得到 π = 16(1/5 - 1/3/5^3 + 1/5/5^5 - ...) - 4(1/239 - 1/3/239^3 ...
用马青公式计算圆周率，Python语言
马青公式: 简单实现 import times_time = time.time() #定义计算的位数 num = 100000 #多计算10位,以防出错 num1 = num + 10 #定义计算小 ...
利用马青公式输出π的后任意位数字
马青公式 π=16arctan15−4arctan1239π=16arctan15−4arctan1239\pi = 16arctan \frac{1}{5} - 4arctan \frac{1}{2 ...
利用Java的BigDecimal与马青公式精确计算π后10000位，
首先给出公式如下: π=16arctan1/5−4arctan1/239: 即是 π=16×(1/(1×5)−1/(3×5的3次方)+1/(5×5的5次方)-)−4×(1/(1×239)−1/(3×2 ...
用马青公式计算π的c程序
/*一个计算π的c程序 ====================================== 由于大多数计算机内置算法有一定精度限制,你想计算那么多位就会产生问题. 这里有一个c程序,允许计算 ...

python求圆周率马青公式_Python 实现丘德诺夫斯基（Chudnovsky）法計算高精度圓周率...

python求圆周率马青公式_Python 实现丘德诺夫斯基（Chudnovsky）法計算高精度圓周率...相关推荐

最新文章

热门文章