login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A320840
Smallest N such that A092391(k) >= n for all k >= N.
2
0, 1, 1, 2, 3, 3, 5, 5, 6, 7, 9, 9, 10, 11, 11, 13, 13, 14, 17, 17, 18, 19, 19, 21, 21, 22, 23, 25, 25, 26, 27, 27, 29, 29, 33, 33, 34, 35, 35, 37, 37, 38, 39, 41, 41, 42, 43, 43, 45, 45, 46, 49, 49, 50, 51, 51, 53, 53, 54, 55, 57, 57, 58, 59, 59, 61, 65, 65
OFFSET
0,4
COMMENTS
For n >= 2, a(n) <= n - 1, and is exactly n - 1 for all n = 2^t + 2.
Consider the diverging sum Sum_{k>=0} 4^k/k!. For k >= a(n), v(4^k/k!, 2) = A092391(k) >= n. As a result, the sum contains only finitely many nonzero terms (and thus converges) modulo 2^n for all n, that is, it converges in the 2-adic field. Here v(k, 2) is the 2-adic valuation of k.
EXAMPLE
a(33) = 29 because A092391(28) = 31 < 33, A092391(29) = 33, A092391(30) = 34, A092391(31) = 36 and A092391(32) = 33. The smallest N such that A092391(k) >= 33 for all k >= N is N = 29.
MATHEMATICA
a[n_] := Module[{i = n-1-Boole[n >= 2]}, While[i+Total[IntegerDigits[i, 2]] >= n, i--]; i+1]; a[0]=0; Table[a[n], {n, 0, 67}] (* Jean-François Alcover, Nov 23 2018, from PARI *)
PROG
(PARI) a(n) = if(n, my(i=n-1-(n>=2)); while(i+hammingweight(i)>=n, i--); i+1, 0)
CROSSREFS
Sequence in context: A335599 A227065 A010761 * A161172 A093505 A238527
KEYWORD
nonn,base
AUTHOR
Jianing Song, Oct 22 2018
STATUS
approved