login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 60th year, we have over 367,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Other ways to Give
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A145898 Ending prime: where number of consecutive squares m^2 satisfy r = p + 4*m^2, prime. 2
37, 151, 29, 23, 293, 107, 263, 83, 113, 107, 113, 131, 1607, 197, 239, 233, 313, 311, 317, 353, 383, 401, 443, 461, 499, 523, 503, 617, 659, 677, 743, 773, 773, 887, 857, 863, 881, 887, 911, 953, 983, 1013, 1283, 1129, 1277, 1283, 1301, 1319, 1619, 1433 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Farideh Firoozbakht noticed the unusually high number of 19 primes ending in a(13)=1607 in A145741 where she specified an m sequence of 1..10. My m sequence is open. 1,2,3,...,m+1, all squared.
LINKS
EXAMPLE
a(1)=37 because when m is 3, the first prime is 5 and the ending prime is 37: r=1+4*1^1=5, prime; and r=1+4*2^2=17, prime; and r=1+4*3^2=37, prime (and the next value of r does not produce a prime).
PROG
(UBASIC)
10 'p464
20 N=1
30 A=3:S=sqrt(N)
40 B=N\A
50 if B*A=N then 100
60 A=A+2
70 if A<=S then 40
80 M=M+1:R=N+4*M^2:if R=prmdiv(R) and M<100 then print N; R; M:goto 80
90 if M>=1 then stop
100 M=0:N=N+2:goto 30
CROSSREFS
Sequence in context: A251106 A141968 A142656 * A096700 A158596 A142688
KEYWORD
easy,nonn
AUTHOR
Enoch Haga, Oct 25 2008
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified December 11 01:45 EST 2023. Contains 367717 sequences. (Running on oeis4.)