A064371 Zero, together with positive numbers k such that prime(k) + k is a square.

0, 5, 12, 86, 105, 176, 214, 230, 241, 412, 503, 696, 1065, 1147, 1170, 1273, 1334, 2021, 2455, 2600, 2660, 2772, 3299, 3332, 3365, 4417, 4861, 6478, 6572, 8115, 8858, 8905, 9229, 9380, 9590, 9692, 9749, 10501, 10829, 11338, 11633, 11690, 12099

Offset

1,2

Links

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

Mathematica

Join[{0}, Transpose[Select[Table[{n, Prime[n]}, {n, 15000}], IntegerQ[ Sqrt[ Total[#]]]&]][[1]]] (* Harvey P. Dale, May 26 2011 *)

Prog

(PARI) j=[]; for(n=0, 20000, if(n==0 || issquare(prime(n)+n), j=concat(j, n))); j
(PARI) { n=0; default(primelimit, 20000000); for (m=0, 10^9, if (m==0 || issquare(prime(m) + m), write("b064371.txt", n++, " ", m); if (n==300, break)) ) } \\ Harry J. Smith, Sep 14 2009
(Magma) [0] cat [n: n in [1..12500] | IsSquare(NthPrime(n)+n)]; // Bruno Berselli, May 26 2011

Crossrefs

Cf. A014688.

Sequence in context: A215866 A219288 A368074 * A249478 A009414 A009426

Adjacent sequences: A064368 A064369 A064370 * A064372 A064373 A064374

Keyword

nonn,changed

Author

Jason Earls, Sep 26 2001

Extensions

Edited by Harry J. Smith, Sep 14 2009

Status

approved