用python计算圆周率

一、算法

马青公式:马青公式由英国天文学教授约翰·马青(John Machin ,1686 –1751)于1706年发现。马青公式每计算一项可以得到1.4位的十进制精度。因为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机上编程实现。

π=16arctan15−4arctan1239π=16arctan15−4arctan1239 arctanx=x−x33+x55−......arctanx=x−x33+x55−...... π=16×(11×5−13×53+......)−4×(11×239−13×2393+......)π=16×(11×5−13×53+......)−4×(11×239−13×2393+......)
 

 

二、python代码

根据提示输入要计算的长度,计算完成自动生成pai.txt文件,要多长有多长。

经过测试生成10w位大概需要23秒左右。


# -*- coding: utf8 -*-
import time
def comput():n =int(input('请输入要计算的长度:'))start_time = time.time()w = n+10b = 10**wx1 = b*4//5x2 = b// -239he = x1+x2n *= 2for i in range(3,n,2):x1 //= -25x2 //= -57121x = (x1+x2) // ihe += xpai = he*4pai //= 10**10end_time = time.time()run_time = str(end_time - start_time)paistr=str(pai)paistr=paistr[:1] + '.' + paistr[1:]f=open('pai.txt','w')f.write(paistr)f.close()print ('运行时间:' + run_time )#print ('计算结果:',pai)print ('\n'*1) comput()
comput()

  代码github地址:https://github.com/JackieZheng/Python-Pai

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/285876.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【Pix4d精品教程】垂直摄影空三加密生成DOM和DSM,并按10m间距提取高程点,生成等高线

在航测项目中,通常使用pix4d软件跑空三,生成数字正射影像DOM和数字表面模型DSM,当然还有点云数据。那么,在此基础之上,怎样根据航测要求,生成等间距5m或者10m的高程点和等高线呢,如上图所示?本文就来借助ArcGIS软件完成此波操作。 文章目录 一、Pix4d生成DOM和DSM二、A…

介绍一款python类型检查工具pyright

介绍 近日,微软在 Github 上开源了一个 Python 静态类型检查工具:pyright ,引起了社区内的多方关注。 微软在开源项目上的参与力度是越来越大了,不说收购 Github 这种大的战略野心,只说它家开源的 VS Code 编辑器&am…

Eclipse利用Maven2搭建SpringMVC框架的Web工程

一、准备工作: 下载apache-maven--> 配置Maven_home -->下载Eclipse Maven插件 二、新建工程: 选择新建Maven Project archetype选择webapp-->输入group ID (src下包名)和Artifact ID (工程名) 新建Maven工程目录如上图…

【ArcGIS风暴】ArcGIS10.6获取栅格影像边界范围的三种方法案例详解

