A375039 The maximum exponent in the prime factorization of 2*n-1.

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

Offset

1,5

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Index entries for sequences computed from exponents in factorization of n.

Formula

a(n) = A051903(2*n-1).

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=2} (1 - 1/((1-1/2^k) * zeta(k))) = 1.25979668632898014495... .

Mathematica

a[n_] := Max[FactorInteger[2*n - 1][[;; , 2]]]; a[1] = 0; Array[a, 100]

Prog

(PARI) a(n) = if(n == 1, 0, vecmax(factor(2*n-1)[, 2]));

Crossrefs

Bisection of A051903.

Cf. A033150, A375040.

Sequence in context: A333355 A354523 A116531 * A101871 A101875 A081387

Adjacent sequences: A375036 A375037 A375038 * A375040 A375041 A375042

Keyword

nonn,easy

Author

Amiram Eldar, Jul 28 2024

Status

approved