OFFSET
1,2
COMMENTS
No more terms < 6700.
No more terms < 10^4. - Daniel Suteu, Mar 20 2020
PROG
(Python)
from sympy.ntheory.generate import primorial, isprime
A333390_list = [n for n in range(1, 100) if isprime(primorial(2*n, nth=False)*2**n-1)] # Chai Wah Wu, Mar 18 2020
(PARI) is(k) = ispseudoprime(lcm(primes([2, 2*k]))*2^k-1); \\ Jinyuan Wang, Mar 20 2020
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
L. Joris Perrenet, Mar 18 2020
EXTENSIONS
a(31) from Jinyuan Wang, Mar 20 2020
a(32) from Daniel Suteu, Mar 20 2020
a(33)-a(36) from Michael S. Branicky, Jun 26 2024
STATUS
approved