Разложение на множители

def delit(a) :
    res = []
    i = 1
    while i * i < a + 1 :
        if a % i == 0 :
            temp = [i, a // i]
            res.append(temp)
        i += 1
    return res 
    
n = int(input())
npr = [2,3,5,7,11,13,17,19,23,29]
 
# первичное разложение на парные делители
fact = delit(n)
 
flen = len(fact)
while True :
    for j in range(flen) :
        zero = 1
 
        # если в очередном списке делителей 
        #есть "1", пропустить
        if 1 not in fact[j] :
            zero = 0
            n1 = list(fact[j])
 
            # "вытаскиваем" последний делитель из списка
            # и ищем все попарные разложения на множители
            factj = delit(n1.pop())
 
            # каждую полученную пару 
            #объединяем со списком
            for k in range(len(factj)) :
                factj[k].extend(n1)
                factj[k].sort()
                fact.append(factj[k])
    flen = len(fact)
 
    # если zero не изменилось во всех списках,
    # все варианты найдены
    if zero == 1:
        break
 
# удаляем 1 из всех вариантов
for ik in fact :
    ik.sort(reverse=True)
    if ik[-1] == 1 :
        ik.pop()
 
# удаляем одинаковые списки
j =[0]
resj = []
for i in fact :
    if i != j :
        resj.append(i)
    j=i
 
######################################

fact.sort()

npr = [2,3,5,7,11,13,17,19,23,29]

from functools import reduce
from operator import mul
npr = [2,3,5,7,11,13,17,19,23,29]
 
 
def prime_factors(n):
    primes_iter = iter(npr)
    remainder = n
    result = []
 
    while remainder > 1:
        factor = next(primes_iter)
        while remainder % factor == 0:
            result.append(factor)
            remainder //= factor
 
    # lets check:
    assert reduce(mul, result) == n
    return result
 
 
print(prime_factors(720))