|
|
A087910
|
|
Exponent of the greatest power of 2 dividing the numerator of 2^1/1 + 2^2/2 + 2^3/3 + ... + 2^n/n.
|
|
2
|
|
|
1, 2, 2, 5, 8, 5, 5, 13, 9, 10, 10, 12, 12, 12, 12, 22, 17, 18, 18, 21, 22, 21, 21, 27, 25, 26, 26, 27, 27, 27, 27, 40, 33, 34, 34, 37, 39, 37, 37, 48, 41, 42, 42, 44, 44, 44, 44, 54, 49, 50, 50, 53, 54, 53, 53, 58, 57, 59, 62, 58, 58, 58
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Problem 9 of the 2002 Sydney University Mathematical Society Problems competition asked for a proof that a(n) tends to infinity with n. While this is immediate from the theory of the 2-adic logarithm, elementary proofs are available.
a(n) tends to infinity with n implies that log(-1) = 0 in the 2-adic field, by setting x = 2 in -log(1-x) = Sum_{k>=1} x^k/k. - Jianing Song, Aug 05 2019
|
|
REFERENCES
|
A. M. Robert, A Course in p-adic Analysis, Springer, 2000; see p. 278.
|
|
LINKS
|
|
|
FORMULA
|
Sum_{k=1..n} 2^k/k = (2^n/n)*Sum_{k=0..n-1} 1/binomial(n-1,k), so a(n) >= n - v(n,2) - max_{k=0..n-1} v(binomial(n-1,k),2) = n - A007814(n) - A119387(n) = n - floor(log_2(n)), where v(n,2) is the 2-adic valuation of n. It seems that the equality holds if and only if n = 2^m - 1 for some m. - Jianing Song, Feb 22 2020
|
|
EXAMPLE
|
a(5) = 8 as 2^1/1 + 2^2/2 + 2^3/3 + 2^4/4 + 2^5/5 = 256/15 whose numerator is divisible by 2^8 but not by 2^9.
|
|
MAPLE
|
S:= 0:
for n from 1 to 100 do
S:= S + 2^n/n;
a[n]:= padic:-ordp(numer(S), 2);
od:
|
|
MATHEMATICA
|
s[n_] := -2^(n + 1) LerchPhi[2, 1, n + 1] - I Pi;
a[n_] := IntegerExponent[Numerator[Simplify[s[n]]], 2];
|
|
PROG
|
(PARI) a(n) = valuation(sum(k=1, n, 2^k/k), 2) \\ Jianing Song, Feb 22 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|