OFFSET
1,2
COMMENTS
1# = 2 2# = 2*3 = 6 3# = 2*3*5 = 30.
No more terms < 5000. - L. Joris Perrenet, Mar 17 2020
EXAMPLE
a(1) = 1 because 1#*2^1 - 1 = 3 is prime.
a(2) = 2 because 2#*2^2 - 1 = 23 is prime.
MATHEMATICA
For[n = 1, n < 60, n++, If[PrimeQ[2^n*Product[Prime[i], {i, 1, n}] - 1], Print[n]]] (* Stefan Steinerberger, Feb 06 2006 *)
PROG
(PARI) pp(n)= s=1; for(i=1, n, s=s*prime(i)); return(s);
f(n)=pp(n)!*2^n -1;
for (i=1, 500, if(isprime(f(i)), print1(i, ", ")))
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 02 2004
EXTENSIONS
a(17) from Stefan Steinerberger, Feb 06 2006
a(18)-a(22) from L. Joris Perrenet, Mar 17 2020
a(23)-a(24) from Michael S. Branicky, Aug 28 2024
STATUS
approved