PRPs written as PARI/GP expressions

Fermat and Lucas PRP (PFGW)

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

Also see http://oeis.org/A104499 (indices of prime terms of A001945)

To be completed...