Python分解质因数脚本源码纠错求助，python脚本源码, #Num(ori

Python分解质因数脚本源码纠错求助，python脚本源码, #Num(ori

    #Num(original number); Lg(length of original number); Lf(length file);#Pl(primes list); Pf(prime file); Ps(primes); d(list of answers);#length(open txt file 'length'); Pb(open txt file 'primebase'); Qpt()Num=raw_input('Please input a number!\nNumber: ')Lg=len(Num)length=open('length.txt','r+')Lf=length.readlines()Pb=open('primebase.txt','r+')Pl=Pb.readlines()if int(Lg)>int(Lf[0]):    print 'No enough primes in database!\nUpdating the database...'    for i in range(int(Lf[0])*10,int(Lg)*10):        Qpt=0        for x in range(2,int(i*0.5)+1):            if i%x==0:                Qpt+=1                if Qpt!=0:                    break        if Qpt==0 and i!=0 and i!=1:            Pl.append('%s\n'%str(i))    Pf=open('primebase.txt','w+')    Pf.writelines(Pl)    Pf.close()    Lf[0]=str(Lg)    file=open('length.txt','w+')    file.writelines(Lf)    file.close()    print 'Finish updating database!'    d=[]    Num=int(Num)    while Num!=1:        for i in Pl:            i=int(i)            if Num%i==0:                Num=Num/i                d.append(i)    if d==[]:        print 'Prime number!'    else:        print 'Prime factors: '        for i in d:            print ielse:    d=[]    Num=int(Num)    while Num!=1:        for i in Pl:            i=int(i)            if Num%i==0:                Num=Num/i                d.append(i)    if d==[]:        print 'Prime number!'    else:        print 'Prime factors: '        for i in d:            print i

本人py新手。我的具体思路是，搞两个txt文件分别储存质数库和目前质数库的范围。利用len（）函数来判断用户输入的数字是几位数，如果位数超过了数据库的范围，就自行更新数据库。举例：数据库里面有2-10的所有质数，然后我现在输入20，超出范围，程序自行更新数据库到2-100的所有质数，再运算。这样就避免了在没有数据库的情况下因为数字太大运算效率低（牺牲第一次的运算效率来计算质数库）。但是实际使用的时候，我发现我输入100，然后程序开始自主更新数据库，但是更新出来的并不是2-100的所有质数，仅仅到了29。同理，在输入很大的数字的时候，也只能更新到300多左右的质数，也就是说更新的不够多，但是我检查了很多遍代码都没搞明白问题出在哪里，。。

*代码中的注释内容基本解释了所有的变量，方便大家看代码。。。

import mathdef isPrime(number):    assert number >= 0    if number == 0 or number == 1:        return False    sqrtNumber = int(math.sqrt(number)) + 1    for i in range(sqrtNumber, 1, -1):        if not number % i:            return False    return Truedef getAllPrimes(number):    retval = []    for i in range(number):        if isPrime(i):            retval.append(i)    return retvalif __name__ == '__main__':    import time    start = time.time()    getAllPrimes(10 * 10000)    print time.time() - start  在我的机子上用了1.5秒求了10w以内的素数。当然，求素数的算法有很多种，你这种将数据缓存起来的思路是对的，不过算法写的还是应该可以加强些，起码应该可读性更强些。

编橙之家文章，