批处理计算圆周率任意位

plp626

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

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

复制代码

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

yslyxqysl

最后少了个<nul&pause.

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

plp626

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

lockmove

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

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

batman

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

plp626

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

BBCC

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

Batcher

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

zqz0012005

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

523066680

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

plp626

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

CrLf

这个牛，顶一下古董兄

zz100001

为什么叫“古董兄”？

▄︻┻═┳一

13# zz100001


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

applba

关键是效率……