圆周率小数点后1000位为：
1415926535 8979323846 2643383279 5028841971 6939937510
5820974944 5923078164 0628620899 8628034825 3421170679
8214808651 3282306647 0938446095 5058223172 5359408128
4811174502 8410270193 8521105559 6446229489 5493038196
4428810975 6659334461 2847564823 3786783165 2712019091
4564856692 3460348610 4543266482 1339360726 0249141273
7245870066 0631558817 4881520920 9628292540 9171536436
7892590360 0113305305 4882046652 1384146951 9415116094
3305727036 5759591953 0921861173 8193261179 3105118548
0744623799 6274956735 1885752724 8912279381 8301194912
9833673362 4406566430 8602139494 6395224737 1907021798
6094370277 0539217176 2931767523 8467481846 7669405132
0005681271 4526356082 7785771342 7577896091 7363717872
1468440901 2249534301 4654958537 1050792279 6892589235
4201995611 2129021960 8640344181 5981362977 4771309960
5187072113 4999999837 2978049951 0597317328 1609631859
5024459455 3469083026 4252230825 3344685035 2619311881
7101000313 7838752886 5875332083 8142061717 7669147303
5982534904 2875546873 1159562863 8823537875 9375195778
1857780532 1712268066 1300192787 6611195909 2164201989

不是故弄玄虚，这应该就是所谓的混乱代码，国外有混乱代码大赛的

我还是看不懂，看懂的能不谁能讲解一下27楼的程序？

是用的这个公式吗？
[img]http://www.jason314.com/palgorithm.htm[/img]

27楼的仁兄
你的程序好像不行哟
编译都有问题
改了之后 结果全是0
(不知是不是我改错了)

34楼的仁兄
你的公式怎么无法显示
能不能再写一边

改写了一下，计算结果正确，不过依然看不明白。

#include<stdlib.h>
#include<stdio.h>

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

main(){
    int i;
    for(i=0;i<2801;i++)
        f[i]=2000;
    for(c=2800;c!=0;c-=14){
        d=0;
        g=c*2-1;
        b=c;
        while(b){
            d+=f[b]*a;
            f[b]=d%g;
            d/=g;
            g-=2;
            if(--b>0)
                d*=b;
        }
        printf("%.4d",e+d/a);
        e=d%a;
    }
}

看了一会儿，因为不长，可是看了无数遍，依然看不明白。

第一个循环里面把数组f[2801]里面每个数都赋予2000。

第二个循环有输出，而且看的出来是4位4位的输出，a的值一直是10000，也能看出来，输入的是e和d两个数，e就相当于汇编里的溢出，起进位的作用，d是输出的主要部分，并且把d/a的余数再赋予e。

循环里面的循环就是为输出这4个数做的计算，e没有改变，所以关键是看d的变化，循环里还有其他的变量c,f[],g,b，很乱啊，没有什么头绪。知道的是，c和b分别是为了历遍两个循环而用的index，并且一个递减14，一个递减1。那么剩下的只有g和f[]了，g所用到的大小一直是b*2-1。让人很容易联想到这个公式 pi/4=1-1/3+1/5-1/7+1/9……不过似乎又没什么联系。f[]也是很难考虑的东西，因为一直是在变化中（每个大循环修改一次，不过少修改14位）。最后看到头痛了还是不明白。

然后又想这个人是怎么想的呢？应该是刚开始写了很多、也写了很长，然后一次有一次改良，最后拿出了变态的来吓唬我们这些没有经验的人。我想如果可以，我们也应该自己先写一个，然后对照着再看他写的，不过似乎很难啊……

大家可以去这里看看：

http://www.cppfans.com/articles/basecalc/c_pi_10000.asp

http://www.apfloat.org/apfloat/
比小日本的 super pi 还要快的多