OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 11 because 1^4 + 1^4 = 2 which is prime.
a(10) = 89 because 8^4 + 9^4 = 10657 which is prime.
MATHEMATICA
upto=500; Select[Prime[Range[upto]], PrimeQ[Total[IntegerDigits[#]^4]]&] (* Paolo Xausa, Nov 23 2021 *)
PROG
(Python)
from sympy import isprime, primerange
def ok(p): return isprime(sum(int(d)**4 for d in str(p)))
def aupto(limit): return [p for p in primerange(1, limit+1) if ok(p)]
print(aupto(1433)) # Michael S. Branicky, Nov 23 2021
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Chuck Seggelin (barkeep(AT)plastereddragon.com), Jan 03 2004
STATUS
approved