A276330 a(n) = largest term of A001563 that divides n, a(0) = 0.

0, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 18, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 18, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 18, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 18, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, 1, 18, 1, 4, 1, 1, 1, 96

Offset

0,5

Links

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

Formula

a(0) = 0; for n >= 1, a(n) = A001563(A276329(n)).

Other identities. For all n >= 0:

A276331(n) = n - a(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[# #! &[Length@ TakeWhile[Reverse@ f@ n, # == 0 &] + 1], {n, 120}] (* Michael De Vlieger, Aug 31 2016 *)

Prog

(Scheme) (define (A276330 n) (if (zero? n) n (A001563 (A276329 n))))

Crossrefs

Cf. A001563, A258199, A276328, A276329, A276331.

Sequence in context: A366524 A366523 A365635 * A234957 A327939 A365837

Adjacent sequences: A276327 A276328 A276329 * A276331 A276332 A276333

Keyword

nonn

Author

Antti Karttunen, Aug 30 2016

Status

approved