OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
If 2m+1 is prime and is in the sequence, 2m+1 = gcd(m!+1, 2^m+1).
MATHEMATICA
Table[GCD[n!+1, 2^n+1], {n, 0, 4000}]/.(1->Nothing) (* Harvey P. Dale, Aug 03 2018 *)
PROG
(PARI) for(k=0, 3000, d=gcd(k!+1, 2^k+1); if(d<>1, print1(d, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Feb 03 2002
EXTENSIONS
Corrected and extended by Rick L. Shepherd, May 20 2002
Offset corrected by Amiram Eldar, Jun 06 2022
STATUS
approved