OFFSET
1,2
COMMENTS
a(11) > 10^9. - Giovanni Resta, Mar 29 2017
EXAMPLE
sigma(1) = 1. Its digit reverse is again 1 and sigma(1) = 1 = 1 * 1;
sigma(16) = 31 and sigma(61) = 62 = 2 * 31;
sigma(15) = 24 and sigma(51) = 72 = 3 * 24;
sigma(17) = 18 and sigma(71) = 72 = 4 * 18; etc.
MAPLE
with(numtheory): R:=proc(w) local x, y, z; x:=w; y:=0; for z from 1 to ilog10(x)+1 do y:=10*y+(x mod 10); x:=trunc(x/10); od; y; end:
P:=proc(q) local k, n; for k from 1 to q do for n from 1 to q do
if k*sigma(n)=sigma(R(n)) then print(n); break; fi; od; od; end: P(10^9);
MATHEMATICA
rev[n_] := FromDigits@ Reverse@ IntegerDigits@ n; s[n_] := DivisorSigma[1, n]; a[n_] := Block[{k, v}, For[k = 1, Mod[v = s@rev@k, n] > 0 || v != n s@k, k++]; k]; Array[a, 10] (* Giovanni Resta, Mar 29 2017 *)
CROSSREFS
KEYWORD
nonn,more,base
AUTHOR
Paolo P. Lava, Mar 28 2017
STATUS
approved