太郁闷了，python求质数怎么比VBA还慢！！

forceps

在我的机器上求40000以下的所有质数, python要61秒，Excel vba只要4秒多。

-------------------- python代码 --------------------
#!/usr/bin/env python
y=3
myfile=open('primes.txt','w')
while y<40000:
    x=y/2
    while x>1:
        if y%x==0:break
        x-=1
    else:
        print >>myfile,y,'\n'
    y+=1
myfile.close()
print 'End'
-------------------- End Python --------------------



-------------------- Excel VBA --------------------
Sub find_primer()
        Dim Arry1D_Primer(1 To 30000, 0 To 1)
        Dim y As Long, x As Long, m As Long, n As Long
        Dim i As Integer, k As Integer
        k = 0
        m = 1
        For y = 3 To 40000
                i = Int(y / 2)
                For x = i To 2 Step -1
                        If (y Mod x) = 0 Then
                                k = 1
                                Exit For
                        End If
                Next x
                If k = 0 Then
                        Arry1D_Primer(m, 0) = y
                        m = m + 1
                End If
                k = 0
        Next y
        Range("a1:b30000").Value = Arry1D_Primer
End Sub
-------------------- End VBA --------------------

千年沉寂

#!/usr/bin/env python y=3 myfile=open('primes.txt','w') while y<40000:     x=2     while x<y**0.5:         if y%x==0:break         x+=1     else:         print >>myfile,y,     y+=1 myfile.close() print 'End'

Python代码这样写就快多了，1秒就够了

forceps

确实快了许多。但更郁闷的事情在后头，我按你的python代码改写了vba，现在vba只要0.235秒。还是比python快很多。

Sub find_prime()     Dim Arry1D_prime(1 To 30000, 0 To 1)     Dim y As Long, x As Long, m As Long, n As Long     Dim i As Integer, k As Integer     k = 0     m = 1     For y = 3 To 40000         i = Int(Sqr(y))         For x = i To 2 Step -1             If (y Mod x) = 0 Then                 k = 1                 Exit For             End If         Next x         If k = 0 Then             Arry1D_prime(m, 0) = y             m = m + 1         End If         k = 0     Next y     Range("a1:b30000").Value = Arry1D_prime End Sub

[ 本帖最后由 forceps 于 2009-4-12 07:59 编辑 ]

千年沉寂

运行这段代码差不多0.35S左右。Python和VBA都是解释执行的，Python的优势并非执行效率而是开发效率。如果需要很高效率的话直接用C写扩展就行了。

import time y=3 start=time.time() myfile=open('primes.txt','w') while y<40000:     x=2     limit=y**0.5     while x<=limit:         if y%x==0:break         x+=1     else:         print >>myfile,y,     y+=2 myfile.close() end=time.time() print "End!\nUse Time:%f s" % (end-start)

forceps

很有道理。只是太不甘心了。
另外，你是什么机器，真快！我是1.8G的Intel Core 2 duo, 运行cygwin下的python，你上面的python在我的机器上运行时间是1.67s, 是vba(0.23s)的7.2倍。继续改：

#!/usr/bin/env python myfile=open('primes.txt','w') for y in range(3,500001):     i = int(y**0.5)+1     for x in range(2,i):         if y%x == 0 : break         x -= 1     else:         print >> myfile,y     y += 1 myfile.close()

这一次算到50万以下，耗时19.4s，是Excel vba(5.7s)的3.4倍。

efhilt

你应该在Linux下面，用native的Python和Wine下面的Excel的VBA比较。这样比较公平。

千年沉寂

机子配置差不多吧，也可能是你电脑运行太多程序的缘故，我是windows下Python直接运行的。

gawk

楼主你python有往文件里写入的，vba好像没有，不熟vba，你看看是不是这个原因

forceps

文件没存盘，但是写到单元格里了。见附件。

gawk

楼主，要不你写一个2个都不文件的也不写到单元格的比较一下，看看时间差距是多少