A195086 Numbers k such that (number of prime factors of k counted with multiplicity) less (number of distinct prime factors of k) = 2.

8, 24, 27, 36, 40, 54, 56, 88, 100, 104, 120, 125, 135, 136, 152, 168, 180, 184, 189, 196, 225, 232, 248, 250, 252, 264, 270, 280, 296, 297, 300, 312, 328, 343, 344, 351, 375, 376, 378, 396, 408, 424, 440, 441, 450, 456, 459, 468, 472, 484, 488

Offset: 1

Author

Harvey P. Dale, Sep 08 2011

Keywords

• nonn

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=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.
• Index entries for sequences computed from exponents in factorization of n

Crossrefs

Subsequence of A048108.

Subsequences: A030078 and A085986.

Cf. A001221, A001222, A025487, A057521, A060687, A195069, A195087, A195088, A195089, A195090, A195091, A195092, A195093, A046660, A257851, A261256, A264959.

Programs

Haskell

a195086 n = a195086_list !! (n-1)
a195086_list = filter ((== 2) . a046660) [1..]
-- Reinhard Zumkeller, Nov 29 2015

Mathematica

Select[Range[500],PrimeOmega[#]-PrimeNu[#]==2&]

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

Formula

A001222(a(n)) - A001221(a(n)) = 2.

A046660(a(n)) = 2. - Reinhard Zumkeller, Nov 29 2015