I think this is the case gcd(a, p**n) = p**n. This comment and one below should be interchanged.