| Current Path : /usr/local/share/doc/pari/examples/ |
FreeBSD hs32.drive.ne.jp 9.1-RELEASE FreeBSD 9.1-RELEASE #1: Wed Jan 14 12:18:08 JST 2015 root@hs32.drive.ne.jp:/sys/amd64/compile/hs32 amd64 |
| Current File : //usr/local/share/doc/pari/examples/lucas.gp |
lucas(p) =
{
local(u,q); u=4; q=1<<p - 1;
for(k=3,p, u = (sqr(u)-2) % q);
u == 0
}