OFFSET
1,1
MATHEMATICA
lst={}; Do[If[Plus@@Last/@FactorInteger[Prime[n]^2+1]==2, AppendTo[lst, n]], {n, 0, 6!}]; lst
Select[Range[400], PrimeOmega[Prime[#]^2+1]==2&] (* Harvey P. Dale, Jun 14 2011 *)
PROG
(PARI) n=0; forprime(p=2, 1e4, n++; if(bigomega(p^2+1)==2, print1(n", "))) \\ Charles R Greathouse IV, Jun 14 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Mar 24 2009
STATUS
approved