According to line 240 and 242, it seems like `_sqrt_mod_prime_power()` is used only when `a % p != 0`.
So I think a % p == 0 here should be `a % p != 0`.