作者DrStein (啤酒肚)
看板C_and_CPP
標題Re: [問題] 有人知道這pi公式的名字嗎?
時間Sat Dec 4 19:16:31 2010
※ 引述《tropical72 (藍影)》之銘言:
: 看到這個不禁想到目前史上最強的 pi... 三行算千位
: 這程序是對岸寫的
: #include <stdio.h>
: int main()
: {
: long a=10000,b=0,c=2800,d=0,e=0,f[2801],g=0;
: for(;b-c;) f[b++]=a/5;
: for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
: for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
: return 0;
: }
: 到現在還是完全看不懂它是怎麼算的.
: 不過我知道它是用 J.Marchin 原理下去跑的,
: 網路上有人示範 excel 怎麼拉 1000 位的 pi,
: 和這份原始碼跑的方式很像..
: http://www.esproc.com/exercises/pi1000.html
: 不禁又要長嘆, 對岸強者怎麼這麼多..
題外話, <math.h>中的函數
是用什麼公式去算(ex:sin, cos, atan, log , pow)的啊?
感覺很快又精確。
有人知道的嗎?
或是有沒有這方面的資料可以查哩
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 58.115.143.194
推 purpose:基本上,你自己反組譯看指令就知道大概了 12/04 20:07
推 tropical72:內建的函數我覺得似乎沒很快耶..fabs隨便寫都比它快吧. 12/04 23:21
推 tropical72:我覺得可以去找找power series, nonlinear root 方面. 12/04 23:51