pub fn generate_prime(bits: u32) -> BigUint
Generate a random probable prime (Miller-Rabin, 20 rounds).