%I #14 Sep 04 2022 17:00:37
%S 1,3,5,7,9,2,4,6,8,10,12,14,16,9,2,4,6,8,10,12,14,16,18,11,4,6,8,10,
%T 12,14,16,18,20,4,6,8,10,12,14,16,18,20,4,6,8,10,12,14,16,18,20,13,6,
%U 8,10,12,14,16,18,20,13,15,8,10,12,14,16,18,20,13,15,17,10,12,14,16,18,20,13
%N Sum of digits of the n-th Self-number (or Colombian number), A003052(n).
%H Amiram Eldar, <a href="/A163125/b163125.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A007953(A003052(n)).
%e a(6) = 2 + 0 = 2;
%e a(7) = 3 + 1 = 4.
%p A007953 := proc(n) add( d,d= convert(n,base,10) ); end:
%p isA003052 := proc(n) local k ; for k from 0 to n do if k+A007953(k) = n then RETURN(false): fi; od: RETURN(true) ; end:
%p A003052 := proc(n) option remember; if n = 1 then 1; else for a from procname(n-1)+1 do if isA003052(a) then RETURN(a) ; fi; od; fi; end:
%p A163125 := proc(n) A007953( A003052(n)) ; end: seq(A163125(n),n=1..100) ; # _R. J. Mathar_, Jul 27 2009
%t sumdig[n_] := Plus @@ IntegerDigits[n]; m = 500; sumdig /@ Complement[Range[m], Table[n + sumdig[n], {n, 1, m}]] (* _Amiram Eldar_, Nov 28 2020 *)
%Y Cf. A003052, A007953.
%K nonn,base
%O 1,2
%A _Juri-Stepan Gerasimov_, Jul 21 2009
%E Values after a(51) corrected by _R. J. Mathar_, Jul 27 2009
|