Function factor_pk 

pub fn factor_pk<P: Prime>(p: P, n: u32) -> (u32, u32)

Factor $n$ as $p^k m$. Returns $(k, m)$.