OFFSET
1,3
FORMULA
a(n) = A003052(n) - n.
EXAMPLE
a(6) = 20 - 6 = 14.
a(7) = 31 - 7 = 24.
MAPLE
A007953 := proc(n) add(d, d=convert(n, base, 10)) ; end:
isA003052 := proc(n) for k from 1 to n do if k+A007953(k) = n then RETURN(false) ; fi; od: 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:
for n from 1 to 100 do printf("%d, ", A163128(n)) ; od: # R. J. Mathar, Jul 31 2009
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Juri-Stepan Gerasimov, Jul 21 2009
EXTENSIONS
Entries checked by R. J. Mathar, Jul 31 2009
STATUS
approved