OFFSET
1,2
COMMENTS
a(n)<=2^n; see A057838 for the indices n where a(n)=2^n.
LINKS
Ray Chandler, Table of n, a(n) for n=1..4438
EXAMPLE
a(5) = 30 because it is the first integer which has five nonprime divisors (1, 6, 10, 15 and 30; the divisors 2, 3 and 5 are prime).
a(35) = 2^35 = 34359738368.
a(71) = 2^71 = 2361183241434822606848.
a(191) = 2^191 = 3138550867693340381917894711603833208051177722232017256448.
MATHEMATICA
a = Table[0, {100} ]; Do[ c = Count[ PrimeQ[ Divisors[ n ] ], False]; If[ c < 101 && a[[ c ]] == 0, a[[ c ]] = n], {n, 2, 10077696} ];
Table[SelectFirst[Table[{n, Count[Divisors[n], _?(!PrimeQ[#]&)]}, {n, 10000}], #[[2]]==k&], {k, 34}][[;; , 1]] (* The program generates the first 34 terms of the sequence. *) (* Harvey P. Dale, Mar 04 2024 *)
PROG
(PARI) sme(n) = {k = 1; while (sumdiv(k, d, ! isprime(d)) != n, k++); k; } \\ Michel Marcus, Dec 13 2013
(Haskell)
a055079 n = head [x | x <- [1..], a033273 x == n]
-- Reinhard Zumkeller, Dec 16 2013
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Labos Elemer, Jun 13 2000
EXTENSIONS
More terms from Robert G. Wilson v, Nov 20 2000
Edited by Ray Chandler, Aug 12 2010
STATUS
approved