中国DOS联盟

给大家出了个题目^_^，比较难，

没有任何实用价值，只为提供挑战平台。

[ Last edited by plp626 on 2009-11-30 at 14:33 ]

好久没来现在附代码如下：此贴已结！
@echo off&title 计算圆周率 by plp626@cn-dos.net setlocal enabledelayedexpansion if not %1.==. (set c=%1) else set c=100 set/a c=(c*100/3)+70,cc=c/10,count=0 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_!b!=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:************，

[ Last edited by plp626 on 2010-10-5 at 15:31 ]

说实话，我不知道圆周率是怎样计算的^_^

呵呵，我也不知道。

其实，没有所谓的纯批处理，所谓的纯只是限定用系统内置的命令（包括外部命令）。所以，只要将相关的外部命令（例如汇编程序）用批处理写入磁盘，然后调用，也就是纯批处理。这个应该不是很难吧？

[ Last edited by Climbing on 2008-4-30 at 02:12 PM ]

......

[ Last edited by plp626 on 2008-5-1 at 06:45 PM ]

Originally posted by Climbing at 2008-4-30 02:11 PM:
其实，没有所谓的纯批处理，所谓的纯只是限定用系统内置的命令（包括外部命令）。所以，只要将相关的外部命令（例如汇编程序）用批处理写入磁 ...

那个可能更复杂了，你不觉得吗？set位运行，/a参数不是现成的工具吗？

大学里的微积分知识全都还给老师了，看到就头晕。

我都不知道圆周率是什么东西，小的中学毕业，唉，惭愧

网上找到的一段VBS的圆周率计算说明

Atn 函数返回数值的反正切值。

Atn(number)

number 参数可以是任意有效的数值表达式。

说明
Atn 函数计算直角三角形两个边的比值 (number) 并返回对应角的弧度值。此比值是该角对边的长度与邻边长度之比。 结果的范围是从 -pi/2 到 pi/2 弧度。

弧度变换为角度的方法是将弧度乘以 pi/180。反之，角度变换为弧度的方法是将角度乘以180/pi 。

下面的示例利用 Atn 来计算 pi 的值:

Dim pipi = 4 * Atn(1) ' 计算 pi 的值。注意 Atn 是 Tan（将角作为参数返回直角三角形两边的比值）的反三角函数。不要混淆 Atn 与余切（正切的倒数 (1/tangent)）函数。

这是一段计算圆周率的VBS，不过很费时，我的机子花了30多秒才计算了6位（Abs(t) > 0.0000001，小数点后有多少个0，就精确到多少位）
k = 1: s = 1: t = 1: m = 1 time_ = timer While Abs(t) > 0.0000001 k = k + 2 m = -m t = m / k s = s + t Wend tt = timer - time_ MsgBox tt & vbCrLf & "圆周率π=" & 4 * s

C 语言的：
这个强悍，一瞬间（不到1秒）输出800位（已更正）：
#include < stdio.h> long a=10000, b, c=2800, d, e, f, g; main(){ for(;b-c;)f=a/5; for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a) for(b=c;d+=f*a,f=d%--g,d/=g--,--b;d*=b); }

原帖：http://wangcong.org/articles/puzzle.html

（上面代码和这个一样）这个大家能看懂（用的就是方案2）：
#include < stdio.h> long b, c=2800, d, e, f; main(){ while(b-c!=0){ f=2000; b++; } while(c!=0){ b=c; d=f*10000; f=d%(b*2-1); d=d/(b*2-1); --b; while(b!=0){ d=d*b+f*10000; f=d%(b*2-1); d=d/(b*2-1); --b; } c-=14; printf("%.4d",e+d/10000); e=d%10000; } }

要真想算100位小数就得有好算法，即使高级语言也不例外，C能提供的实型数据最高精度也才19位（查书），这个代码是把每步除法的的余数存放在数组里，保证不会有截断误差，
那这么看来，只要P能实现任意位数, 的两个整数除法，我们就可以用第二个方案算出Pi.
-------------------------------------------------
baidu批处理吧搜的分数化小数（或循环小数即任意精度），对分子分母不超过2^30的分数可行：

http://tieba.baidu.com/f?ct=335675392&tn=baiduPostBrowser&sc=1180921912&z=137039382&pn=0&rn=50&lm=0&word=%C5%FA%B4%A6%C0%ED#1180921912

咱论坛不得不爱版主有个算整数除法的代码，支持被除数任意位数，除数不超过8位数的(102楼)：

http://www.cn-dos.net/forum/viewthread.php?tid=23568&fpage=1&highlight=%E9%99%A4%E9%98%BF%E6%8B%89%E4%BC%AF%E6%95%B0%E5%AD%97%E5%A4%96&page=7

似乎快达到我们的要求？

[ Last edited by plp626 on 2008-4-30 at 05:04 PM ]

C 语言的：
这个强悍，一瞬间（不到1秒）输出2800位：


代码不知道你调试了没？反正看着是有错误。代码输出的是余数？貌似全是0
for的那部分有错误。。

如果对的话，还想改成VBS看看速度呢。



zh159的VBS不要用while，直接用for i=0 to 1000000....next

这个才3秒多就14位了。

改大点就好了。

不过……貌似还是14位……还是数据类型不好呀。

沒有看懂, 大學就沒有去學微積分~~~

Originally posted by slore at 2008-4-30 04:35 PM:
C 语言的：
这个强悍，一瞬间（不到1秒）输出2800位：


代码不知道你调试了没？反正看着是有错误。代码输出的是余数？貌似全是0
for的那部分有 ...

调试了，有源代码，exe文件在此。
http://upload.cn-dos.net/img/345.rar

你没加头文件#include <stdio.h>？

我粘贴有点匆忙，现在补上了。

[ Last edited by plp626 on 2008-4-30 at 05:16 PM ]

你看FOR语句就不对^b f[b]

不过找到了代码了。
#include <stdio.h>

long a=10000,b,c=2800,d,e,f[2801],g;
void main()
{
for(;b-c;)
f[b++]=a/5;
for(;d=0,g=c*2;c-=14,printf("%04d",e+d/a),e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
}

但是这个是(2800/14)*4位……