%I #8 Dec 15 2017 17:35:39
%S 7,61,293,919,2693,9857,13763,16033,18367,39419,44789,64433,132511,
%T 157307,195407,213289,241513,243589,258331,293989,332573,436673,
%U 462067,478637,523777,583367,976933,983557,1329673,1481099,1582069,1753963
%N Numbers n such that p(n) + p(n+1) is a square and n is prime.
%e 7 is in the sequence because the seventh and eighth primes are 17 and 19. Added together, they make 36 which is the square of 6.
%t p = q = 1; Do[ p = q; q = Prime[n + 1]; If[ PrimeQ[n] && IntegerQ[ Sqrt[p + q]], Print[n]], {n, 1, 10^7/4} ]
%o (PARI) for(n=1,10^6,x=prime(n)+prime(n+1); if(issquare(x) && isprime(n),print(n)))
%K nonn
%O 1,1
%A _Jason Earls_, Sep 29 2001
%E More terms from _Robert G. Wilson v_, Sep 30 2001
|