The OEIS is supported by the many generous donors to the OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A177891 Numbers n such that sum of proper (or aliquot) divisors of n is a semiprime. 1
 6, 9, 14, 15, 16, 18, 20, 22, 25, 33, 36, 38, 45, 46, 51, 52, 62, 68, 70, 72, 75, 80, 86, 87, 91, 93, 95, 99, 104, 105, 110, 116, 117, 118, 119, 130, 136, 141, 142, 143, 144, 145, 148, 154, 158, 159, 160, 162, 165, 166, 169, 183, 195, 196, 200 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS This is to A037020 as semiprimes A001358 are to primes A000040. The first four values are themselves semiprime. Contains k^2 if k is in A005383. - Robert Israel, Feb 16 2020 LINKS Robert Israel, Table of n, a(n) for n = 1..10000 FORMULA A001065(a(n)) is in A001358. EXAMPLE a(2) = 9 because the aliquot divisors of 9 are 1 and 3, whose sum is 4 = 2*2, semiprime. a(5) = 16 because the aliquot divisors of 16 are 1, 2, 4, and 8, whose sum is 15 = 3*5, semiprime. MAPLE filter:= proc(n) uses numtheory; bigomega(sigma(n)-n) = 2 end proc: select(filter, [\$1..1000]); # Robert Israel, Feb 16 2020 MATHEMATICA semiPrimeQ[x_] := Plus @@ Last /@ FactorInteger@ x == 2; fQ[n_] := semiPrimeQ[ DivisorSigma[1, n] - n]; Select[ Range@ 200, fQ] PROG (PARI) isok(n) = bigomega(sigma(n)-n) == 2; \\ Michel Marcus, Apr 05 2015 CROSSREFS Cf. A001065, A001358, A005383, A037020. Sequence in context: A232553 A155705 A267369 * A228804 A315975 A315976 Adjacent sequences: A177888 A177889 A177890 * A177892 A177893 A177894 KEYWORD nonn,easy AUTHOR Jonathan Vos Post, Dec 14 2010 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified November 30 21:14 EST 2023. Contains 367462 sequences. (Running on oeis4.)