%I #2 Mar 30 2012 17:34:19
%S 7,13,19,47,53,61,73,83,131,359,389,401,449,479,491,691,709,787,811,
%T 829,919,1021,1039,1051,1153,1171,1249,1399,1471,1627,1699,1723,1801,
%U 1879,2017,2029,2053,2069,2089,2099,2143,2297,2399,2447,2579,2621
%N a(n) = If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 4, Prime[n]]
%t a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 4, Prime[n], {}], {n, 1, 400}]]
%Y Cf. A001359.
%K nonn
%O 0,1
%A _Roger L. Bagula_, Jun 27 2005
|