PRPs written as PARI/GP expressions

For n up to 250k, norm(Mod(x,x^4+x^3-x+1)^n-1)/2 is prime for n =

5, 7, 13, 17, 19, 23, 59, 83, 107, 229, 271, 367, 389, 491, 547, 691, 997, 1049, 1699, 1753, 1871, 2293, 3617, 5521, 7537, 9803, 10093, 10711, 14401, 14753, 56983, 163741, 181273, 207341

Search for larger PRPs of this form is ongoing.

Additional Lucas and Fermat PRP checks made with PFGW

To be completed...