OFFSET
1,1
COMMENTS
Numbers of the form m*p^2, p prime and m squarefree (A005117). [Corrected by Peter Munn, Nov 17 2020]
The asymptotic density of this sequence is (6/Pi^2)*Sum_{n>=1} 1/prime(n)^2 = 0.274933... (A222056). - Amiram Eldar, Jul 07 2020
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
A046951(a(n)) = 2.
MATHEMATICA
Select[Range[2, 200], MemberQ[{2, 3}, (e = Sort[FactorInteger[#][[;; , 2]]])[[-1]]] && (Length[e] == 1 || e[[-2]] == 1) &] (* Amiram Eldar, Jul 07 2020 *)
PROG
(PARI) is(n)=my(f=vecsort(factor(n)[, 2], , 4)); #f && f[1]>1 && f[1]<4 && (#f==1 || f[2]==1) \\ Charles R Greathouse IV, Oct 16 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Apr 08 2003
STATUS
approved