OFFSET
1,4
COMMENTS
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = i+j for 3-smooth numbers n = 2^i*3^j (A003586).
MATHEMATICA
smoothNumbers[p_, max_] := Module[{a, aa, k, pp, iter}, k = PrimePi[p]; aa = Array[a, k]; pp = Prime[Range[k]]; iter = Table[{a[j], 0, PowerExpand @ Log[pp[[j]], max/Times @@ (Take[pp, j-1]^Take[aa, j-1])]}, {j, 1, k}]; Table[Times @@ (pp^aa), Sequence @@ iter // Evaluate] // Flatten // Sort]; PrimeOmega /@ smoothNumbers[3, 10^5] (* Jean-François Alcover, Nov 11 2016 *)
PROG
(Haskell)
a069352 = a001222 . a003586 -- Reinhard Zumkeller, May 16 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Mar 18 2002
EXTENSIONS
Edited by N. J. A. Sloane, Oct 27 2008 at the suggestion of R. J. Mathar.
STATUS
approved