def smoothness(n):

    if n == 1:
        return (1, 1)  # not prime, but otherwise this causes headaches
    facs = factorint(n)
    return max(facs), max(m**facs[m] for m in facs)
Return the B-smooth and B-power smooth values of n.
dataset/ETHPy150Open sympy/sympy/sympy/ntheory/factor_.py/smoothness