OFFSET
1,3
COMMENTS
a(10) > 10^4 if it exists. - Michael S. Branicky, Sep 25 2025
a(10) > 12600. - Brenden W. Lamp, Feb 02 2026
a(10) > 25000. - Michael S. Branicky, Feb 05 2026
EXAMPLE
5 is a term because 5^5 + 6^4 = 4421, which is prime.
PROG
(Python)
from sympy import isprime
def ok(n): return n == 0 or isprime(n**n + (n+1)**(n-1))
print([k for k in range(100) if ok(k)])
(PARI) isok(k) = ispseudoprime(k^k + (k+1)^(k-1)); \\ Michel Marcus, Sep 25 2025
CROSSREFS
KEYWORD
nonn,more,hard
AUTHOR
Brenden W. Lamp, Sep 25 2025
STATUS
approved
