Obtaining the Prime Factors of a Semiprime

semiprime_factors.py

import math


def obtain_pq(n):
    p = 0.1
    k = math.ceil(math.sqrt(n))
    while not p.is_integer():
        p = k - math.sqrt((k ** 2) - n)
        k += 1

    q = int(n / p)
    p = int(p)

    return p, q