大数阶乘的计算(四)

王朝vb·作者佚名  2006-01-09
窄屏简体版  字體: |||超大  

看了northwolves 的文章大数阶乘的计算(三) http://dev.csdn.net/article/28/28432.shtm 很有启发,联想起以前自己编写过一个类似的函数,方法与其差不多,都是动态数组的增位保存方式。效率上比northwolves兄要差一些,大约相差50%;当时以为是最快的阶乘算法了,所以看到这篇文章很是钦佩!

不过,我总觉得,这个算法还有可优化的余地,于是,又重新做了点新尝试,不想大获成功,效率比大数阶乘的计算(三)一下子又提高了5倍多,而且算法简单得惊人,故不敢独吞,赶紧贴出与大家共享。

Option Explicit

Private Function cacl(num As Long) As String

Dim numlen As Long, last As Long, x As Long

Dim i As Long, m As Long, n As Long

Dim result() As Long, starttime As Single

numlen = 1

starttime = Timer

ReDim result(1 To numlen)

result(1) = 1

x = 1

Do While x <= num

last = 0

For i = 1 To numlen

m = result(i) * x + last

result(i) = m Mod 10

last = m \ 10

Next

If last > 0 Then

n = Len(CStr(last))

ReDim Preserve result(1 To numlen + n)

For i = 1 To n

result(numlen + i) = last Mod 10

last = last \ 10

Next

numlen = UBound(result)

End If

x = x + 1

Loop

ReDim s(1 To numlen)

For i = 1 To numlen

s(i) = result(numlen + 1 - i)

Next

cacl = Join(s, "")

Debug.Print num & "! : 用时 "; Timer - starttime & " 秒, 结果 " & numlen & " 位"

End Function

Private Sub Command1_Click()

Dim i As Long

For i = 1 To 10

cacl i * 100

Next

'For i = 1 To 10

'cacl i * 1000

'Next

End Sub

输出结果:(我的机器在运行“计算(三)”时,与其结果相当,所以应是与northwolves兄,同档次的机器)

100! : 用时 0 秒, 结果 158 位

200! : 用时 1.171875E-02 秒, 结果 375 位

300! : 用时 3.515625E-02 秒, 结果 615 位

400! : 用时 4.296875E-02 秒, 结果 869 位

500! : 用时 8.203125E-02 秒, 结果 1135 位

600! : 用时 .1054688 秒, 结果 1409 位

700! : 用时 .15625 秒, 结果 1690 位

800! : 用时 .1992188 秒, 结果 1977 位

900! : 用时 .2617188 秒, 结果 2270 位

1000! : 用时 .3632813 秒, 结果 2568 位

2000! : 用时 1.53125 秒, 结果 5736 位

3000! : 用时 3.75 秒, 结果 9131 位

4000! : 用时 7.171875 秒, 结果 12674 位

5000! : 用时 11.45313 秒, 结果 16326 位

6000! : 用时 16.41016 秒, 结果 20066 位

7000! : 用时 22.83984 秒, 结果 23878 位

8000! : 用时 30.375 秒, 结果 27753 位

9000! : 用时 41.75 秒, 结果 31682 位

10000! : 用时 54.57813 秒, 结果 35660 位

这次输出慢了些,因为我在第一次测试10000!时,只用了47秒。不知会不会还有更快的算法。

 
 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
 
 
© 2005- 王朝網路 版權所有 導航