A064504 Numbers k such that prime(k) + prime(k+1)*2 is a square.

9, 96, 461, 548, 565, 876, 1053, 1202, 1302, 1426, 1596, 2040, 3040, 4392, 5710, 6960, 7749, 7988, 8630, 8818, 9323, 9505, 10465, 11617, 11828, 12674, 13119, 17580, 18507, 18954, 22922, 23197, 25821, 33502, 35877, 37065, 37321

Offset

1,1

Links

Harry J. Smith, Table of n, a(n) for n = 1..175

Mathematica

okQ[list_]:=IntegerQ[Sqrt[First[list]+2Last[list]]];
PrimePi/@Transpose[Select[Partition[Prime[Range[40000]], 2, 1], okQ]][[1]]
(* Harvey P. Dale, Mar 05 2011 *)

Prog

(PARI) for(n=1, 10^6, if(issquare(prime(n)+prime(n+1)*2), print1(n, " ")))
(PARI) { n=0; default(primelimit, 21000000); for (m=1, 10^9, if (issquare(prime(m) + 2*prime(m + 1)), write("b064504.txt", n++, " ", m); if (n==175, break)) ) } \\ Harry J. Smith, Sep 16 2009

Crossrefs

Sequence in context: A065456 A052389 A197665 * A090448 A264219 A005545

Adjacent sequences: A064501 A064502 A064503 * A064505 A064506 A064507

Keyword

nonn,changed

Author

Jason Earls, Oct 06 2001

Extensions

Offset changed from 0 to 1 by Harry J. Smith, Sep 16 2009

Status

approved