OFFSET
1,1
COMMENTS
From Amiram Eldar, Nov 07 2020: (Start)
Numbers whose powerful part (A057521) is either a cube of a prime (A030078) or a square of a squarefree semiprime (A085986).
The asymptotic density of this sequence is (6/Pi^2) * (Sum_{p prime} 1/(p^2*(p+1)) + Sum_{p<q primes} 1/(p*(p+1)*q*(q+1))) = (1/zeta(2)) * (2*P(3) + Sum_{k>=4} (-1)^(k+1)*(k-1)*P(k) + (Sum_{k>=2} (-1)^k*P(k))^2))/2 = 0.0963023158..., where P is the prime zeta function. (End)
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Alfred Rényi, On the density of certain sequences of integers, Publications de l'Institut Mathématique, Vol. 8 (1955), pp. 157-162.
FORMULA
A046660(a(n)) = 2. - Reinhard Zumkeller, Nov 29 2015
MATHEMATICA
Select[Range[500], PrimeOmega[#]-PrimeNu[#]==2&]
PROG
(PARI) is(n)=bigomega(n)-omega(n)==2 \\ Charles R Greathouse IV, Sep 14 2015
(PARI) is(n)=my(f=factor(n)[, 2]); vecsum(f)==#f+2 \\ Charles R Greathouse IV, Aug 01 2016
(Haskell)
a195086 n = a195086_list !! (n-1)
a195086_list = filter ((== 2) . a046660) [1..]
-- Reinhard Zumkeller, Nov 29 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Harvey P. Dale, Sep 08 2011
STATUS
approved