外星人计算pi

x671 · #1

请问该代码是通过什么原理计算pi值的？

代码：全选

#!/usr/bin/python
def pi():
   a=10000
   c=50960
   e=0
   f=[a/5]*(c+1)
   while c>50000:
      d=0
      g=c*2
      b=c
      while 1:
         d+=f[b]*a
         g-=1
         d,f[b]=divmod(d,g)
         g-=1
         b-=1
         if b==0:break
         d*=b
      c-=14
      print str(e+d/a).zfill(4),
      e=d%a
   return
if __name__ == '__main__':
   pi()

$ ./pi.py
3141 5926 5358 9793 2384 6264 3383 2795 0288 4197 1693 9937 5105 8209 7494 4592 3078 1640 6286 2089 9862 8034 8253 4211 7067 9821 4808 6513 2823 0664 7093 8446 0955 0582 2317 2535 9408 1284 8111 7450 2841 0270 1938 5211 0555 9644 6229 4895 4930 3819 6442 8810 9756 6593 3446 1284 7564 8233 7867 8316 5271 2019 0914 5648 5669 2346 0348 6104 5432

stlxv · #2

虽然没看内容，不过直觉告诉我，原理就是数学上的某个

pi的展开式

forrid · #3

π的展开式？？？
看来数学是基础啊。。。

三翻领 · #4

${L_SMILIES_VERY_HAPPY}$
祖冲之看到这个，一口血喷死...........

ttand · #5

f(1)=1
f(n+1)=f(n)*n/(2*n+1)
PI/2=f(1)+f(2)+f(3)+.....+f(n)
n->无穷大。

PI={f(1)+f(2)+f(3)+.....+f(n)}*2

ttand · #6

c语言的

#include <stdio.h>
long a=10000,b,c=2800,d,e,f[2801],g;
int main(){
while(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*a,f=d%--g,d/=g--,--b;d*=b);
}

还有一个无法形容的bb
bt 的宏＋bt的算法

代码：全选

Here's another good one:

/*
 * Program to compute an approximation of pi
 * by Brian Westley, 1988
 * (requires pcc macro concatenation; try gcc -traditional-cpp)
 */

#define _ -F<00||--F-OO--;
int F=00,OO=00;
main(){F_OO();printf("%1.3f\n",4.*-F/OO/OO);}F_OO()
{
            _-_-_-_
       _-_-_-_-_-_-_-_-_
    _-_-_-_-_-_-_-_-_-_-_-_
  _-_-_-_-_-_-_-_-_-_-_-_-_-_
 _-_-_-_-_-_-_-_-_-_-_-_-_-_-_
 _-_-_-_-_-_-_-_-_-_-_-_-_-_-_
_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
 _-_-_-_-_-_-_-_-_-_-_-_-_-_-_
 _-_-_-_-_-_-_-_-_-_-_-_-_-_-_
  _-_-_-_-_-_-_-_-_-_-_-_-_-_
    _-_-_-_-_-_-_-_-_-_-_-_
        _-_-_-_-_-_-_-_
            _-_-_-_
}

cozo · #7

不是有一个一行代码算pi的程序的吗？比这个还外星人。