cp's OEIS Frontend

A114985 Numbers whose sum of distinct prime factors is semiprime.

%I A114985 #11 Sep 14 2015 10:36:01
%S A114985 14,21,26,28,30,33,38,46,52,56,57,60,62,63,69,70,74,76,85,90,92,93,94,
%T A114985 98,99,102,104,105,106,112,120,124,129,133,134,140,145,147,148,150,
%U A114985 152,166,171,174,177,178,180,182,184,188,189,190,195,196,204,205,207,208
%N A114985 Numbers whose sum of distinct prime factors is semiprime.
%C A114985 This is the semiprime analog of A114522 "numbers n such that sum of distinct prime divisors of n is prime." See also A110893 "numbers with a semiprime number of prime divisors (counted with multiplicity)."
%F A114985 {k such that A008472(k) is an element of A001358}. {k such that sopf(k) is an element of A001358}. {k = Product(Prime(j)^e_j) such that Sum(Prime(j)) is in A001358}.
%e A114985 a(1) = 14 because 14 = 2 * 7 and 2 + 7 = 9 = 3^2 is semiprime.
%e A114985 a(2) = 21 because 21 = 3 * 7 and 3 + 7 = 10 = 2 * 5 is semiprime.
%e A114985 a(3) = 26 because 26 = 2 * 13 and 2 + 13 = 15 = 3 * 5 is semiprime.
%e A114985 a(4) = 28 because 28 = 2^2 * 7 and 2 + 7 = 9 = 3^2 is semiprime.
%e A114985 a(5) = 30 because 30 = 2 * 3 * 5 and 2 + 3 + 5 = 10 = 2 * 5 is semiprime.
%e A114985 a(6) = 33 because 33 = 3 * 11 and 3 + 11 = 14 = 2 * 7 is semiprime.
%e A114985 a(7) = 38 because 38 = 2 * 19 and 2 + 19 = 21 = 3 * 7 is semiprime.
%p A114985 with(numtheory): a:=proc(n) local A,s,B: A:=factorset(n): s:=sum(A[j],j=1..nops(A)): B:=factorset(s): if nops(B)=2 and B[1]*B[2]=s or nops(B)=1 and B[1]^2=s then n else fi end: seq(a(n),n=2..250); # _Emeric Deutsch_, Mar 07 2006
%t A114985 Select[Range[250],PrimeOmega[Total[Transpose[FactorInteger[#]][[1]]]]==2&] (* _Harvey P. Dale_, May 06 2013 *)
%o A114985 (PARI) is(n)=bigomega(vecsum(factor(n)[,1]))==2 \\ _Charles R Greathouse IV_, Sep 14 2015
%Y A114985 Cf. A001358, A008472, A110893, A114522.
%K A114985 easy,nonn
%O A114985 1,1
%A A114985 _Jonathan Vos Post_, Feb 22 2006
%E A114985 Corrected and extended by _Emeric Deutsch_, Mar 07 2006