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...