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

forceps

原帖由 gawk 于 2009-4-12 11:06 发表
楼主，要不你写一个2个都不文件的也不写到单元格的比较一下，看看时间差距是多少

改了，不写文件也不写单元格，还是Cygwin python vs Excel vba

Cygwin python: 19.3s
Excel vba: 5.5s
19.3s / 5.5s = 3.5

gawk

的确是慢，我也试了一下

yaksavage

有意思啊，太有意思了
怎么都在讨论python的速度啊
看来python真的火了

附：
一个简单的排序程序 c++的居然没有python 很郁闷 这是为什么呢
https://groups.google.com/group/pongba/browse_thread/thread/840dd3bba4ce2762/9f8ba8fc9ff12f63#

ps:搞不懂，为什么要用Cygwin python，不是有原装的msc编译的吗

[ 本帖最后由 yaksavage 于 2009-4-12 12:32 编辑 ]

gawk

楼主，你的python代码中为什么要x-=1 y+=1?

gawk

楼主你把range换成xrange试试
我测了一下快好多了

caesarok

楼主,你还需要调用个模块
import psyco
psyco.full()

加上这段绝对快!

gawk

1. 
   
2. import time
   
3. y=3
   
4. start=time.time()
   
5. myfile=open('primes.txt','w')
   
6. for y in xrange(3,500001):
   
7.     i = int(y**0.5)+1
   
8.     for x in xrange(2,i):
   
9.         if y%x == 0 : break
   
10.     else:
    
11.         print >> myfile,y
    
12. myfile.close()
    
13. end=time.time()
    
14. print "End!\nUse Time:%f s" % (end-start)
    

复制代码

forceps

原帖由 gawk 于 2009-4-12 12:51 发表
楼主，你的python代码中为什么要x-=1 y+=1?

对啊。
先前是用while x>1, 这样x-=1递减x，改成for x in range[]就不用了。我是新手，多谢！

另外回13楼，其实我折腾的对象主要是debian，对我这样蜻蜓点水，今天弄弄gcc，明天试试python，偶尔还要试试新东西，比如有兄弟提到的psyco模块，用debian可能更方便一些。办公机器，只能模拟一下linux，我用coLinux装的Debian Lenny，但coLinux还是要稍许慢一些。同样求质数要多花至少1.5倍的时间。听说cygwin速度已经很接近windows下直接执行的速度了，而相比用windows下的python，用cygwin又和Debian更接近一点。我懒得重新适应 windows下的python，所以和vba比较速度的时候就用相对coLinux更快一点的cygwin python计时。我也是瞎折腾，和诸位大虾靠python吃饭不一样。

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

forceps

原帖由 gawk 于 2009-4-12 13:34 发表

import time
y=3
start=time.time()
myfile=open('primes.txt','w')
for y in xrange(3,500001):
    i = int(y**0.5)+1
    for x in xrange(2,i):
        if y%x == 0 : break
    else:
        ...

果真又快了，现在求50万以下的质数耗时9.58s，换用xrange()快了一倍。但还是比vba慢。9.58s/5.5s = 1.74。

下一步试试psyco。

千年沉寂

用筛选法求质数，1000000以内的用了0.765s.

import time start_time=time.time() myfile=open('primes.txt','w') NUM =1000000 prime_num = [i for i in xrange(NUM+1)] for i in xrange(2,int(NUM**0.5)+1):   if prime_num[i]:       start = i**2       step  = i       prime_num[start::step] = ((NUM - start)/step + 1)*[0] print >> myfile,[i for i in prime_num if prime_num[i]!=0] myfile.close() end_time=time.time() print "End!\nUse Time:%f s" % (end_time-start_time)