Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #20 Sep 08 2022 08:45:33
%S 991,2143,3331,4591,4663,5563,5653,8623,9181,9343,9631,12781,13033,
%T 13411,14851,15643,16183,16363,20161,20521,20593,23311,24373,25471,
%U 26641,27073,28711,30133,30763,32401,32491,34381,35803,37783,38431,38611,39043,39241
%N Prime numbers p such that p +- ((p-1)/3) are primes.
%H Vincenzo Librandi, <a href="/A137703/b137703.txt">Table of n, a(n) for n = 1..1000</a>
%e 991 is in the sequence because 991+(990/3) = 1321 and 991-(990/3) = 661, which are primes;
%e 2143 is in the sequence because 2143+(2142/3) = 2857 and 2143-(2142/3) = 1429, which are primes.
%t Select[Prime[Range[50000]], PrimeQ[# + (# - 1) / 3]&&PrimeQ[# - (# - 1) / 3]&] (* _Vincenzo Librandi_, Jun 15 2013 *)
%o (Magma) [p: p in PrimesInInterval(5, 40000)| IsPrime((4*p-1) div 3 ) and IsPrime((2*p+1) div 3)]; // _Vincenzo Librandi_, Jun 15 2013
%K nonn,easy
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Apr 27 2008
%E More terms from _Vincenzo Librandi_, Jun 15 2013