A276331 a(0) = 0 and for n >= 1, a(n) = n minus the largest term of A001563 that divides n.

0, 0, 1, 2, 0, 4, 5, 6, 4, 8, 9, 10, 8, 12, 13, 14, 12, 16, 0, 18, 16, 20, 21, 22, 20, 24, 25, 26, 24, 28, 29, 30, 28, 32, 33, 34, 18, 36, 37, 38, 36, 40, 41, 42, 40, 44, 45, 46, 44, 48, 49, 50, 48, 52, 36, 54, 52, 56, 57, 58, 56, 60, 61, 62, 60, 64, 65, 66, 64, 68, 69, 70, 54, 72, 73, 74, 72, 76, 77, 78, 76, 80, 81, 82, 80

Offset

0,4

Links

Antti Karttunen, Table of n, a(n) for n = 0..4320

Formula

a(n) = n - A276330(n).

Mathematica

f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], (# #!) &[# - i]]], {i, 0, # - 1}] &@ NestWhile[# + 1 &, 0, (# #!) &[# + 1] <= n &]; Rest[a][[All, 1]]]; {0}~Join~Table[n - # #! &[Length@ TakeWhile[Reverse@ f@ n, # == 0 &] + 1], {n, 84}] (* Michael De Vlieger, Aug 31 2016 *)

Prog

(Scheme) (define (A276331 n) (- n (A276330 n)))

Crossrefs

Cf. A001563, A276330, A276332.

Sequence in context: A111677 A326186 A326055 * A049271 A004178 A068333

Adjacent sequences: A276328 A276329 A276330 * A276332 A276333 A276334

Keyword

nonn

Author

Antti Karttunen, Aug 30 2016

Status

approved