Мультивселенная. Выбрать значения, у которых ровно два делителя, не считая единицу и само число

E = int(input())
F = int(input())
A = []
for B in range(E, F + 1):
    C = 0
    for D in range(B - 1, 1, -1):
        if B % D == 0: C += 1; A.append(D)
    if C == 2: print(min(A[0], A[1]), max(A[0], A[1]))
    A = []

E = int(input())
F = int(input())
A = []
for B in range(E, F + 1):
    C = 0
    for D in range(B // 2, 1, -1):
        if B % D == 0: C += 1; A.append(D)
    if C == 2: print(min(A[0], A[1]), max(A[0], A[1]))
    A = []

def prns(n):
    list_true = [True] * n
    for o in range(3, (int(n ** 0.5) + 1), 2):
        if list_true[o]: list_true[o * o::2 * o] = [False] * ((n - o * o - 1) // (2 * o) + 1)
    return [2] + [u for u in range(3, n, 2) if list_true[u]]
 
 
def check(n):
    if n & 1 == 0: return False
    d = 3
    while d * d <= n:
        if n % d == 0: return False
        d = d + 2
    return True
 
 
def solution():
    e = int(input());
    f = int(input());
    k = 1
    if e > f:
        k = -1
    prn = prns(int(max(e, f) ** 0.5))
    for n in range(e, f + k, k):
        for i in prn:
            if (i*i < n and not n % i and
             check(n // i)):
                 print(i, n // i)
                 break
 
solution()

def prns(n):
    list_true = [True] * n
    for o in range(3, (int(n ** 0.5) + 1), 2):
        if list_true[o]: list_true[o * o::2 * o] = [False] * ((n - o * o - 1) // (2 * o) + 1)
    return [2] + [u for u in range(3, n, 2) if list_true[u]]
 
 
def check(n):
    if n == 2: return True
    if n & 1 == 0: return False
    d = 3
    while d * d <= n:
        if n % d == 0: return False
        d = d + 2
    return True
 
 
def solution():
    e = int(input());
    f = int(input());
    k = 1
    if e > f:
        k = -1
    prn = prns(int(max(e, f) ** 0.5))
    for n in range(e, f + k, k):
        for i in prn:
            if (i*i <= n and not n % i and
             check(n // i)):
                 print(i, n // i)
                 break
 
solution()

def prns(n):
    list_true = [True] * n
    for o in range(3, (int(n ** 0.5) + 1), 2):
        if list_true[o]: list_true[o * o::2 * o] = [False] * ((n - o * o - 1) // (2 * o) + 1)
    return [2] + [u for u in range(3, n, 2) if list_true[u]]
 
 
def check(n):
    if n & 1 == 0: return False
    d = 3
    while d * d <= n:
        if n % d == 0: return False
        d = d + 2
    return True
 
 
def solution():
    e = int(input());
    f = int(input());
    k = 1
    if e > f:
        k = -1
    prn = prns(int(max(e, f) ** 0.5))
    for n in range(e, f + k, k):
        for i in prn:
            if (i*i < n and not n % i and
             (check(n // i) or n == i*i*i)):
                 print(i, n // i)
                 break
 
solution()

i*i < n