OFFSET
1,1
COMMENTS
This sequence was arrived at by doing the arithmetic on primes in their increasing order.
The sequence treats zero terms as generating a "1" digit; i.e., 0^0 = 1. - Harvey P. Dale, Dec 16 2025
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
Dario Alpern, Factorization using the Elliptic Curve Method
Chris Caldwell, The First 10,000 Primes
EXAMPLE
31 from 23: 2^2 + 3^3 = 31 (prime);
823799 from 47: 4^4 + 7^7 = 823799 (prime);
826669 from 157: 1^1 + 5^5 + 7^7 = 826669 (prime).
MATHEMATICA
Select[Table[With[{c=IntegerDigits[p]/.(0->1)}, Total[c^c]], {p, Prime[Range[500]]}], PrimeQ] (* Harvey P. Dale, Dec 16 2025 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Parthasarathy Nambi, Aug 20 2009
EXTENSIONS
More terms from Zak Seidov, Aug 21 2009
STATUS
approved
