一个计算圆周率任意精度的spigot算法研究

plp626

38# liangbch


不知到为什么，你的代码在我的gcc上运行提示找不到lib-gmp10.dll，

精简版vc6.0编译报错误：
--------------------Configuration: p - Win32 Debug--------------------
Compiling...
p.c
d:\documents_and_settings\plp2\desktop\p.c(11) : warning C4244: 'function' : conversion from '__int64 ' to 'double ', possible loss of data
d:\documents_and_settings\plp2\desktop\p.c(14) : error C2063: 'DWORD' : not a function
d:\documents_and_settings\plp2\desktop\p.c(15) : error C2063: 'DWORD' : not a function
d:\documents_and_settings\plp2\desktop\p.c(26) : error C2059: syntax error : 'type'
d:\documents_and_settings\plp2\desktop\p.c(29) : error C2275: 'INT64' : illegal use of this type as an expression
        d:\documents_and_settings\plp2\desktop\p.c(6) : see declaration of 'INT64'
d:\documents_and_settings\plp2\desktop\p.c(29) : error C2146: syntax error : missing ';' before identifier 't2'
d:\documents_and_settings\plp2\desktop\p.c(29) : error C2065: 't2' : undeclared identifier
d:\documents_and_settings\plp2\desktop\p.c(29) : warning C4244: '=' : conversion from '__int64 ' to 'int ', possible loss of data
d:\documents_and_settings\plp2\desktop\p.c(32) : error C2065: 't1' : undeclared identifier
d:\documents_and_settings\plp2\desktop\p.c(32) : warning C4244: '/=' : conversion from '__int64 ' to 'int ', possible loss of data
d:\documents_and_settings\plp2\desktop\p.c(33) : warning C4244: '/=' : conversion from '__int64 ' to 'int ', possible loss of data
d:\documents_and_settings\plp2\desktop\p.c(39) : error C2059: syntax error : 'type'
d:\documents_and_settings\plp2\desktop\p.c(43) : error C2059: syntax error : 'else'
d:\documents_and_settings\plp2\desktop\p.c(47) : error C2065: 'low_p' : undeclared identifier
d:\documents_and_settings\plp2\desktop\p.c(47) : error C2065: 'high_p' : undeclared identifier
d:\documents_and_settings\plp2\desktop\p.c(47) : error C2099: initializer is not a constant
d:\documents_and_settings\plp2\desktop\p.c(48) : error C2065: 'low_q' : undeclared identifier
d:\documents_and_settings\plp2\desktop\p.c(48) : error C2065: 'high_q' : undeclared identifier
d:\documents_and_settings\plp2\desktop\p.c(48) : error C2099: initializer is not a constant
d:\documents_and_settings\plp2\desktop\p.c(49) : error C2059: syntax error : '}'

plp626

a,p作为分母超过40万(2^32/10000)了就可以舍弃了

plp626

这个问题要是能给出数学上的解法 该有多迷人啊。。。。

plp626

31# liangbch

sqrt(N)=(q/p)sqrt(1-1/a), 显然p,q是sqrt(N)的一个有理数逼近

这个思路确实不错。。

wayne

20# ▄︻┻═┳一‥

额，建议▄︻┻═┳一‥  童鞋 不要盗用我们老大的图像，。。。。

gxqcn

45# wayne

按 Ctrl+F5 强制刷新一下，头像才会是最新的。

wayne

46# gxqcn
呵呵，原来我这的信息落后了啊

wayne

41# plp626
第一，liangbch 兄的代码似乎并没有用到GMP库
第二，把print_pq(3,1000000I64); 里面的I64去掉，用 g++编译
第三，VC6.0 跟IE6一样，早就该进坟墓了

liangbch

48# wayne

  也许我比较保守吧。我一直用VC6.现在我的电脑中同时装有VC6和VC2008, 但VC6用的更多些。极少使用GCC编译自己的程序，尽管公司用的编译器主要是GCC。

wayne

49# liangbch
过谦了，你这叫做专业，