OFFSET
1,1
MAPLE
with(numtheory): P:=proc(q) local a, b, d, i, j, k, n, t;
for n from 1 to q do t:=0; a:=convert(n, base, 10);
for b from 1 to nops(a) do d:=0;
for i from 1 to nops(a) do if i<>b then d:=d*10+a[-i]; fi; od; k:=0; j:=d;
if d>0 then while j mod 2<>1 do k:=k+1; j:=j/2; od;
t:=sigma(2*d+1)+sigma(2*d-1)+sigma(d/2^k)*2^(k+1)-6*d-2+t; fi; od;
if n=t then print(n); fi; od; end: P(10^8);
CROSSREFS
KEYWORD
nonn,base,more
AUTHOR
Paolo P. Lava, Jul 28 2024
STATUS
approved