pub fn factorize(n: u128) -> Vec<u128>
Return the sorted list of prime factors of n (with repetition). For n <= 1, returns an empty vector.
Uses trial division for small primes, then Pollard’s rho for large factors.