OFFSET
1,1
COMMENTS
All the terms are divisible by 3 by definition.
Šalát (1994) proved that the asymptotic density of this sequence is 0.287106... (A336069).
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
3 is a term since A007949(3) = 1 is a divisor of 3.
MATHEMATICA
Select[Range[200], Mod[#, 3] == 0 && Divisible[#, IntegerExponent[#, 3]] &]
PROG
(PARI) isok(m) = if (!(m%3), (m % valuation(m, 3)) == 0); \\ Michel Marcus, Jul 08 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 07 2020
STATUS
approved