login
A094634
A094371(n) / A094372(n).
3
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.
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, 2, 40320}], L/Denominator[M] (revised by Eric W. Weisstein, May 17 2004, and by N. J. A. Sloane, Jul 03 2017)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved