OFFSET
1,1
COMMENTS
All the terms are even by definition.
If m is a term then m*(2*k+1) is a term for all k>=1.
Šalát (1994) proved that the asymptotic density of this sequence is 0.435611... (A336067).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Tibor Šalát, On the function a_p, p^a_p(n) || n (n > 1), Mathematica Slovaca, Vol. 44, No. 2 (1994), pp. 143-151.
EXAMPLE
2 is a term since A007814(2) = 1 is a divisor of 2.
MATHEMATICA
Select[Range[2, 150, 2], Divisible[#, IntegerExponent[#, 2]] &]
PROG
(PARI) isok(m) = if (!(m%2), (m % valuation(m, 2)) == 0); \\ Michel Marcus, Jul 08 2020
(Python)
from itertools import count, islice
def A336066_gen(startvalue=2): # generator of terms >= startvalue
return filter(lambda n:n%(~n&n-1).bit_length()==0, count(max(startvalue+startvalue&1, 2), 2))
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 07 2020
STATUS
approved