|
%I #11 Apr 15 2019 14:00:02
%S 5,241,503,3299,4861,9749,10501,11633,13799,14057,23003,40493,56569,
%T 74323,90023,90641,185177,224633,307187,367097,648887,732761,766679,
%U 954221,997103,1004567,1070423,1095491,1100311,1252997,1284083,1434757,1483397,1695139,1716521,1720151,2115823,2176817
%N Primes p such that p + prime(p) is a square.
%H Zak Seidov, <a href="/A242625/b242625.txt">Table of n, a(n) for n = 1..100</a>
%e p = 5, q = prime(5) = 11, p + q = 16 = 4^2.
%e p = 241, q = prime(241) = 1523, p + q = 1764 = 42^2.
%t Select[Prime[Range[200000]],IntegerQ[Sqrt[#+Prime[#]]]&] (* _Harvey P. Dale_, Apr 15 2019 *)
%o (PARI) s=[]; forprime(p=2, 3000000, if(issquare(p+prime(p)), s=concat(s, p))); s \\ _Colin Barker_, May 19 2014
%Y Cf. A073136.
%K nonn
%O 1,1
%A _Zak Seidov_, May 19 2014
|
