OFFSET
1,1
MATHEMATICA
With[{c=(2^61-1)^2}, Select[Range[232000], AllTrue[c #+{1, -1}, PrimeQ]&]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Jun 17 2021 *)
PROG
(PARI) isok(n) = my(p); ispseudoprime(p=n*((2^61-1)^2) - 1) && ispseudoprime(p+2); \\ Michel Marcus, Sep 23 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Oct 07 2004
EXTENSIONS
More terms from Michel Marcus, Sep 23 2019
STATUS
approved