pub fn prime_factorize(n: usize) -> BTreeMap<usize, usize>