OFFSET
1,1
COMMENTS
The powers of 2 (A000079), except for 1, are all terms.
The product of any two terms (not necessarily distinct) is also a term.
This sequence is a disjoint union of {2} and the subsequences of numbers m of the form 2^k*o where o = A000265(m), the odd part of m, is a k-free number, for k >= 2. These subsequences include, for k = 2, numbers of the form 4*o where o is an odd squarefree number (A056911); for k = 3, numbers of the form 8*o where o is an odd cubefree number; etc.
The asymptotic density of this sequence is Sum_{k>=2} 1/(zeta(k)*2*(2^k-1)) = 0.222707226888193809... .
The asymptotic mean of the exponent of 2 in the prime factorization of the terms of this sequence is Sum_{k>=2} k/(zeta(k)*2*(2^k-1)) = 0.69116459363299090827... .
This sequence is a subsequence of A360015 and the asymptotic density of this sequence within A360015 is exactly 1/2.
Also even numbers whose multiset of prime factors has unique mode 2. - Gus Wiseman, Jul 10 2023
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = 2*A360015(n). - Gus Wiseman, Jul 10 2023
EXAMPLE
From Gus Wiseman, Jul 09 2023: (Start)
108 = 2*2*3*3*3 is missing because its mode is not 2.
180 = 2*2*3*3*5 is missing because 2 is not the unique mode.
120 = 2*2*2*3*5 is present because its unique mode is 2.
The terms together with their prime factorizations begin:
2 = 2
4 = 2*2
8 = 2*2*2
12 = 2*2*3
16 = 2*2*2*2
20 = 2*2*5
24 = 2*2*2*3
28 = 2*2*7
32 = 2*2*2*2*2
40 = 2*2*2*5
44 = 2*2*11
48 = 2*2*2*2*3
52 = 2*2*13
56 = 2*2*2*7
60 = 2*2*3*5
64 = 2*2*2*2*2*2
(End)
MATHEMATICA
q[n_] := Module[{e = IntegerExponent[n, 2], m}, m = n/2^e; (m == 1 && e > 0) || AllTrue[FactorInteger[m][[;; , 2]], # < e &]]; Select[Range[256], q]
PROG
(PARI) is(n) = {my(e = valuation(n, 2), m = n >> e); (m == 1 && e > 0) || (m > 1 && vecmax(factor(m)[, 2]) < e)};
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Jan 21 2023
STATUS
approved