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.
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
KEYWORD
easy,nonn
AUTHOR
Enoch Haga, Oct 25 2008
STATUS
approved