OFFSET
1,1
COMMENTS
Subsequence of A000961.
The corresponding values of the primes p are 2, 3, 5, 7, 7, 17, 31, 43, 53, 71, 97, 107, 103, 107, 127, 107, ... .
The corresponding values of the powers m are 1, 1, 1, 1, 4, 3, 7, 7, 7, 9, 10, 11, 13, 13, 14, 15, ... .
EXAMPLE
2401 is in the sequence because 2401 = 7^4 and 2+4+0+1 = 7.
MAPLE
nn:=2000:
for p from 1 to 15 do :
for m from 1 to nn do:n:=ithprime(m)^p:
y:=convert(n, base, 10):n1:=nops(y):
s1:=sum('y[i]', 'i'=1..n1):
if s1=ithprime(m)
then
printf(`%d, `, n):
else
fi:
od:
od:
MATHEMATICA
a252357[n_Integer] :=
Select[Range[n],
PrimeQ[Plus @@ IntegerDigits[#]] &&
IntegerQ[Log[Plus @@ IntegerDigits[#], #]] &]; a252357[10^5] (* Michael De Vlieger, Dec 17 2014 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Dec 17 2014
STATUS
approved