OFFSET
0,2
EXAMPLE
n=9: a(9) = 41361, center = 1710732321, radius = 22; the nine primes in the zone are {1710732299, 1710732307, 1710732311, 1710732313, 1710732319, 1710732323, 1710732329, 1710732337, 1710732343}.
MATHEMATICA
f[n_] := (PrimePi[n^2 + Ceiling[ Log[n^2]]] - PrimePi[n^2 - Ceiling[ Log[n^2]] - 1]); t = Table[0, {15}]; Do[a = f[n]; If[a < 15 && t[[a + 1]] == 0, t[[a + 1]] = n], {n, 10^5}] (* Robert G. Wilson v, Jul 14 2004 *)
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Labos Elemer, Jul 14 2004
EXTENSIONS
Offset corrected and a(11) from Donovan Johnson, Jul 11 2010
STATUS
approved