A109560 - OEIS

%I #2 Mar 30 2012 17:34:19

%S 89,139,181,241,467,509,619,773,839,863,953,1409,1831,1847,1933,2113,

%T 2221,2251,2593

%N a(n) = If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 6, Prime[n]]

%t a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 6, Prime[n], {}], {n, 1, 400}]]

%Y Cf. A001359.

%K nonn

%O 0,1

%A _Roger L. Bagula_, Jun 27 2005