A276333 The most significant digit in greedy A001563-base (A276326): a(n) = floor(n/A258199(n)), a(0) = 0.

0, 1, 2, 3, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 1

Offset

0,3

Links

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

Formula

a(0) = 0; for n >= 1, a(n) = floor(n/A258199(n)).

Mathematica

Table[Floor[n/(# #!)] &@ NestWhile[# + 1 &, 0, # #! <= n &[# + 1] &], {n, 96}] (* or *)
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[First@ f@ n, {n, 96}] (* Michael De Vlieger, Aug 31 2016 *)

Prog

(Scheme) (define (A276333 n) (if (zero? n) n (floor->exact (/ n (A258199 n)))))

Crossrefs

Cf. A001563, A276326, A258199, A276334, A276335.

Sequence in context: A327462 A255824 A122587 * A347380 A308638 A356320

Adjacent sequences: A276330 A276331 A276332 * A276334 A276335 A276336

Keyword

nonn

Author

Antti Karttunen, Aug 30 2016

Status

approved