OFFSET
1,1
COMMENTS
FORMULA
EXAMPLE
a(1) = 14 because 14 = 2 * 7 and 2 + 7 = 9 = 3^2 is semiprime.
a(2) = 21 because 21 = 3 * 7 and 3 + 7 = 10 = 2 * 5 is semiprime.
a(3) = 26 because 26 = 2 * 13 and 2 + 13 = 15 = 3 * 5 is semiprime.
a(4) = 28 because 28 = 2^2 * 7 and 2 + 7 = 9 = 3^2 is semiprime.
a(5) = 30 because 30 = 2 * 3 * 5 and 2 + 3 + 5 = 10 = 2 * 5 is semiprime.
a(6) = 33 because 33 = 3 * 11 and 3 + 11 = 14 = 2 * 7 is semiprime.
a(7) = 38 because 38 = 2 * 19 and 2 + 19 = 21 = 3 * 7 is semiprime.
MAPLE
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
MATHEMATICA
Select[Range[250], PrimeOmega[Total[Transpose[FactorInteger[#]][[1]]]]==2&] (* Harvey P. Dale, May 06 2013 *)
PROG
(PARI) is(n)=bigomega(vecsum(factor(n)[, 1]))==2 \\ Charles R Greathouse IV, Sep 14 2015
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, Feb 22 2006
EXTENSIONS
Corrected and extended by Emeric Deutsch, Mar 07 2006
STATUS
approved