|
%I #13 Aug 03 2014 14:01:19
%S 5,7,11,23,47,83,107,263,347,467,503,863,887,1283,1487,1823,1907,2027,
%T 2063,2447,2903,3203,3623,4007,4127,4547,4703,4787,5387,5807,7523,
%U 7703,8147,8423,11423,11483,11807,12107,12227,12647,12983,13043,13163,14087,14207
%N Safe primes (A005385) (p and (p-1)/2 are primes) such that 6*p+1 is also prime.
%H David Consiglio, Jr. and T. D. Noe, <a href="/A075705/b075705.txt">Table of n, a(n) for n = 1..1000</a>
%e 47 is prime, so is (47-1)/2=23 and also 6*47+1=283; 83 is a prime, (83-1)/2=41 and 6*83+1=499, ...
%p ts_sg_var_pras := proc(nmax) local i,tren,atek; tren := 0: for i from 1 to nmax do atek := numtheory[safeprime](i): if (atek > tren) then if (isprime(atek)='true' and isprime(6*atek+1)='true') then tren := atek: fi; fi; od; end: seq(ts_sg_var_pras(i), i=1..3000);
%t Select[Range[20000], PrimeQ[#] && PrimeQ[(#-1)/2] && PrimeQ[6#+1] &] (* _T. D. Noe_, Nov 07 2011 *)
%t Select[Prime[Range[1700]],And@@PrimeQ[{(#-1)/2,6#+1}]&] (* _Harvey P. Dale_, Feb 28 2013 *)
%Y Cf. A005384, A005385, A059455, A007693.
%K nonn
%O 1,1
%A _Jani Melik_, Oct 02 2002
|
