Напечатать все пары «близнецов» из отрезка

def isPrime(m):
    import numpy as np
    n = 2*m + 1
    pr = np.ones(n, dtype = 'int')
    for p in range(3, int(n**0.5) + 1, 2) :
        if pr[p] :
            pr[np.arange(p*p, n, 2*p)] = 0
    k = np.array(np.nonzero(pr))[0]
    k = k[(k%2 == 1) & (k>=m)]
    k = np.intersect1d(k, k + 2)
    k = np.union1d(k, k - 2)
    return k.reshape(k.size//2, 2)
 
n = int(input())
for nums in isPrime(n):
    print(*nums)

def search(n):
    def divisor(num):
        divisors = []
        for d in range(1, num + 1):  # divisors
            if num % d == 0:
                divisors.append(d)
        return divisors
    
    twins = []
    for i in range(n, n*2 + 1):         # segment
        if len(divisor(i)) == 2:
            if i + 2 <= n*2:
                if len(divisor(i + 2)) == 2:
                    twins.append(f'{i} и {i + 2}')
 
    if len(twins) != 0:
        print(', '.join(twins))
    else:
        print(f'Близнецы в диапазоне {n} -- {n * 2} отсутствуют')
 
 
while True:
    # noinspection PyBroadException
    try:
        number = int(input('Число для отрезка: '))
        if number <= 2:
            print('Введите число, большее 2')
            continue
    except Exception:
        print('Введите число, большее 2')
    else:
        search(number)
        print()