OFFSET
1,1
COMMENTS
a(17) > 52000.
a(18) > 10^5. - Michael S. Branicky, Aug 20 2025
EXAMPLE
a(5) = 12 because 1 + 12^2 * 2^(12) + 12^3 * 3^(12) = 918919873 is prime.
MATHEMATICA
p = 31000
ParallelTable[
If[PrimeQ[1 + n^2*2^n + n^3*3^n], n, Nothing], {n, 0, p}]
PROG
(PARI) is(n)=ispseudoprime(1 + n^2*2^n + n^3*3^n)
(Magma) [k: k in [1..31000] | IsPrime(1 + k^2*2^k + k^3*3^k)];
(Python)
from sympy import isprime
print([k for k in range(52000) if isprime(1 + (2**k)*(k**2) + (3**k)*(k**3))])
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Enrico Masina, Sep 11 2022
EXTENSIONS
a(17) from Michael S. Branicky, Aug 20 2025
STATUS
approved
