[数值计算] 批处理计算圆周率任意位

Rank: 8 Rank: 8

帖子: 640
积分: 2838
技术: 101
捐助: 0
注册时间: 2008-1-17

1楼 跳转到 » 倒序看帖

打印

字体大小: tT

发表于 2009-2-22 21:50 | 只看该作者

[数值计算] 批处理计算圆周率任意位

本帖最后由 plp626 于 2014-5-16 23:11 编辑

直接双击:默认计算100位,
命令行下:第一个参数(正整数)为计算圆周率的位数

 
@echo off&title 计算圆周率 by plp626
setlocal enabledelayedexpansion
if %1.==. (set c=100) else set/a c=%1
set/a c=(c*100/3)+70,cc=c/10,n=cc
for /l %%a in (1 1 %cc%)do set/a f_%%a=2000
for /l %%a in (%c% -132 100)do (set/a n=%%a/10,m=2*n-1
   set/a "d=f_!n!*10000,f_!n!=d%%m,d=d/m,n-=1"
   for /l %%b in (!n! -1 1)do (set/a n=%%b,m=2*n-1
      set/a "d=d*n+f_!n!*10000,f_!n!=d%%m,d=d/m,n-=1"
   )
   set/a "an=e+d/10000,e=d%%10000"
   if !an! lss 1000 set an=000!an!
   set/p=!an:~-4!<nul
)
echo\&echo/&set/p=数学爱好者加我QQ:275928264，验证号:plp626
复制代码

所用公式:
pi=2*1+2*1!/3!!+2*2!/5!!+2*3!/7!!+...+2*k!/(2*k+1)!!+...

6 评分人数

老刘1号: 6666666666666666666666666666666666666666技术 + 1
523066680: 今天回来看帖，膜拜！PB + 30 技术 + 10
zm900612: 保持队形技术 + 1
随风: 牛PB + 11
batman: 数学高人，多谢分享！技术 + 1

yslyxqysl

荣誉版主

Rank: 8 Rank: 8

帖子: 342
积分: 3936
技术: 9
捐助: 0
注册时间: 2008-8-27

2楼

发表于 2009-2-22 21:56 | 只看该作者

最后少了个<nul&pause.

[ 本帖最后由 yslyxqysl 于 2009-2-22 21:58 编辑 ]

Rank: 8 Rank: 8

帖子: 640
积分: 2838
技术: 101
捐助: 0
注册时间: 2008-1-17

3楼

发表于 2009-2-22 22:00 | 只看该作者

回复 2楼的帖子

no,no,no
不需要....

lockmove

六级士官

Rank: 4

帖子: 45
积分: 377
技术: 0
捐助: 0
注册时间: 2009-1-17

4楼

发表于 2009-2-22 22:39 | 只看该作者

有 set/p= 就停了吧.
pi=2*1+2*1!/3!!+2*2!/5!!+2*3!/7!!+...+2*k!/(2*k+1)!!+...

这公式在后面越来越慢,计算也有限,不过挺好了.

何时才有自己的头像?

batman

荣誉版主

Rank: 8 Rank: 8

帖子: 2702
积分: 14562
技术: 177
捐助: 630
注册时间: 2008-6-16

5楼

发表于 2009-2-23 08:17 | 只看该作者

plp626，兄弟，我们有好久没见了，我是cn-dos的ZW750516，你能来到我们论坛，我真是太高兴了。。。

＊＊＊共同提高＊＊＊

Rank: 8 Rank: 8

帖子: 640
积分: 2838
技术: 101
捐助: 0
注册时间: 2008-1-17

6楼

发表于 2009-2-23 17:01 | 只看该作者

zw兄，还记得俺，感觉很亲切，其实我早来这里了，只是没时间上网罢了，

BBCC

中尉

Rank: 5 Rank: 5

帖子: 245
积分: 914
技术: 1
捐助: 0
注册时间: 2008-10-25

7楼

发表于 2009-2-23 18:47 | 只看该作者

昨晚正想问你公式是啥的...
之前见est用python写过一个,好像速度更快...

for /f "delims=" %%a in ('%0') do (echo %%a)

Batcher

管理员

Rank: 12

帖子: 14935
积分: 46142
技术: 857
捐助: 745
注册时间: 2008-6-9

8楼

发表于 2009-2-23 19:57 | 只看该作者

回复 7楼的帖子

python的优势之一是它的数据类型可以存放N大的数据^_^

我帮忙写的代码不需要付钱。如果一定要给，请在微信群或QQ群发给大家吧。
【微信公众号、微信群、QQ群】http://bbs.bathome.net/thread-3473-1-1.html
【支持批处理之家，加入VIP会员！】http://bbs.bathome.net/thread-67716-1-1.html

zqz0012005

荣誉版主

Rank: 8 Rank: 8

帖子: 1069
积分: 12123
技术: 38
捐助: 0
注册时间: 2008-5-10

9楼

发表于 2009-2-24 18:16 | 只看该作者

plp626这个“老古董”居然重新冒出来了。

命令行参考：hh.exe ntcmds.chm::/ntcmds.htm
求助者请拿出诚心，别人才愿意奉献热心！
把查看手册形成条件反射！

523066680

版主

Rank: 7 Rank: 7 Rank: 7

帖子: 3151
积分: 6455
技术: 317
捐助: 70
注册时间: 2008-8-3

10楼

发表于 2009-2-26 16:59 | 只看该作者

核潜水艇浮起来拉～
话说这里成为迎宾大道了欢迎～

FuniCode 编程论坛

Rank: 8 Rank: 8

帖子: 640
积分: 2838
技术: 101
捐助: 0
注册时间: 2008-1-17

11楼

发表于 2011-4-25 11:34 | 只看该作者

一楼代码有一处笔误，导致后面精度有所丢失，两年多了，今儿才发现，现在已改正过来。

技术娱乐 3802

CrLf

论坛巡查

Rank: 8 Rank: 8

帖子: 6388
积分: 18843
技术: 982
捐助: 100
注册时间: 2010-10-9

12楼

发表于 2011-4-25 13:03 | 只看该作者

这个牛，顶一下古董兄

链接：在线第三方命令行工具下载 bat、vbs、js 原生混编

zz100001

六级士官

Rank: 4

帖子: 243
积分: 376
技术: 24
捐助: 50
注册时间: 2010-7-25

13楼

发表于 2011-4-25 14:25 | 只看该作者

为什么叫“古董兄”？

▄︻┻═┳一

列兵

Rank: 1

帖子: 29
积分: 37
技术: 0
捐助: 0
注册时间: 2011-4-24

14楼

发表于 2011-4-25 14:28 | 只看该作者

13# zz100001

你看他的注册日期就明白了。

applba

中校

Rank: 6 Rank: 6

帖子: 551
积分: 2799
技术: 39
捐助: 0
注册时间: 2011-4-22

15楼

发表于 2011-4-25 15:19 | 只看该作者

关键是效率……