在使用python解决问题之前,我们先说一下,什么是素因子分解
所谓素因子分解就是,先找这个数的所有约数(约数即:a%b == 0,也就是a可以被b整除)
例如:20的约数集合为 [1, 2, 5, 10, 20]
那么素因子分解呢?
就是从最小的素数约数开始除,也就是这个除数要满足两个条件,一是约数,二是素数
那么这里20的最小的素约数是2,所以我们从2开始除,并且一直除到不能被整出为止:
num = 20
num = num / 2
num = 10(这里num依旧可以被2整除,所以再来一次)
num = num / 2
num = 5 (num很明显除以2不能整除)
所以接下来被除数需要向后走,即5(再取下一个数为除数之前,要先判断是否为素数,这里5为素数,如果不是素数比如4则需要继续向后取)
num = num / 5
num = 1
至此,运算结束。所以,我们得到20的素因子集合为:[2, 2, 5]
那么接下来就是使用python实现的部分了
通过之前的解析过程,我们可以看到我们需要这么两个小模块:
1. 判断数字是否是素数,如果数字本身是素数,那么直接返回1和它本身就可以了,因为素数不能被分解了,即 [1, num]
2. 我们需要获得这个数的所有约数,以便于我们进行作除法
下面是这两个部分的实现:
# 判断是否是素数
def isprime(num):
count = num / 2
while count > 1:
if num % count == 0:
return False
count -= 1
else:
return True
# 得到所有的约数
def getfactors(num):
return [x for x in range(1, num) if num % x == 0]
通过上面两个部分,我们就可以得到所有的基础了,接下来就是运算了
def primefactor(num):
if isprime(num):
return [1,num]
factors = getfactors(num)
retList = []
consult = num
for i in range(1,len(factors)):
if consult == 1:
break
while True:
if consult % factors[i] != 0:
break
if isprime(factors[i]):
consult /= factors[i]
retList.append(factors[i])
else:
break
return retList
ok,最后调用primefactor就可以得到结果了,注意返回的是个list对象哦