pub fn factorize(n: u128) -> Vec<u128>
Return the sorted list of prime factors of n (with repetition).