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

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()
 
# удаляем одинаковые списки
fact.sort()
print(len(fact))
j =[0]
resj = []
for i in fact :
    if i != j :
        resj.append(i)
    j=i
print(len(resj))
######################################
 
minc = int(1e300) + 1
for i in resj :
    j = 0
    ss = 1
    for ii in i:
        ss *= npr[j] ** (ii-1)
        if ss > int(1e300) + 1 :
            ss = int(1e300) + 1
            break
        j += 1
    minc = min(minc,ss)
print(minc)

def prime_divisors(n):
    a = []
    i = 2
    while i*i <= n:
        if n % i == 0:
            a.append(i)
            n //= i
        else:
            i += 1
    a.append(n)
    a.sort(reverse=True)
    return a
 
 
def solve(n):
    pr = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31]
    res = 1
    dv = prime_divisors(n)
    for x, p in zip(pr, dv):
        res *= x**(p-1)
    return res
 
print(solve(int(input())))

from itertools import combinations
#import numpy as np
 
def prime_divisors(n):
    a = []
    i = 2
    while i*i <= n:
        if n % i == 0:
            a.append(i)
            n //= i
        else:
            i += 1
    a.append(n)
    a.sort(reverse=True)
    return a
    
n = int(input())
npr = [2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61]
 
# поиск всех списков-вариантов множителей "n"
res =[]
# разложение на простые делители "n"
delit = prime_divisors(n)
tt =[n,1]
res.append(tt)
res.append(delit)
for i in range(2,len(delit)) :
    # все комбинации перебора делителей,
    # начиная с двух делителей до всех делителей
    data = list(combinations(delit, i))
    
    for j in data:
        jj = list(j)
        
        # здесь пытался применить prod от numpy,
        # но ниже при нахождении "ss" шло 
        # переполнение (????)
        
        # добавление к набору делителей
        # частного от "n" к произведению делителей
        temp = 1
        for k in jj:
            temp *= k
        temp1 = n//temp
        jj.append(temp1)
        jj.sort(reverse=True)
        
        # добавление набора "ttemp" из произведения 
        # делителей и частного от "n" к произведению
        # и удаление одинаковых списков-наборов
        if jj not in res:
            res.append(jj)
        ttemp = sorted([temp,temp1],reverse=True)
        if ttemp not in res:
            res.append(ttemp)
            
######################################
minc = int(1e300) + 1
for i in res :
    j = 0
    ss = 1
    for ii in i:
        ss *= npr[j] ** (ii-1)
        if ss > int(1e300) + 1 :
            ss = int(1e300) + 1
            break
        j += 1
    minc = min(minc,ss)
print(minc)

n = int(input('n = ')) 
dev = []  #Массив ВСЕХ делителей числа. Заполняется в другом куске кода
d = int(n / 2)
mas = list(range(d+1))  #Список чисел до n/2 так как делителей > n/2 нет
primes = []  # Массив простых делителей
 
 
def prosto():
    global dev, primes
    c = len(dev)
 
    for i in range(1, c, 1): # Поиск простых чисел начиная с первого простого делителя
        if dev[i] != 0:
            primes.append(dev[i])
            for help1 in range(i, c): #Само решето. Сравнивает значения массива ВСЕХ делителей с на кратность простым числам.
 
                for j in primes:
                    if dev[help1] % j == 0:  
                        dev[help1] = 0
 
    return primes