|
阅读：4684回复：1
大整数的阶乘计算
<P>n!=1*2*3*....*n → lg(n!)=lg(1*2*3*....*n)=lg(1)+lg(2)+ lg(3)+..+lg(n) → n!=10^(lg(1)+lg(2)+ lg(3)+..+lg(n))</P>
<P>考虑到当n很大时，输出具体的结果并无多大实际意义，遂另写代码如下：</P> <P>Sub calcfactorial(ByVal num As Long, Optional ByRef POWER As String, Optional ByRef FIRSTNUM As String)<BR>Dim I As Long, TEMP As Double, temp2 As Long, STIME As Single<BR>TEMP = 0<BR>temp2 = 0<BR>STIME = Timer<BR>For I = 1 To num<BR>TEMP = TEMP + Log(I) / Log(10)<BR>If TEMP > 1000000 Then<BR>temp2 = temp2 + 1<BR>TEMP = TEMP - 1000000<BR>End If<BR>Next<BR>POWER = Format(Int(TEMP + 1), "000000")<BR>TEMP = TEMP - Val(POWER)<BR>POWER = temp2 ; POWER<BR>If Val(POWER) < 10 ^ 6 Then POWER = Val(POWER)<BR>FIRSTNUM = Left(Replace(10 ^ (TEMP), ".", ""), 10)<BR>Debug.Print Right(Space(9) ; num, 9) ; "！ ： 用时 " ; Right(Space(8) ; Format(Timer - STIME, "0.00000"), 8) ; " 秒, 结果 " ; Right(Space(10) ; POWER, 10) ; " 位,前10位为 " ; FIRSTNUM<BR>End Sub</P> <P><BR>Private Sub Command1_Click()<BR>Dim I As Long, J As Long<BR>For J = 2 To 7<BR>For I = 1 To 10<BR>calcfactorial I * 10 ^ J<BR>Next<BR>Debug.Print<BR>Next<BR>End Sub</P> <br> <P>输出结果：</P> <P>' 100！ ： 用时 0.00000 秒, 结果 158 位,前10位为 9332621544<BR>' 200！ ： 用时 0.00000 秒, 结果 375 位,前10位为 7886578673<BR>' 300！ ： 用时 0.00000 秒, 结果 615 位,前10位为 3060575122<BR>' 400！ ： 用时 0.00000 秒, 结果 869 位,前10位为 6403452284<BR>' 500！ ： 用时 0.00000 秒, 结果 1135 位,前10位为 1220136825<BR>' 600！ ： 用时 0.00000 秒, 结果 1409 位,前10位为 1265572316<BR>' 700！ ： 用时 0.01563 秒, 结果 1690 位,前10位为 2422040124<BR>' 800！ ： 用时 0.00000 秒, 结果 1977 位,前10位为 7710530113<BR>' 900！ ： 用时 0.00000 秒, 结果 2270 位,前10位为 6752680220<BR>' 1000！ ： 用时 0.00000 秒, 结果 2568 位,前10位为 4023872600<BR>'<BR>' 1000！ ： 用时 0.00000 秒, 结果 2568 位,前10位为 4023872600<BR>' 2000！ ： 用时 0.01563 秒, 结果 5736 位,前10位为 3316275092<BR>' 3000！ ： 用时 0.00000 秒, 结果 9131 位,前10位为 4149359603<BR>' 4000！ ： 用时 0.00000 秒, 结果 12674 位,前10位为 1828801951<BR>' 5000！ ： 用时 0.00000 秒, 结果 16326 位,前10位为 4228577926<BR>' 6000！ ： 用时 0.01563 秒, 结果 20066 位,前10位为 2683999765<BR>' 7000！ ： 用时 0.00000 秒, 结果 23878 位,前10位为 8842007956<BR>' 8000！ ： 用时 0.00000 秒, 结果 27753 位,前10位为 5184181060<BR>' 9000！ ： 用时 0.01563 秒, 结果 31682 位,前10位为 8099589986<BR>' 10000！ ： 用时 0.00000 秒, 结果 35660 位,前10位为 2846259680<BR>'<BR>' 10000！ ： 用时 0.00000 秒, 结果 35660 位,前10位为 2846259680<BR>' 20000！ ： 用时 0.00000 秒, 结果 77338 位,前10位为 1819206320<BR>' 30000！ ： 用时 0.04688 秒, 结果 121288 位,前10位为 2759537248<BR>' 40000！ ： 用时 0.01563 秒, 结果 166714 位,前10位为 2091692420<BR>' 50000！ ： 用时 0.03125 秒, 结果 213237 位,前10位为 3347320501<BR>' 60000！ ： 用时 0.01563 秒, 结果 260635 位,前10位为 1564137704<BR>' 70000！ ： 用时 0.03125 秒, 结果 308760 位,前10位为 1176812415<BR>' 80000！ ： 用时 0.03125 秒, 结果 357507 位,前10位为 3097722262<BR>' 90000！ ： 用时 0.03125 秒, 结果 406799 位,前10位为 1580119155<BR>' 100000！ ： 用时 0.04688 秒, 结果 456574 位,前10位为 2824229408<BR>'<BR>' 100000！ ： 用时 0.04688 秒, 结果 456574 位,前10位为 2824229408<BR>' 200000！ ： 用时 0.10938 秒, 结果 973351 位,前10位为 1420225354<BR>' 300000！ ： 用时 0.12500 秒, 结果 1512852 位,前10位为 1477391532<BR>' 400000！ ： 用时 0.23438 秒, 结果 2067110 位,前10位为 2534486091<BR>' 500000！ ： 用时 0.21875 秒, 结果 2632342 位,前10位为 1022801596<BR>' 600000！ ： 用时 0.35938 秒, 结果 3206318 位,前10位为 2234878214<BR>' 700000！ ： 用时 0.40625 秒, 结果 3787566 位,前10位为 6491076363<BR>' 800000！ ： 用时 0.43750 秒, 结果 4375040 位,前10位为 5684678869<BR>' 900000！ ： 用时 0.54688 秒, 结果 4967957 位,前10位为 3990191306<BR>' 1000000！ ： 用时 0.54688 秒, 结果 5565709 位,前10位为 8263932189<BR>'<BR>' 1000000！ ： 用时 0.51563 秒, 结果 5565709 位,前10位为 8263932189<BR>' 2000000！ ： 用时 1.10938 秒, 结果 11733475 位,前10位为 3776821110<BR>' 3000000！ ： 用时 1.56250 秒, 结果 18128484 位,前10位为 9038572324<BR>' 4000000！ ： 用时 1.93750 秒, 结果 24671066 位,前10位为 5467878127<BR>' 5000000！ ： 用时 2.07813 秒, 结果 31323382 位,前10位为 2294768200<BR>' 6000000！ ： 用时 2.53125 秒, 结果 38063145 位,前10位为 2506385954<BR>' 7000000！ ： 用时 2.90625 秒, 结果 44875629 位,前10位为 5350767028<BR>' 8000000！ ： 用时 3.35938 秒, 结果 51750368 位,前10位为 7786539372<BR>' 9000000！ ： 用时 3.75000 秒, 结果 58679537 位,前10位为 1330564471<BR>' 10000000！ ： 用时 4.20313 秒, 结果 65657060 位,前10位为 1202423617<BR>'<BR>' 10000000！ ： 用时 4.17188 秒, 结果 65657060 位,前10位为 1202423617<BR>' 20000000！ ： 用时 8.46875 秒, 结果 137334715 位,前10位为 2112611265<BR>' 30000000！ ： 用时 14.59375 秒, 结果 211284808 位,前10位为 2099630957<BR>' 40000000！ ： 用时 16.75000 秒, 结果 286710625 位,前10位为 3776665124<BR>' 50000000！ ： 用时 23.31250 秒, 结果 363233781 位,前10位为 2345380960<BR>' 60000000！ ： 用时 27.28125 秒, 结果 440631411 位,前10位为 2494533092<BR>' 70000000！ ： 用时 30.10938 秒, 结果 518756254 位,前10位为 2451374102<BR>' 80000000！ ： 用时 34.85938 秒, 结果 597503645 位,前10位为 5724408058<BR>' 90000000！ ： 用时 37.65625 秒, 结果 676795327 位,前10位为 7152504088<BR>'100000000！ ： 用时 46.21875 秒, 结果 756570557 位,前10位为 1617205312<BR></P> |
|
|
|
1楼#
发布于：2007-11-14 00:50
<P>好强啊，我以前听说有用数组实现的，没想到还有这样的算法，多谢了！<img src="images/post/smile/dvbbs/em02.gif" /></P>
|
|
|
