OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..1001
EXAMPLE
a(4) = 17 is a term because 17+1+7 = 25 = 5^2.
MAPLE
filter:= proc(n) local s, F;
s:= n + convert(convert(n, base, 10), `+`);
F:= ifactors(s)[2];
nops(F)=1 and F[1][2]>1
end proc:
select(filter, [$1..20000]);
PROG
(PARI) isok(k) = isprimepower(k + sumdigits(k)) > 1; \\ Michel Marcus, Mar 22 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
J. M. Bergot and Robert Israel, Mar 21 2021
STATUS
approved