OFFSET
1,1
COMMENTS
Also numbers greater than the square of their smallest prime-factor: a(n)>A020639(a(n))^2=A088377(a(n));
Numbers with at least 4 divisors. - Franklin T. Adams-Watters, Jul 28 2006
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = n + O(n/log n). - Charles R Greathouse IV, Sep 14 2015
EXAMPLE
8=2*2*2 and 10=2*5 are terms; 4=2*2 is not a term.
From Gus Wiseman, Jul 03 2019: (Start)
The sequence of terms together with their prime indices begins:
6: {1,2}
8: {1,1,1}
10: {1,3}
12: {1,1,2}
14: {1,4}
15: {2,3}
16: {1,1,1,1}
18: {1,2,2}
20: {1,1,3}
21: {2,4}
22: {1,5}
24: {1,1,1,2}
26: {1,6}
27: {2,2,2}
28: {1,1,4}
30: {1,2,3}
32: {1,1,1,1,1}
(End)
MATHEMATICA
Select[Range[100], PrimeNu[#]>1||PrimeOmega[#]>2&] (* Harvey P. Dale, Jul 23 2013 *)
PROG
(Haskell)
a080257 n = a080257_list !! (n-1)
a080257_list = m a024619_list a033942_list where
m xs'@(x:xs) ys'@(y:ys) | x < y = x : m xs ys'
| x == y = x : m xs ys
| x > y = y : m xs' ys
-- Reinhard Zumkeller, Apr 02 2012
(PARI) is(n)=omega(n)>1 || isprimepower(n)>2
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Feb 10 2003
EXTENSIONS
Definition clarified by Harvey P. Dale, Jul 23 2013
STATUS
approved