OFFSET
1,1
MAPLE
with(numtheory): P:=proc(q) local a, d, k, i, n, t;
for n from 1 to q do t:=0; a:=convert(n, base, 10);
for k from 1 to nops(a) do d:=0; for i from 1 to nops(a) do
if i<> k then d:=d*10+a[-i]; fi; od; t:=sigma(d)-d+t; od;
if n=t then print(n); fi; od; end: P(10^7)
CROSSREFS
KEYWORD
nonn,base,more
AUTHOR
Paolo P. Lava, Jul 03 2024
EXTENSIONS
a(6)-a(13) by Giovanni Resta, Jul 06 2024
STATUS
approved