基于ArcGIS平台有多种办法可以提取栅格影像边界,常见的方法有3种: 栅格范围(Raster Domain)栅格转面(Raster to Polygon)创建轮廓(BuildFootprints)/构建边界(BuildBoundary)文章目录 1. 栅格范围(Raster Domain)2. 栅格转面(Raster to Polygon)3. 创建轮廓(Bui…

IDEA中使用数据库可视化操作工具

文章目录 1.入门介绍2. 没有数据库驱动3. 准备&测试连接3.1测试报错 4.连接5.编写SQL语句 1.入门介绍 在IDEA的专业版的右侧工具栏应该会有DataBase按钮如果没有的同学可以这样操作(必须是IDEA专业版) 新建数据库 2. 没有数据库驱动 如果提示: missing driver files ,…

WPF效果第一百九十篇之再耍ListBox

前面一篇效果基于Expander和ListBox实现了一下所需要的效果;今天再次实现点底部不一样的效果;最终实现的效果:1、ItemContainerStyle我是比较简单粗暴直接分了二行:ListBoxCanvas实现:<ControlTemplate TargetType"{x:Type ListBoxItem}"><Grid Background&…

Rabbitmq~对Vhost的配置

rabbitmq里有一些概念我们要清楚&#xff0c;如vhost,channel,exchange,queue等&#xff0c;而前段时间在部署rabbitmq环境时启用了虚拟主机vhost&#xff0c;感觉他主要是起到了消息隔离的作用,下面分别再说一下它们的知识。 VHost vhost去做第一层的区分&#xff0c;虚拟主机…

表单元素 开篇

今天开始讲述表单这个重要模块 可以说,JS 最早是为表单而发明的, 因此在没有JS之前,所有操作都需要提交后端验证,发现有误再重定向回原页面, 加上之前1,2KB的网速,这用户体验真是奇差无比.因此JS最初发明出来&#xff0c;就是做表单验证的&#xff0e; 围绕表单&#xff0c;添加…

如何画出一张合格的技术架构图?

阿里妹导读&#xff1a;技术传播的价值&#xff0c;不仅仅体现在通过商业化产品和开源项目来缩短我们构建应用的路径&#xff0c;加速业务的上线速率&#xff0c;也体现在优秀工程师在工作效率提升、产品性能优化和用户体验改善等经验方面的分享&#xff0c;以提高我们的专业能…

.NET 发布和支持计划介绍

.NET 发布和支持计划介绍Intro对于 .NET 的发布&#xff0c;大多数童鞋都知道现在每年发布一个版本&#xff0c;针对 .NET 的发布&#xff0c;最近有些更新&#xff0c;Current 版本将改为 STS 版本&#xff0c;所以写一篇文章介绍一下每年 11 月都会发布新的 .NET 主要版本&am…

【spring boot】8.spring boot的日志框架logback使用

在继续上一篇的Debug调试之后&#xff0c;把spring boot的日志框架使用情况逐步蚕食。 参考&#xff1a;http://tengj.top/2017/04/05/springbo 开篇之前&#xff0c;贴上完整application.properties日志相关配置 简介&#xff1a;spring boot的默认日志框架Logback SLF4J——…

通过iscsi配置在aix上挂载存储设备

本文中我们利用starwind虚拟存储进行设置&#xff0c;以下为实验环境说明&#xff1a;Windows环境&#xff1a;win7&#xff0c;ip address:10.3.5.7&#xff0c;iscsi initiator name &#xff1a;iqn.2008-08.com.starwindsoftware:joker-pc-aixAix环境&#xff1a;ip addres…

原生js声音播放代码

最终测试页页面 测试页面html代码(test.html) <!doctype html> <html lang"en"><head><meta charset"UTF-8"><meta name"Generator" content"EditPlus"><meta name"Author" content"…

写给 Kubernetes 工程师的 mTLS 指南

本文翻译节选自 A Kubernetes engineer’s guide to mTLS[1]&#xff0c;为了便于读者理解&#xff0c;笔者对原文做了一点修改 &#xff08;本文删除了原文中的与主题关系不大的 Linkerd 安装的部分&#xff0c;将 Twillio 替换成国内读者比较熟悉的 GitHub&#xff09;。因为…

二、solidity 基础进阶(2.1)—— library 库合约《实战NFT web3 solidity(新版本0.8.+)》

《web3 solidity0.8.版本&#xff08;持续更新新版本内容&#xff09; 基础到实战NFT开发》会及时更新新版本 solidity 内容&#xff0c;以及完成最终的 NFT 实战商业项目部分。 注&#xff1a;由于是付费专栏内容&#xff0c;若有错误请及时联系1_bit&#xff0c;博客链接&am…

paper 17 : 机器学习算法思想简单梳理

前言&#xff1a; 本文总结的常见机器学习算法&#xff08;主要是一些常规分类器&#xff09;大概流程和主要思想。 朴素贝叶斯&#xff1a; 有以下几个地方需要注意&#xff1a; 1. 如果给出的特征向量长度可能不同&#xff0c;这是需要归一化为通长度的向量&#xff08;这里以…

BZOJ1179 Atm //缩点+spfa

1179: [Apio2009]Atm Description Input 第一行包含两个整数N、M。N表示路口的个数&#xff0c;M表示道路条数。接下来M行&#xff0c;每行两个整数&#xff0c;这两个整数都在1到N之间&#xff0c;第i1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行&#xff0c…

基于Spring Boot和Spring Cloud实现微服务架构学习

目录 Spring 顶级框架 Spring cloud子项目 WHAT - 什么是微服务 微服务简介 微服务的具体特征 SOA vs Microservice HOW - 怎么具体实践微服务 客户端如何访问这些服务&#xff1f; 服务之间如何通信&#xff1f; 这么多服务&#xff0c;怎么找? 这么多服务&#x…

ArcGIS实验教程——实验三十七:基于ArcGIS的太阳辐射分析案例教程

ArcGIS实验视频教程合集:《ArcGIS实验教程从入门到精通》(附配套实验数据)》 文章目录 一、太阳辐射的基本概念1. 视域2. 太阳图3. 星空图二、太阳辐射ArcGIS案例实现1. 对该区域进行太阳辐射区域分析2. 对单个点的太阳辐射进行分析太阳辐射是地球上各种物理过程和生物过程的…

ArcGIS实验教程——实验三十八:基于ArcGIS的等高线、山体阴影、山顶点提取案例教程

ArcGIS实验视频教程合集:《ArcGIS实验教程从入门到精通》(附配套实验数据)》 文章目录 1. 加载DEM2. 提取等高距为15m的等高线3. 提取等高距为75m的等高线4. 生成山体阴影5. 生成三维等高线6. 提取山顶点7. 实验数据下载地址山顶点指那些在特定邻域分析范围内,该点都比周围…