Add large Perrin pseudoprime
Fix negative is_power for defined a, still need a=0 case