gis
gis
管理员
管理员
  • 注册日期2003-07-16
  • 发帖数15945
  • QQ554730525
  • 铜币25337枚
  • 威望15352点
  • 贡献值0点
  • 银元0个
  • GIS帝国居民
  • 帝国沙发管家
  • GIS帝国明星
  • GIS帝国铁杆
阅读:4441回复:1

大整数的阶乘计算

楼主#
更多 发布于:2005-08-29 14:48
<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>
喜欢0 评分0
jinyu1123
路人甲
路人甲
  • 注册日期2006-07-05
  • 发帖数67
  • QQ
  • 铜币262枚
  • 威望0点
  • 贡献值0点
  • 银元0个
1楼#
发布于:2007-11-14 00:50
<P>好强啊,我以前听说有用数组实现的,没想到还有这样的算法,多谢了!<img src="images/post/smile/dvbbs/em02.gif" /></P>
迷茫啊! Email : jinyu1123@163.com.
举报 回复(0) 喜欢(0)     评分
游客

返回顶部