Looks like I got it wrong. Requirements for `_sqrt_mod_prime_power()` are either `k == 1` or `a % p != 0 ` if `k > 1`.
So I still think `a % p == 0` should be changed to `a % p != 0`