A355038 -id:A355038 - cp's OEIS frontend

A036537 Numbers whose number of divisors is a power of 2.

1, 2, 3, 5, 6, 7, 8, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 24, 26, 27, 29, 30, 31, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 46, 47, 51, 53, 54, 55, 56, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 88, 89, 91, 93, 94, 95, 97, 101, 102

Offset: 1

Labos Elemer

nonn

Primes and A030513(d(x)=4) are subsets; d(16k+4) and d(16k+12) have the form 3Q, so x=16k+4 or 16k-4 numbers are missing.
A number m is a term if and only if all its divisors are infinitary, or A000005(m) = A037445(m). - Vladimir Shevelev, Feb 23 2017
All exponents in the prime number factorization of a(n) have the form 2^k-1, k >= 1. So it is an S-exponential sequence (see Shevelev link) with S={2^k-1}. Using Theorem 1, we obtain that a(n) ~ C*n, where C = Product((1-1/p)*(1 + Sum_{i>=1} 1/p^(2^i-1))). - Vladimir Shevelev Feb 27 2017
This constant is C = 0.687827... . - Peter J. C. Moses, Feb 27 2017
From Peter Munn, Jun 18 2022: (Start)
1 and numbers j*m^2, j squarefree, m >= 1, such that all prime divisors of m divide j, and m is in the sequence.
Equivalently, the nonempty set of numbers whose squarefree part (A007913) and squarefree kernel (A007947) are equal, and whose square part's square root (A000188) is in the set.
(End)

			383, 384, 385, 386 have 2, 16, 8, 4 divisors, respectively, so they are consecutive terms of this sequence.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Vladimir Shevelev, S-exponential numbers, Acta Arithm., 175(2016), 385-395.
Eric Weisstein's World of Mathematics, Square Part, Squarefree Part

A005117, A030513, A058891, A175496, A336591 are subsequences.
Complement of A162643; subsequence of A002035. - Reinhard Zumkeller, Jul 08 2009
Subsequence of A162644, A337533.
Cf. A000005, A000188, A007913, A007947, A036538, A352780.
The closure of the squarefree numbers under application of A355038(.) and lcm.

a036537 n = a036537_list !! (n-1)
a036537_list = filter ((== 1) . a209229 . a000005) [1..]
-- Reinhard Zumkeller, Nov 15 2012

bi[ x_ ] := 1-Sign[ N[ Log[ 2, x ], 5 ]-Floor[ N[ Log[ 2, x ], 5 ] ] ]; ld[ x_ ] := Length[ Divisors[ x ] ]; Flatten[ Position[ Table[ bi[ ld[ x ] ], {x, 1, m} ], 1 ] ]
Select[Range[110],IntegerQ[Log[2,DivisorSigma[0,#]]]&] (* Harvey P. Dale, Nov 20 2016 *)

is(n)=n=numdiv(n);n>>valuation(n,2)==1 \\ Charles R Greathouse IV, Mar 27 2013

isok(m) = issquarefree(m) || (omega(m) == omega(core(m)) && isok(core(m,1)[2])); \\ Peter Munn, Jun 18 2022

from itertools import count, islice
from sympy import factorint
def A036537_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:all(map(lambda m:not((k:=m+1)&-k)^k,factorint(n).values())),count(max(startvalue,1)))
A036537_list = list(islice(A036537_gen(),30)) # Chai Wah Wu, Jan 04 2023

A209229(A000005(a(n))) = 1. - Reinhard Zumkeller, Nov 15 2012
a(n) << n. - Charles R Greathouse IV, Feb 25 2017
m is in the sequence iff for k >= 0, A352780(m, k+1) | A352780(m, k)^2. - Peter Munn, Jun 18 2022

A335988 Cubefull exponentially odd numbers: numbers whose prime factorization contains only odd exponents that are larger than 1.

Original entry on oeis.org

1, 8, 27, 32, 125, 128, 216, 243, 343, 512, 864, 1000, 1331, 1944, 2048, 2187, 2197, 2744, 3125, 3375, 3456, 4000, 4913, 6859, 7776, 8192, 9261, 10648, 10976, 12167, 13824, 16000, 16807, 17496, 17576, 19683, 24389, 25000, 27000, 29791, 30375, 31104, 32768, 35937

Offset: 1

Amiram Eldar, Jul 03 2020

nonn

This sequence is a permutation of A355038.
This sequence is also a permutation of the exponentially odd numbers (A268335) multiplied by the square of their squarefree kernel (A007947).
a(n)/rad(a(n)) is a permutation of the squares.
a(n)/rad(a(n))^2 is a permutation of the exponentially odd numbers.

			8 = 2^3 is a term since the exponent of its prime factor 2 is 3 which is odd and larger than 1.

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

Intersection of A001694 and A268335.
Intersection of A036966 and A268335.
A355038 in ascending order.
A030078, A050997, A092759, A179665, A079395 and A138031 are subsequences.
Cf. A007947, A065487.

Join[{1}, Select[Range[10^5], AllTrue[Last /@ FactorInteger[#], #1 > 1 && OddQ[#1] &] &]]

from math import isqrt, prod
from sympy import factorint
def afind(N): # all terms up to limit N
    cands = (n**2*prod(factorint(n**2)) for n in range(1, isqrt(N//2)+2))
    return sorted(c for c in cands if c <= N)
print(afind(4*10**4)) # Michael S. Branicky, Jun 16 2022

Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + 1/(p*(p^2-1))) = 1.2312911... (A065487).

A367406 The exponentially odd numbers (A268335) multiplied by their squarefree kernels (A007947).

Original entry on oeis.org

1, 4, 9, 25, 36, 49, 16, 100, 121, 169, 196, 225, 289, 361, 441, 484, 529, 144, 676, 81, 841, 900, 961, 64, 1089, 1156, 1225, 1369, 1444, 1521, 400, 1681, 1764, 1849, 2116, 2209, 2601, 2809, 324, 3025, 784, 3249, 3364, 3481, 3721, 3844, 4225, 4356, 4489, 4761

Offset: 1

Amiram Eldar, Nov 17 2023

Analogous to A355038, with the exponentially odd numbers instead of the square numbers (A000290).

This sequence is a permutation of the square numbers.

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

Cf. A000290, A007947, A053143, A064549, A065463, A268335, A355038, A367407, A367417, A367418.

s[n_] := n * Times @@ FactorInteger[n][[;;, 1]]; s /@ Select[Range[100], AllTrue[FactorInteger[#][[;; , 2]], OddQ] &]

b(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i,2]%2, f[i,1]^(f[i,2]+1), 0));}
lista(kmax) = {my(b1); for(k = 1, kmax, b1 = b(k); if(b1 > 0, print1(b1, ", ")));}

a(n) = A064549(A268335(n)).
a(n) = A268335(n)*A367417(n).
a(n) = A367407(n)^2.
a(n) = A268335(n)^2/A367418(n).
Sum_{k=1..n} a(k) = c * n^3 / 3, where c = (Pi^2/(15*d^3)) * Product_{p prime} (1 - 1/(p^3*(p+1))) = 1.78385074227198915372..., and d = A065463 is the asymptotic density of the exponentially odd numbers.
a(n) = A053143(A268335(n)). - Amiram Eldar, Nov 30 2023

cp's OEIS Frontend

A036537 Numbers whose number of divisors is a power of 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Python

Formula

A335988 Cubefull exponentially odd numbers: numbers whose prime factorization contains only odd exponents that are larger than 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A367406 The exponentially odd numbers (A268335) multiplied by their squarefree kernels (A007947).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula