#include <iostream>

using namespace std;

long a=10000,b=0,c=2800,d,e=0,f[2801],g;

void main()

{

for(;b-c;)f[b++]=a/5;

for(;d=0,g=c*2;c-=14,cout<<e+d/a,e=d%a)

for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);

}

它可以输出π的前781位。

楼上这位大哥，你这样改能行么？

cout<<e+d/a,

人原本写的是printf("%.4d",e+d/a),

看懂算法再改，看不懂，改了白改。

当然能行了！knocker 你如不懂C++，就不要说，自己试试就知道了

想不到能把for循环搞成那样，请详细解释一下算法，谢谢！

哪位能说一下计算圆周率的算法啊，谢谢

楼主的程序很好呀，不过少了第32位，我背过PI的应该是50288
#include <iostream>
using namespace std;
long a=10000,b=0,c=2800,d,e=0,f[2801],g;
int main(){
while(b-c)
f[b++]=a/5;
for(;d=0,g=c*2;c-=14,cout<<e+d/a,e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
}
程序稍微改了一下，这个我敢肯定是从C语言移植过来的，把for用成这样只能是学C的人干的