login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A036538 Number of integers m <= 2^n such that d(m) = 2^k for some k = 0, 1, 2, 3, ... 2
2, 3, 7, 12, 23, 45, 89, 178, 356, 707, 1409, 2822, 5639, 11273, 22546, 45088, 90165, 180315, 360637, 721258, 1442491, 2884973, 5769941, 11539858, 23079721, 46159395, 92318705, 184637321, 369274467, 738548867, 1477097749, 2954195452, 5908390605, 11816780739 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

a(n+1)/a(n) is very close to 2; a(n)/2^n is near 0.7.

As n goes to infinity, lim a(n)/2^n = 0.687827... (see comments in A036537). - Vladimir Shevelev, Feb 28 2017

LINKS

Table of n, a(n) for n=1..34.

FORMULA

a(n) = number of 1s in f(tau(k)) mapped across k = 1..2^n, with f(x):= 1-sign(log_2 x - floor( log_2 x )). - Michael De Vlieger, Mar 01 2017

EXAMPLE

Of the numbers 1 .. 2^4 = 16, only 4, 9, 12 and 16 are not in A036537, so a(4) = 16 - 4 = 12.

MAPLE

IversonBrackets := expr -> subs(true=1, false=0, expr):

A := proc(n) option remember; {seq(2^k, k=0..n)} end:

h := proc(n) option remember; add(evalb(numtheory:-tau(j) in A(n)), j=2^(n-1) + 1..2^n); IversonBrackets(%) end:

a := n -> 1 + add(h(k), k=1..n); seq(a(n), n=1..17); # Peter Luschny, May 14 2018

MATHEMATICA

Table[Count[#, 1] &@ Table[1 - Sign[# - Floor@ #] &@Log[2, #] &@ DivisorSigma[0, x], {x, 1, 2^m}], {m, 1, 20}] (* original program edited by Michael De Vlieger, Mar 01 2017, or *)

1 + Accumulate@ Table[Count[Range[2^(n - 1) + 1, 2^n], k_ /; IntegerQ@ Log2@ DivisorSigma[0, k]], {n, 20}] (* Michael De Vlieger, Feb 28 2017 *)

PROG

(PARI) a(n) = sum(k=1, 2^n, d = numdiv(k); (d<=2) || (ispower(d, , &p) && (p==2))); \\ Michel Marcus, May 14 2018

CROSSREFS

Cf. A000005, A036537, A036539, A037992.

Sequence in context: A054176 A289977 A000630 * A108742 A018240 A090596

Adjacent sequences:  A036535 A036536 A036537 * A036539 A036540 A036541

KEYWORD

nonn

AUTHOR

Labos Elemer

EXTENSIONS

a(20)-a(26) from Michael De Vlieger, Feb 28 2017

a(27)-a(34) from Giovanni Resta, May 14 2018

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 17 08:36 EDT 2019. Contains 328107 sequences. (Running on oeis4.)