|
%I #7 Mar 19 2020 15:18:22
%S 54,66,90,156,240,270,306,474,570,576,636,750,780,1080,1320,1350,2034,
%T 2154,2406,2700,2760,3204,3240,3306,3480,3516,3756,3774,3984,4056,
%U 4086,4140,4146,4176,4716,4734,4794,5154,5370,5424,5550,5664,5700,5850,5856
%N Numbers m such that m^2+1 is prime and m^2-7 = prevprime(m^2) (= A007917(m^2)).
%C The corresponding primes are listed in A157183 resp. A157935.
%H Harvey P. Dale, <a href="/A157934/b157934.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Range[6000],PrimeQ[#^2+1]&&#^2-7==NextPrime[#^2,-1]&] (* _Harvey P. Dale_, Mar 19 2020 *)
%o (PARI) forstep(m=2,9999,2, isprime(m^2+1) & precprime(m^2)==m^2-7 & print1(m,","))
%Y Cf. A157183, A028883, A028882, A005574, A157935.
%K nonn
%O 1,1
%A _M. F. Hasler_, Mar 18 2009
|
