大数阶乘的计算(三)

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

下面的算法在 http://www.csdn.net/Develop/Edit.asp?id=28308基础上改进,比其至少快10%:

Sub calcfactorial(ByVal n As Integer)

Dim XYS() As Integer, x() As Integer, y() As Integer, result() As String, i As Long, j As Long, k As Long, TEMP As Long, stimer As Double

If n < 0 Then Exit Sub

ReDim x(1)

ReDim XYS(1)

x(1) = 1

XYS(1) = 1

k = 1

stimer = Timer

Do While k <= n

TEMP = Len(CStr(k))

ReDim y(1 To TEMP)

For i = 1 To TEMP

y(i) = Val(Mid(k, TEMP + 1 - i, 1))

Next

ReDim XYS(1 To UBound(x) + UBound(y))

For i = 1 To UBound(x)

For j = 1 To UBound(y)

XYS(i + j - 1) = XYS(i + j - 1) + x(i) * y(j)

Next

Next

For i = 1 To UBound(x) + UBound(y) - 1

TEMP = XYS(i) \ 10

XYS(i) = XYS(i) Mod 10

XYS(i + 1) = XYS(i + 1) + TEMP

Next

x = XYS

If x(UBound(x)) = 0 Then ReDim Preserve x(1 To UBound(x) - 1)

k = k + 1

Loop

ReDim result(1 To UBound(x))

For i = 1 To UBound(x)

result(i) = x(UBound(x) + 1 - i)

Next

factorial = Join(result, "")

Debug.Print k - 1 & "! : 用时 "; Timer - stimer & " 秒, 结果 " & UBound(x) & " 位"

'Debug.Print factorial

Erase x()

Erase y()

Erase XYS()

Erase result()

End Sub

Private Sub Command1_Click()

For i = 1 To 9

calcfactorial i * 100

Next

For i = 1 To 10

calcfactorial i * 100

Next

End Sub

结果输出:

100! : 用时 1.17187501018634E-05 秒, 结果 158 位

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

300! : 用时 .1090087890625 秒, 结果 615 位

400! : 用时 .219002929687576 秒, 结果 869 位

500! : 用时 .344011718749925 秒, 结果 1135 位

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

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

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

900! : 用时 1.281005859375 秒, 结果 2270 位

1000! : 用时 1.59400292968758 秒, 结果 2568 位

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

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

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

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

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

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

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

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

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

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