PRPs written as PARI/GP expressions

Fermat and Lucas PRP (PFGW)

norm(Mod(x,x^3-x-1)^n+1) is A001945(2*n)/A001945(n) in OEIS: http://oeis.org/A001945

To be completed...