OFFSET
1,1
COMMENTS
Smallest prime p(0) such that the n-chain governed by recurrence p(i+1)=(p(i)^2 + 1)/2 are all primes. Equivalently, least prime p(0) that generates a sequence of n 2-prime triangles, where p(k) is the hypotenuse of the k-th triangle and the leg of the (k+1)-th triangle.
For n>2, the last digit of a(n) is 1 or 9. - Ya-Ping Lu, May 17 2025
REFERENCES
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 258.
LINKS
C. K. Caldwell, The Prime Glossary, Pythagorean triples
H. Dubner, Posting to Number Theory List
H. Dubner and T. Forbes, Prime Pythagorean Triangles
H. Dubner and T. Forbes, Journal of Integer Sequences, Vol. 4(2001) #01.2.3, Prime Pythagorean triangles
T. Forbes, Posting to Number Theory List
T. Forbes, Posting to Number Theory List
EXAMPLE
5 is a(1) because (5^2+1)/2 = 13 is prime, but (13^2+1)/2 = 85 is not.
PROG
(Python)
from sympy import isprime, nextprime; m = lambda x: (x*x+1)//2; p = 2; D = {}
while p < 2185103796349763249:
p = nextprime(p); q = m(p); n = 1
while isprime(q) and isprime(m(q)): n += 1; q = m(q)
if n not in D: D.update({n: p})
[print(k, end =', ') for key, k in sorted(D.items())] # Ya-Ping Lu, May 17 2025
CROSSREFS
KEYWORD
hard,more,nonn
AUTHOR
Lekraj Beedassy, Apr 26 2005
EXTENSIONS
a(1) added by T. D. Noe, Jan 29 2011
STATUS
approved
