OFFSET
1,1
COMMENTS
Conjecture: There are infinitely many numbers m such that the two polynomials m^2 + 1 and m^2 - m + 1 are prime.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Paul T. Bateman and Roger A. Horn, A heuristic asymptotic formula concerning the distribution of prime numbers, Mathematics of Computation, Vol. 16, No. 79 (1962), pp. 363-367.
EXAMPLE
6 is in the sequence because 6^2+1 = 37 and 37-6 = 31 is prime.
10 is not in the sequence because 10^2+1 = 101 is prime but 101-10 = 91 = 7*13 is not prime.
MAPLE
nn:=2800:
for m from 2 by 2 to nn do :
p:=m^2+1:if isprime(p) and isprime(p-m) then printf("%d, ", m):
else fi:
od:
MATHEMATICA
Select[Range[3000], And @@ PrimeQ[#^2 + {1-#, 1}] &] (* Amiram Eldar, May 15 2026 *)
PROG
(PARI) isok(m) = isprime(m^2+1) && isprime(m^2-m+1); \\ Michel Marcus, May 16 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, May 15 2026
STATUS
approved
