A163128 - OEIS

%I #9 Feb 10 2019 23:03:00

%S 0,1,2,3,4,14,24,34,44,54,64,74,84,94,95,105,115,125,135,145,155,165,

%T 175,185,186,196,206,216,226,236,246,256,266,276,277,287,297,307,317,

%U 327,337,347,357,367,368,378,388,398,408,418,428,438,448,458,459,469,479,489,499

%N a(n) is the n-th self-number minus n.

%F a(n) = A003052(n) - n.

%e a(6) = 20 - 6 = 14.

%e a(7) = 31 - 7 = 24.

%p A007953 := proc(n) add(d,d=convert(n,base,10)) ; end:

%p isA003052 := proc(n) for k from 1 to n do if k+A007953(k) = n then RETURN(false) ; fi; od: 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 A163128 := proc(n) A003052(n)-n ; end:

%p for n from 1 to 100 do printf("%d,",A163128(n)) ; od: # _R. J. Mathar_, Jul 31 2009

%Y Cf. A000027, A003052.

%K nonn,base

%O 1,3

%A _Juri-Stepan Gerasimov_, Jul 21 2009

%E Entries checked by _R. J. Mathar_, Jul 31 2009