OFFSET
1,1
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..107
EXAMPLE
1583704193^(1+5+8+3+7+0+4+1+9+3)=1583704193 (mod 10^10) so 1583704193 is
in the sequence.
It is interesting that each of the four numbers 751^(7+5+1), 751^(7*5*1),
751^pi(751) and 751^prime(751) ends with 751.
MATHEMATICA
Do[n=Prime[m]; a=IntegerDigits[n]; If[PowerMod[n, Apply[Plus, a], 10^Length[a]]
==n, Print[n]], {m, 100000000}]
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Farideh Firoozbakht, Apr 28 2010
EXTENSIONS
Terms a(28) onward from Max Alekseyev, Aug 18 2013
STATUS
approved