Fix probable_prime over large shift