OFFSET
1,1
COMMENTS
Note that 2*3*5*7*11*13 = 30030.
EXAMPLE
5945941+5945940/30030 = 5946139 and 5945941-5945940/30030 = 5945743, which are primes.
MATHEMATICA
w=30030; s=""; For[i=1, i<10^3*2, p=Prime[i]; If[PrimeQ[p-((p-1)/w)]&&PrimeQ[p+((p-1)/w)], (*Print[p, ":", p-((p-1)/w), ", ", p+((p-1)/w)]; *)s=s<>ToString[p]<>", "]; i++ ]; Print[s]
Select[Prime[Range[10550000]], PrimeQ[# + (# - 1) / 30030] && PrimeQ[# - (# - 1) / 30030] &] (* Vincenzo Librandi, Jun 15 2013 *)
PROG
(Magma) [p: p in PrimesInInterval(5, 120000000)| IsPrime((30031*p-1) div 30030 ) and IsPrime((30029*p+1) div 30030)]; // Vincenzo Librandi, Jun 15 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Apr 27 2008, corrected May 10 2008
EXTENSIONS
More terms from Vincenzo Librandi, Jun 15 2013
STATUS
approved