OFFSET
1,3
COMMENTS
Include 2 times the primes (once for the real axis, once for the imaginary axis).
More precisely, a(n) includes all Gaussian primes (with the appropriate norms) on the first quadrant's bounding semi-axes. All such Gaussian primes occur in pairs {p, pi} (one real and one imaginary associate), where p is a classical prime of the form 4m + 3 (so p is in A002145) and p <= n. - Rick L. Shepherd, Jun 16 2017
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
MATHEMATICA
nn = 100; t = Select[Flatten[Table[a + b*I, {a, 0, nn}, {b, 0, nn}]], PrimeQ[#, GaussianIntegers -> True] &]; t2 = Table[0, {nn}]; Do[f = Ceiling[Abs[i]]; If[f <= nn, t2[[f]]++], {i, t}]; Accumulate[t2] (* T. D. Noe, Aug 19 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Olivier Gérard, Aug 17 2013
STATUS
approved