![]() ![]() | ![]() |
File: [local] / OpenXM_contrib / pari / examples / Attic / lucas.gp (download)
Revision 1.1.1.1 (vendor branch), Sun Jan 9 17:35:30 2000 UTC (24 years, 7 months ago) by maekawa
Import PARI/GP 2.0.17 beta. |
lucas(p) = { local(u,q); u=4; q=1<<p - 1; for(k=3,p, u = (sqr(u)-2)%q); u == 0 }