A094634 a(n) = A094371(n) / A094372(n).

1, 3, 4, 5, 4, 5, 5, 2, 6, 7, 5, 6, 6, 7, 7, 6, 7, 7, 7, 6, 7, 7, 7, 8, 8, 7, 8, 8, 8, 7, 8, 8, 8, 9, 9, 8, 9, 9, 8, 9, 9, 9, 10, 9, 8, 9, 9, 9, 10, 9, 10, 10, 9, 2, 10, 10, 9, 10, 10, 10, 11, 10, 11, 9, 10, 10, 10, 11, 10, 11, 11, 10, 11, 11, 10, 11, 11, 11, 11, 10, 11, 11, 11, 11, 12, 11, 10

Offset

1,2

Comments

When a(n)=2, A094404(n) is not 1.

Links

Table of n, a(n) for n=1..87.

Mathematica

Kempner[1] := 1; Kempner[n_] := Max[Kempner @@@ FactorInteger[n]]; Kempner[p_, 1] := p; Kempner[p_, alpha_] := Kempner[p, alpha] = Module[{a, k, r, i, nu, k0 = alpha(p - 1)}, i = nu = Floor[Log[p, 1 + k0]]; a[1] = 1; a[n_] := (p^n - 1)/(p - 1); k[nu] = Quotient[alpha, a[nu]]; r[nu] = alpha - k[nu]a[nu]; While[r[i] > 0, k[i - 1] = Quotient[r[i], a[i - 1]]; r[i - 1] = r[i] - k[i - 1]a[i - 1]; i-- ]; k0 + Plus @@ k /@ Range[i, nu]]; L = M = {}; a = 2; Do[ s = Kempner[n]; If[s/n < a, a = s/n; AppendTo[M, a]; AppendTo[L, n]], {n, 1, 10^5}]; L/Denominator[M] (* revised by Eric W. Weisstein, May 17 2004, and by N. J. A. Sloane, Jul 03 2017 *)

Crossrefs

Cf. A002034, A094371, A094372.

Sequence in context: A155078 A332420 A115051 * A114545 A309704 A232702

Adjacent sequences: A094631 A094632 A094633 * A094635 A094636 A094637

Keyword

nonn

Author

Jonathan Sondow and Robert G. Wilson v, May 17 2004

Status

approved