A160054 -id:A160054 - cp's OEIS frontend

A271050 Positive integer k such that k^2 = p^2 + q^2 - 1 where p and q are consecutive primes.

13, 17, 37, 157, 307, 437, 1451, 6487, 60773, 421133, 1445957, 2064493, 15247789, 177075397, 16509255853, 4270704255979, 56635565799013, 124750634536736711, 628179811369719907, 81815870181890275241, 13861061749008806276269, 91566796731172246399571

Offset: 1

Emre APARI, Mar 29 2016

nonn

Prime terms of this sequence are listed in A167276. - Altug Alkan, Mar 30 2016

			         7^2 + 11^2 - 1 = 169 (13^2, k is prime),
        11^2 + 13^2 - 1 = 289 (17^2, k is prime),
        23^2 + 29^2 - 1 = 1369 (37^2, k is prime),
      109^2 + 113^2 - 1 = 24649 (157^2, k is prime),
      211^2 + 223^2 - 1 = 94249 (307^2, k is prime),
      307^2 + 311^2 - 1 = 190969 (437^2, k is semiprime),
    1021^2 + 1031^2 - 1 = 2105401 (1451^2, k is prime),
  42967^2 + 42979^2 - 1 = 3693357529 (60773^2, k is prime).

Cf. A001248, A069484, A160054 (the corresponding primes p), A167276.

p = 2; q = 3; lst = {}; While[p < 10^15, If[ IntegerQ@ Sqrt[p^2 + q^2 - 1], AppendTo[lst, Sqrt[p^2 + q^2 - 1]];
Print[Sqrt[p^2 + q^2 - 1]]]; p = q; q = NextPrime@ q] (* Robert G. Wilson v, Mar 30 2016 *)

list(nn) = {p = 2; forprime(q=3, nn, if (issquare(s = q^2+p^2-1), print1(sqrtint(s), ", ")); p = q;);} \\ Michel Marcus, Mar 29 2016

More terms from Jinyuan Wang, Jan 09 2021

A359414 Primes prime(k) such that prime(k)^2 + prime(k+1)^2 - 1 is the square of a prime.

Original entry on oeis.org

7, 11, 23, 109, 211, 1021, 42967, 297779, 125211211, 11673806759

Offset: 1

Robert Israel, Dec 30 2022

Suggested in an email from J. M. Bergot.

There are no more terms < 10^100 unless the prime gap g = prime(k+1) - prime(k) exceeds 10000. For all known terms, g <= 14. There are no more terms < 10^1000 with g <= 14. - Jon E. Schoenfield, Dec 31 2022

			a(3) = 23 is a term because 23 is prime, the next prime is 29, and 23^2 + 29^2 - 1 = 37^2 where 37 is prime.

Subset of A160054.

R:= NULL: q:= 2:
while q < 2*10^8 do
  p:= q; q:= nextprime(q);
r:= p^2 + q^2 - 1;
  if issqr(r) and isprime(sqrt(r)) then R:= R, p fi
od:
R;

cp's OEIS Frontend

A271050 Positive integer k such that k^2 = p^2 + q^2 - 1 where p and q are consecutive primes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Extensions