OFFSET
1,1
LINKS
Ivan Neretin, Table of n, a(n) for n = 1..10000
Robert Baillie, Summing the curious series of Kempner and Irwin, arXiv:0806.4410 [math.CA], 2008-2015. See p. 18 for Mathematica code irwinSums.m.
FORMULA
Sum_{n>=1} 1/a(n) = 1.386312271262110321181505974797071257205562524228381227122302929089588534920... (calculated using Baillie's irwinSums.m, see Links). - Amiram Eldar, Feb 14 2022
MATHEMATICA
Select[Range[4000], DigitCount[#, 2, 1] == 10 &] (* Amiram Eldar, Feb 14 2022 *)
PROG
(PARI) isok(n) = hammingweight(n) == 10; \\ Michel Marcus, Feb 29 2016
(Python)
from itertools import islice
def A038461_gen(): # generator of terms
yield (n:=1023)
while True: yield (n:=n^((a:=-n&n+1)|(a>>1)) if n&1 else ((n&~(b:=n+(a:=n&-n)))>>a.bit_length())^b)
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
EXTENSIONS
Offset changed to 1 by Ivan Neretin, Feb 28 2016
STATUS
approved
