Пользователь вводит натуральное число N, найти все совершенные числа меньше N

#Проверка является ли число простым
def is_prime_number(p):
    i = 2
    while (i < p//2):
        if (p % i == 0):
            return False
        i += 1
    return True
 
#Ввод натурального числа 
n = -1
while n < 1:
    n = int(input())
    break
 
#Если совершенных чисел нет, программа выводит число 0. 
if n < 6:
    print (0)
else:
    #Нахождение всех совершенных чисел меньше либо равных n.
    p = 1
    k = 1
    while (p * 2 ** k <= n):
        p = 2 ** (k + 1) - 1
        if (is_prime_number(p)):
            print (p * 2 ** k, end = " ")
        k += 1

def f(n):
    res = 2
    for i in range(n-2):
        res *= 2
    res = res * ( res*2 - 1 )
    return res
#==============================================================================
def print_all_perfect_numbers_less_than(n):
    p = [2,3,5, 7,13,17, 19,31,61, 89,107,127, 521,607,1279, 2203,2281,3217, 4253,4423,9689, 9941,11213]
    for i in range( len(p) ):
        k = f( p[i] )
        if k < n:
            print(k)
        else:
            return
#==============================================================================
while True:
    print()
    n = int( input('n = ') )
    print_all_perfect_numbers_less_than(n)

def f(n):
    res = 1
    for i in range(n-1):
        res *= 2
    res = res * ( res*2 - 1 )
    return res
#==============================================================================
def print_all_perfect_numbers_less_than(n):
    p = [2,3,5, 7,13,17, 19,31,61, 89,107,127, 521,607,1279, 2203,2281,3217, 4253,4423,9689, 9941,11213]
    for i in range( len(p) ):
        k = f( p[i] )
        if k < n:
            print( i+1, ')\t', k, sep='' )
        else:
            return
#==============================================================================
while True:
    print()
    n = float( input('n = ') )
    print_all_perfect_numbers_less_than(n)

def sum_all_div(n):
    s=1
    i=2
    while(i*i<=n):
        if n%i==0:
            s=s+i
            k=n//i
            if k != i:
                s=s+k
        i=i+1
    return s
    
def print_perf(a,b):
    for x in range(a,b+1):
        if sum_all_div(x)==x:
            print(x)
            
print_perf(2,20000)

def perfect_numbers(number):
    answer = 2
    for i in range(number-2):
        answer *= 2
    answer = answer * ( answer*2 - 1 )
    return answer
 
def print_all_perfect_numbers_less_than(number):
    prime_numbers = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29,
                     31, 37, 41, 43, 47, 53, 59, 61, 67,
                     71, 73, 79, 83, 89, 97, 101, 103, 107,
                     109, 113, 127, 131, 137, 139, 149, 151,
                     157, 163, 167, 173, 179, 181, 191, 193,
                     197, 199, 211, 223, 227, 229, 233, 239]
    for i in range( len(prime_numbers) ):
        perfect = perfect_numbers( prime_numbers[i] )
        if perfect < number:
            print(perfect, end = ' ')           
        else:
            return
 
def no_perfect_numbers(number):
    if number < 6:
        print(0)    
 
number = 0
while number < 1:
    number = int(input())
    print_all_perfect_numbers_less_than(number + 1)
    no_perfect_numbers(number)
    break

def is_prime_2_pow(n, k):
    s = 4
    for i in range(k-2):
        s = (s*s - 2)%(n-1)
 
    return s==0
 
 
def perfects(n):
    a = 4
    res = [a//2*(a-1)]
    k = 2
    
    while a*(2*a-1) <= n:
        a *= 2
        k += 1
        if is_prime_2_pow(a, k):
            res.append(a//2*(a-1))
        
    print(*res)
 
 
n=int('1'+'0'*800)#int(input())
if n < 6:
    print(0)
else:
    perfects(n)

p = [2,3,5, 7,13,17, 19,31,61, 89,107,127, 521,607,1279, 2203,2281,3217, 4253,4423,9689, 9941,11213]
while True:
    print()
    n = float( input('n = ') )
    for q in p:
        k = 2**(q-1)*(2**q - 1)
        if k < n:
            print('-> ', k)
        else:
            break

a = [1, 2, 4, 6, 12, 16, 18, 30, 60, 88]
 
print([(1 << p) * ((2 << p) - 1) for p in a])