A159803 Number of primes p with (2m+1)^2 - 2m <= p < (2m+1)^2.
1, 1, 2, 2, 1, 3, 2, 3, 4, 4, 3, 5, 3, 5, 4, 4, 5, 2, 6, 4, 4, 7, 3, 8, 5, 7, 6, 5, 7, 8, 10, 5, 8, 7, 10, 8, 7, 10, 9, 7, 10, 9, 13, 10, 11, 11, 11, 11, 11, 12, 9, 9, 11, 14, 12, 11, 12, 12, 11, 15, 12, 11, 14, 12, 12, 14, 15, 12, 15, 14, 17, 18, 20, 18, 17, 14, 18, 12, 15, 15, 15, 14, 21
Offset: 1
Keywords
Examples
1) m=1: 7 <= p < 9 => prime 7: a(1)=1. 2) m=2: 21 <= p < 25 => prime 23: a(2)=1. 3) m=3: 43 <= p < 49 => primes 43, 47: a(3)=2. 4) m=30: 3661 <= p < 3721 => primes 3671,3673,3677,3691,3697,3701,3709,3719: a(30)=8.
References
- L. E. Dickson, History of the Theory of Numbers, Vol, I: Divisibility and Primality, AMS Chelsea Publ., 1999
- R. K. Guy, Unsolved Problems in Number Theory (2nd ed.) New York: Springer-Verlag, 1994
- P. Ribenboim, The New Book of Prime Number Records. Springer. 1996
Programs
-
Maple
A159803 := proc(n) local a,p; a := 0 ; for p from 4*n^2+2*n+1 to 4*n^2+4*n do if isprime(p) then a := a+1 ; fi; od: a ; end: seq(A159803(n),n=1..120) ; # R. J. Mathar, Apr 22 2009
Extensions
More terms from R. J. Mathar, Apr 22 2009
Comments