pub fn is_primitive_root(g: u64, p: u64) -> bool
Check whether g is a primitive root mod p (p must be prime).
g is a primitive root mod p iff g^((p-1)/q) ≢ 1 (mod p) for every prime q | (p-1).