OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
A prime p = abcdef... (each letter being a single digit) belongs to this sequence if sum = a^5+b^5+... is a prime. - Corrected by Robert Israel, Feb 22 2026
EXAMPLE
227 is a term since 2^5+2^5+7^5 = 32+32+16807 = 16871 is a prime.
MAPLE
filter:= proc(n) local t;
isprime(n) and isprime(add(t^5, t=convert(n, base, 10)))
end proc:
select(filter, [seq(i, i=3..10000, 2)]); # Robert Israel, Feb 22 2026
MATHEMATICA
Select[Prime[Range[500]], PrimeQ[Total[IntegerDigits[#]^5]]&] (* Harvey P. Dale, May 25 2011 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Carmine Suriano, Sep 02 2010
STATUS
approved
