%I #2 Mar 31 2012 12:38:14
%S 37,41,43,47,53,59,67,71,73,79,83,89,97,101,103,107,109,113,127,131,
%T 137,139,149,151,157,163,167,173,179,181,191,193,197,199,211,223,227,
%U 229,233,239,241,251,257,263,269,271,277,281,283,293,307,311
%N Prime numbers p such that p +- a, p +- b, p +- c and p +- d are prime numbers, where a, b, c and d are distinct positive integers with p > d, d > c, c > b and b > a.
%t s=""; q=1; For[i=2, i<10^2, p=Prime[i]; For[a=2, a<p, If[PrimeQ[p-a]&&PrimeQ[p+a], For[b=a+2, b<p, If[PrimeQ[p-b]&&PrimeQ[p+b], For[c=b+2, c<p, If[PrimeQ[p-c]&&PrimeQ[p+c], For[d=c+2, d<p, If[p>q&&PrimeQ[p-d]&&PrimeQ[p+d], (*Print[p, ":", a, ", ", b, ", ", c, ", ", d]; *)s=s<>ToString[p]<>", "; q=p]; d=d+2]]; c=c+2]]; b=b+2]]; a=a+2]; i++ ]; Print[s]
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Apr 27 2008; corrected May 02 2008
|