OFFSET
1,1
COMMENTS
MATHEMATICA
T = Range[10^4]*0; q = # (3 # - 1)/2 &@ Range[81]; Do[ s = q[[i]] + q[[j]]; If[s <= 10^4, T[[s]]++], {i, Length[q]}, {j, i}]; Select[ Flatten@ Position[T, 2], PrimeOmega[#] == 2 &] (* Giovanni Resta, Jun 13 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Aug 06 2006
EXTENSIONS
a(14)-a(44) from Giovanni Resta, Jun 13 2016
STATUS
approved