OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(6) = 2 + 0 = 2;
a(7) = 3 + 1 = 4.
MAPLE
A007953 := proc(n) add( d, d= convert(n, base, 10) ); end:
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:
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:
A163125 := proc(n) A007953( A003052(n)) ; end: seq(A163125(n), n=1..100) ; # R. J. Mathar, Jul 27 2009
MATHEMATICA
sumdig[n_] := Plus @@ IntegerDigits[n]; m = 500; sumdig /@ Complement[Range[m], Table[n + sumdig[n], {n, 1, m}]] (* Amiram Eldar, Nov 28 2020 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Juri-Stepan Gerasimov, Jul 21 2009
EXTENSIONS
Values after a(51) corrected by R. J. Mathar, Jul 27 2009
STATUS
approved