OFFSET
1,2
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 1 by definition.
a(2) = 3 as 3 is prime and 1 + 3 = 4 = 2^2.
a(3) = 13 as 13 is prime and 3 + 13 = 16 = 4^2.
a(4) = 23 as 23 is prime and 13 + 23 = 36 = 6^2.
MATHEMATICA
lps[n_]:=Module[{p=NextPrime[n]}, While[!IntegerQ[Sqrt[n+p]], p= NextPrime[ p]]; p]; NestList[lps, 1, 40] (* Harvey P. Dale, Sep 07 2020 *)
PROG
(PARI) {print1(1, ", "); p=1; a=1; for(i=1, 10^4, p=nextprime(p+1); if(issquare(a+p), print1(p, ", "); a=p))}
CROSSREFS
Cf. A000040, A000290, A083016 Rearrangement of primes such that the sum of two consecutive terms is a square.
Cf. A062064 (case a(1)=2). - Zak Seidov, Oct 11 2014
KEYWORD
nonn,easy
AUTHOR
Jonathan Vos Post, Dec 27 2010
STATUS
approved