/* -*- gp-script -*- */ \\% Pocklington-Lehmer test. \\ Author: Joerg Arndt \\ License: GPL version 3 or later \\ online at http://www.jjj.de/pari/ \\ version: 2014-October-16 (18:30) pocklington_lehmer(F, u, c=10000)= { /* Pocklington-Lehmer test for the primality of p=f*u+1. * Return last successful base, else zero. * F must be the factorization of f. * Test bases a=2...c * Must have u