pub fn primitive_roots(n: u64) -> Vec<u64>
Find all primitive roots modulo n (if they exist). Primitive roots exist iff n is 1, 2, 4, p^k, or 2*p^k for odd prime p.