|
MATHEMATICA
|
a = 0; n = 1; t = {}; Do[If[PrimeOmega[x^2 + x + 41]== 2, a = a + 1]; If[Mod[x, n] == 0, n = n*10; AppendTo[t, a]], {x, 1, 100000000}]; t
nn=8; With[{ep=If[PrimeOmega[#]==2, 1, 0]&/@Table[x^2+x+41, {x, 10^nn}]}, Table[ Total[Take[ep, 10^n]], {n, nn}]] (* Harvey P. Dale, Dec 12 2014 *)
|