"1 million factorial"의 두 판 사이의 차이

2017년 6월 19일 (월) 12:19 기준 최신판

ril https://www.quora.com/What-is-an-efficient-algorithm-to-find-the-factorial-of-huge-numbers-which-lie-in-the-range-100000-1000000

just by using bc

rewrote a factorial function not to use recursive calling.

define g(x) {
    answer=1;
    for(i=1;i<x+1;i++) {
        answer *= i;
    }
    return answer;
}

g(1000000);

real 578m28.905s
user 450m56.192s
sys 12m33.006s
(with poor computing power (centOS virtual machine on windows7. i5 processor))

5,565,709 digits starting

82639316883312400623766461031726662911353479789638730451677758855633796110356450844465305113114639733516068042108785885414647469506478361823012109754232995901156417462491737988838926919341417654578323931987280247219893964365444552161533920583519938798941774206240841593987701818807223169252057737128436859815222389311521255279546829742282164292748493887784712443572285950934362117645254493052265841197629905619012120241419002534128319433065076207004051595915117186613844750900755834037427137686877042093751023502633401248341314910217684549431273636399066971952961345733318557782792616690299056202054369409707066647851950401003675381978549679950259346666425613978573559764142083506…

ending with 249,998 zeros(about 4.49% of the whole digits)

more easy way to get the number of trailing zeros(from here)

Every multiple of 5 contributes a zero.
Every multiple of 25 contributes a second zero
Every multiple of 125 contributes a third zero
Every multiple of 625 contributes a fourth zero
etc.

floor(1000000/5) + floor(1000000/25) + ... floor(1000000/5^8)
= 200,000 + 40,000 + 8,000 +1,600 +320 + 64 +12 + 2
= 249,998 trailing zeros

someone can get this in about only 10sec!! here
He's saying he used FFT

ref

http://www.luschny.de/math/factorial/FastFactorialFunctions.htm

http://www.luschny.de/math/factorial/csharp/FactorialPrimeSwing.cs.html

http://answers.google.com/answers/threadview/id/509662.html

http://en.wikipedia.org/wiki/Elliptic_curve_factorization

@@ 1번째 줄: / 1번째 줄: @@
-I'm interested in this position <a href=" http://atameqyrasy.blog4ever.com/blog/lire-article-665231-8666974-little_tits_tgp.html ">tyeen little video</a>  nxy <a href=" http://dumesoqatiq.blog4ever.com/blog/lire-article-665218-8666963-little_doll.html ">little boy butts</a>  68974 <a href=" http://eibeqisynusi.blog4ever.com/blog/lire-article-665243-8666980-little_teen_movies.html ">little tiny nudist</a>  71034 <a href=" http://nyheteqeeqa.blog4ever.com/blog/lire-article-665249-8666988-little_porn_girl.html ">little pusy xxx</a>  4716 <a href=" http://ikehuyjufa.blog4ever.com/blog/lire-article-665213-8666956-naughty_little_teens.html ">little to fuck</a>  hxmgsf <a href=" http://asofeciqukon.blog4ever.com/blog/lire-article-665248-8666986-little_angels_upskirts.html ">sex little</a>  gygem <a href=" http://kibomagaho.blog4ever.com/blog/lire-article-665210-8666955-little_nude_illegal.html ">little angles nude</a>  8DDD <a href=" http://fedenifibula.blog4ever.com/blog/lire-article-665225-8666967-illegal_little_sluts.html ">little art porn</a>  412210 <a href=" http://uhyfatefana.blog4ever.com/blog/lire-article-665205-8666950-little_kate.html ">little boys masturbating</a>  qwj <a href=" http://ukydejeou.blog4ever.com/blog/lire-article-665228-8666971-innocent_little_boys.html ">little monica hentai</a>  461
+#ril https://www.quora.com/What-is-an-efficient-algorithm-to-find-the-factorial-of-huge-numbers-which-lie-in-the-range-100000-1000000
+== just by using [http://www.manpagez.com/man/1/bc/ bc] ==
+rewrote a factorial function not to use recursive calling.
+<syntaxhighlight lang="c">
+define g(x) {
+    answer=1;
+    for(i=1;i<x+1;i++) {
+        answer *= i;
+    }
+    return answer;
+}
+g(1000000);
+</syntaxhighlight>
+real    578m28.905s<br>
+user    450m56.192s<br>
+sys     12m33.006s<br>
+(with poor computing power (centOS virtual machine on windows7. i5 processor))
+,565,709 digits starting <div style='word-wrap: break-word;'>82639316883312400623766461031726662911353479789638730451677758855633796110356450844465305113114639733516068042108785885414647469506478361823012109754232995901156417462491737988838926919341417654578323931987280247219893964365444552161533920583519938798941774206240841593987701818807223169252057737128436859815222389311521255279546829742282164292748493887784712443572285950934362117645254493052265841197629905619012120241419002534128319433065076207004051595915117186613844750900755834037427137686877042093751023502633401248341314910217684549431273636399066971952961345733318557782792616690299056202054369409707066647851950401003675381978549679950259346666425613978573559764142083506&hellip; </div>ending with 249,998 zeros(about 4.49% of the whole digits)
+* more easy way to get the number of trailing zeros(from [http://answers.yahoo.com/question/index?qid=20071017145425AAaGkd2 here])
+<blockquote>
+Every multiple of 5 contributes a zero.<br>
+Every multiple of 25 contributes a second zero<br>
+Every multiple of 125 contributes a third zero<br>
+Every multiple of 625 contributes a fourth zero<br>
+etc.<br>
+<br>
+floor(1000000/5) + floor(1000000/25) + ... floor(1000000/5^8)<br>
+= 200,000 + 40,000 + 8,000 +1,600 +320 + 64 +12 + 2 <br>
+= 249,998 trailing zeros
+</blockquote>
+* someone can get this in about only 10sec!! [http://www.makewebgames.com/showthread.php/24362-Factorial-Program-in-C-1-Million-Factorial here]<br>He's saying he used FFT
+== ref ==
+http://www.luschny.de/math/factorial/FastFactorialFunctions.htm
+http://www.luschny.de/math/factorial/csharp/FactorialPrimeSwing.cs.html
+http://answers.google.com/answers/threadview/id/509662.html
+http://en.wikipedia.org/wiki/Elliptic_curve_factorization