Function power_mod 

pub const fn power_mod(p: u32, b: u32, e: u32) -> u32

Compute b^e mod p. This is a const version of Prime::pow_mod.