|
%I #17 Feb 21 2023 09:58:53
%S 1,2,3,5,6,7,10,13,19,29,32,36,62,63,70,75,78,85,93,96,102,107,109,
%T 119,123,128,145,158,164,174,177,190,192,197,219,241,247,252,280,284,
%U 299,304,318,335,340,344,354,361,374,377,382,385,387,427,434,439,440
%N Numbers n such that 16n^2-2n-1 and 16n^2+2n-1 are both primes.
%H Harvey P. Dale, <a href="/A182048/b182048.txt">Table of n, a(n) for n = 1..3500</a>
%e a(1)=1 because 16*1^2-2*1-1=13 is prime and 16*1^2+2*1-1=17 is prime.
%t Select[Range[500], PrimeQ[16 #^2 - 2 # - 1] && PrimeQ[16 #^2 + 2 # - 1] &] (* _T. D. Noe_, Apr 16 2012 *)
%t Select[Range[500],AllTrue[16#^2-1+{2#,-2#},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Jan 01 2018 *)
%Y Cf. A171748, A089637.
%K nonn
%O 1,2
%A _Gerasimov Sergey_, Apr 08 2012
%E Corrected and extended by _T. D. Noe_, Apr 16 2012
|
