自己用python写了个分解质因数脚本,但是总是出错怎么办?

发布于 2022-09-01 05:33:21 字数 1926 浏览 10 评论 0

    #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 i
else:
    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多左右的质数,也就是说更新的不够多,但是我检查了很多遍代码都没搞明白问题出在哪里,。。

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

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。

评论(2

情域 2022-09-08 05:33:21

先回答你的问题,请检查你的第一个for循环的范围。range(1 * 10,2 * 10) 是从10到20的数字。但你想要的是从10到100的序列。

下面是这段代码更大的问题,而你并不知道这些问题的存在。
作为python新手,你问出这个问题,但实际上无论你用什么语言,这些错误都是存在的。问题不在于对语言的掌握程度,而是对“写程序”这一技能的掌握程度。
你的程序大概做了几件事:
1. 根据一个数字和一个质数列表,得到这个数字的质数因子
2. 输出质数因子
3. 获得从1到目标位数的质数列表
4. 根据已有数字位数,和所需数字位数,产生对应范围的质数列表
5. 将计算结果的位数和质数列表存储于某处
6. 从某处读取目前存储的质数列表和其对应的位数

在这段代码里所有的事情都混杂在一起。这样的坏处是,当一件事情变化时,其它部分也会受影响,即使它们并不真正与之相关。
想象一下有一个网站提供服务给出质数列表。你的代码会怎么改变
想象你的程序的输出是个页面
你能否写一段完全不用修改的代码可以用在上述不同场景?

当程序正确的分块后,你可以分别验证它。就算你自己找不到出错的地方,但你至少知道1,2,3,5,6都可以正常工作,因而查找范围或提问都会有针对性得多。

另有一个小问题,产生质数的代码中优化上限*0.5是不恰当的,应该为开平方。
验证自己所写的每行代码,而不是猜测它们如何工作

渔村楼浪 2022-09-08 05:33:21

这个求素数的算法貌似应该很经典啊。
求素数的算法应该是这样:
1)输入字符,假设为number, 判断number是否大于0
2)求平方根加1,设置为a
3) number对1~a之间的数字判断是否取余为0即可。
比如你输入15,那么取余+1,应该是4. 那么你就依次判断1~4之间的数是否能整除15即可。
python代码

import math


def 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 True


def getAllPrimes(number):
    retval = []
    for i in range(number):
        if isPrime(i):
            retval.append(i)
    return retval


if __name__ == '__main__':
    import time

    start = time.time()
    getAllPrimes(10 * 10000)
    print time.time() - start
  在我的机子上用了1.5秒求了10w以内的素数。
当然,求素数的算法有很多种,你这种将数据缓存起来的思路是对的,不过算法写的还是应该可以加强些,起码应该可读性更强些。
~没有更多了~
我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
原文