Нахождение простых делителей числа

def isPrime(current): #проверяю делители числа на простоту
    for i in range(2, int(current**0.5)+1):
        if current % i == 0:
            return False
    return True
 
def dividers(current): #нахожу делители числа, и если они простые, то добавляю в список
    arr = []
    for i in range(2, int(current**0.5)+1):
        if current % i == 0:
            if isPrime(i):
                arr.append(i)
                if i < current // i and isPrime(current//i):
                    arr.append(current//i)
    return arr
 
count = 0
for i in range(10001, 50001):
    divid = dividers(i)
    print(i, divid)
    if len(divid) == 3:
        count+=1
print(count)

def primes(n):
    sieve = [True] * n 
    for i in range(3, int(n ** 0.5) + 1, 2):
        if sieve[i]:
            sieve[i * i::2 * i] = [False] * ((n - i * i - 1) // (2 * i) + 1)
    return [2] + [i for i in range(3, n, 2) if sieve[i]] 
 
pr = primes(50000//6+1)
cnt = 0
for num in range(10001, 50001):
    ddd= num
    d = 0
    f1 = f2 = f3 = False
    for p in pr:
        if num%p == 0:
            if f2:
                f3 = True
            elif f1:
                f2 = True
            while num%p==0:
                f1 = True
                num //= p 
        if f3:
            if num == 1:
                cnt += 1
            break  
print(cnt)

import numpy as np
def primes(n):
    sieve = [True] * n 
    for i in range(3, int(n ** 0.5) + 1, 2):
        if sieve[i]:
            sieve[i * i::2 * i] = [False] * ((n - i * i - 1) // (2 * i) + 1)
    return [2] + [i for i in range(3, n, 2) if sieve[i]] 
 
pr = primes(50000//6+1)
res = np.zeros(50001)
for p in pr:
    res[p::p] += 1
print(np.nonzero(res[10001:] == 3)[0].shape[0])